UNCLASSIFIED 


AD  NUMBER 


AD859289 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release^  distribution 
unlimited 


EROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors;  Critical 
Technology;  APR  1969.  Other  requests  shall 
be  referred  to  Naval  Air  Systems  Command^ 
ATTN:  AIR-6022,  Washington,  DC  20360. 


AUTHORITY 


USNASC  Itr  dtd  26  Oct  1971 


THIS  PAGE  IS  UNCLASSIEIED 


AD  859289 


Report  No.  MSD-P69-144 


u-1 

Contract  No.  00019-68-0-02^ 


COLLOCATION 

FLUTTER  ANALYSIS 
STUDY 


This  document  is  subject  to  special  export  controls  and  transmittal 
to  foreign  governments  or  foreign  nationals  may  be  made  only  with 
prior  approval  of  the  Naval  Air  Systems  Command  . 


VOLUME  IV. 

COFA  -  COMPUTER  PROGRAM  TO  PERFORM'^JLUTTER 
ANALYSIS  BY  THE  COLLOCATION  METHOD 


APRIL  1969 


0 


MISSILE  SYSTEMS  DIVISION 


I  HUGHES  I 

I _ j 

KUCHCS  A|INC<*AFT  company 


'.yd 


COFA 


COLLOCATION  FLUTTER  ANALYSIS  STUDY 


VOLUME  IV 


COFA  -  COMPUTER  PROGRAM  TO 
PERFORM  FLUTTER  ANALYSIS  BY  THE  COLLOCATION  METHOD 


Prepared  by  Dynamics  &  Environments  Section  Personnel 

Hughes  Aircraft  Company,  Missile  Systems  Division 
Contract  No.  00019-68-L-0247 


APRIL  1969 


This  document  is  subject  to  special  export  controls  and  transmittal 
to  foreign  governments  or  foreign  nationals.- may  be  made  only  with 
prior  approval  of  the  Naval  Air  Systems  Conimand  " 


TABLE  OF  CONTENTS 


Section  — 

ABSTRACT . . 

1  introduction . 1 

2  SIGN  convention . . 

3  DERIVATION  OF  EQUATIONS . 3 

4  REFERENCES . 13 

5  DESCRIPTION  OF  PROGRAM  INPUT  ....  14 

6  description  OF  PROGRAM  output  ...  40 

7  EXAMPLE  PROBLEMS . . 

8  PROGRAM  LISTING . . 

9  FLOW  DIAGRAM  . 101 

10  . . 


i 

i 


I 

I 

i 

I 

i 


ABSTRACT 


A  collocation  solution  of  the  flutter  and  vibration  problems  for  a 
multiple  component  s/stem  is  presented.  The  formulation  utilizes 
structural,  aerodynamic,  and  inertial  characteristics  in  the  form  of 
matrices  of  structural  and  aerodynamic  influence  coefficients  and  a 
mass  matrix,  respectively,  for  each  component.  The  use  of  a  rigid- 
body  modal  matrix  permits  a  general  analysis  for  a  system  free  in 
space  with  vtp  to  six  rigid-body  degrees  of  freedom. 

The  computer  program  provides  the  flutter  or  vibration  solution 
for  a  system  composed  of  as  many  as  20  flexible  components  with  a 
-'xximum  total  of  49  collocation  control  points.  An  option  is  provided 
to  vary  the  density  as  well  as  the  reduced  velocity.  Another  option  is 
provided  to  yield  the  modes  from  a  vibration  analysis  in  a  punched - 
card  format  for  use  in  flutter  analysis  by  modal  methods,  ■  >( 


SECTION  1 


INTRODUCTION 

The  mathematical  formulation  of  the  flutter  problem  results  in  a 
set  of  integral  equations  whose  closed  form  solution  is  impossible  to 
obtain  for  most  practical  problems.  One  of  the  most  useful  approximate 
methods  of  solving  these  equations  is  by  direct  collocation.  A  solution 
by  collocation  is  one  in  which  the  equations  are  satisfied  at  a  finite  num> 
ber  of  selected  points  on  the  structure.  These  points,  known  as  collo¬ 
cation  points,  are  satisfied  simultaneously.  The  collocation  solution 
results  in  a  matrix  formulation  which  when  cast  in  the  canonical  form 
will  yield  eigenvalues  that  are  directly  related  to  the  flutter  stability 
parameters.  This  manual  presents  a  digital  computer  program  that 
will  perform  collocation  flutter  analyses.  The  computer  program, 
which  is  written  in  Fortran  IV,  was  developed  by  Uodden,  Farkas,  and 
Malcpm  in  Reference  !. 

The  collocation  formulation  of  the  flutter  problem  has  been  pre¬ 
sented  in  Reference  2.  The  equations  are  presented  for  analysis  of 
single  component  systems  restrained  (cantilevered)  in  space  and  for 
sysmmetrical  systems  free  in  space  undergoing  either  symmetric  or 
antisymmetric  flutter.  A  method  of  generalizing  the  matrix  equation 
for  free -free  flutter  to  include  up  to  six  rigid -body  degrees  of  freedom 
has  been  given  in  Reference  3.  The  present  program  extends  the  form¬ 
ulation  of  Reference  2  to  include  an  arbitrary  combination  of  rigid -body 
degrees  of  freedom  (Ref;  3),  and  to  consider  more  than  one  flexible 
component.  In  addition,  an  option  has  been  provided  to  vary  the  altitude 
(i.  e.  density)  as  well  as  the  reduced  velocity.  Finally,  options  have 
been  added  to  darry  out  a  vibration  analysis  (which  requires  no  aerody¬ 
namic  data)  and  tf/  provide  vibration  modes  in  punched -card  format  for 
use  in  a  modal  flutter  or  vibration  analysis. 
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Th«  NASA  body  axis  system  with  the  x,  y,  and  k  axes  directed 
forward,  starboard,  and  downward,  respectively,  is  recommended  for 
consistency  with  the  formulation  of  the  static  aeroelastic  problems  in 
Reference  4.  However,  the  usual  flutter  convention  with  the  x,  y  and  z 
axes  directed  aft,  starboard  and  downward,  respectively,  may  be  used 
instead.  In  either  case,  the  vertical  normal  force  and  deflection  are 
positive  downward. 
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SECTION  3 

1 

*  DERIVATION  OF  EQUATIONS 

The  integral  equations  of  aeroelasticity  consist  of  tv^o  basic 
roilationships:  The  first  is  the  relation  between  the  structural  deforma¬ 
tion,  the  structural  influence  function,  and  the  inertial  and  aerodynamic 
loadings:  the  second  is  the  relation  between  the  aerodynamic  disturbance 
(downwash),  the  aerodynamic  influence  (kernel)  function,  and  the  aero¬ 
dynamic  pressure.  A  collocation  formulation  of  the  deformation  inte¬ 
gral  equation  for  a  vehicle  free  in  space  may  be  written  in  matrix  form 
by  requiring  that  the  integral  equation  be  satisfied  at  a  discrete  set  of 
control  points.  We  choose  a  single  type  of  coordinate,  viz. ,  the  deflec¬ 
tion  h,  as  an  adequate  measure  of  both  the  deformation  and  the  frec- 
stream  disturbance,  not  only  for  simplicity  in  the  resulting  equations 
but  also  because  deflections  have  a  more  general  meaning  on  a  cambered 
vehicle  and  deflection  influence  coefficients  are  more  readily  obtained 
from  a  structural  analysis  than  slope  (or  twist)  influence  coefficients. 
The  resulting  deformation  matrix  equation  is 

{Hj}  -  {hj  =  Ktaj  ({Fj}  +  {F^})  (1) 

where  {h|}  is  the  set  of  components  of  the  absolute  deflections  of  the 
control  points, •■{  h^  }  is  the  set  of  components  of  the  deflections  of  the 
control  points  due  to  the  rigid -body  motion  of  some  reference  points, 

[a]  is  the  set  of  structural  influence  coefficients  (SICs,  or  flexibility 
matrix)  for  the  s/stem  cantilevered  from  (or  otherwise  restrained  at) 
the  reference  point,  {  F.  }  is  the  set  of  inertial  force  components  inte¬ 
grated  throughout  the  region  adjacent  to  each  control  point,  {F^^}  is  the 
set  of  aerodynamic  force  components  integrated  over  the  vehicle  surface 
adjacent  to  each  control  point,  and  the  scaler  K  has  been  introduced  as 
c.  factor  to  the  SICs  for  convenience  in  investigating  variations  in  stiff¬ 
ness  levels.  The  inertial  forces  may  be  written  in  terms  of  a  mass 
matrix  [  M]  and  the  control  point  accelerations. 
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(2) 


t 


{Fj}  r  .(1/386)[M]  {hj} 

whefe  the  diagonal  elements  of  the  mass  matrix  are  found  from  inte¬ 
grating  the  structural  mass  density  throughout  the  region  adjacent  to  the  - 
control  points.  (N.  B.  •  the  mass  matrix  need  not  be  diagonal,  and,  in 
general,  will  not  be  so  if  the  elements  must  be  derived  from  a  set  of 
weight  data  previously  lumped  at  a  system  of  control  points  different  from 
those  required  in  the  aeroelastic  analysis.  The  use  of  a  coupled  mass 
matrix  permits  simulation  of  given  inertial  characteristics  at  a  set  of 
control  points  frequently  dictated  by  more  difficult  aerodynamic  consid¬ 
erations.  ) 

A  collocation  formulation  of  the  aerodynamic  integral  equation  is 
more  difficult  than  in  the  case  of  the  deformation  integral  equation  be¬ 
cause  of  the  singularities  in  the  aerodynamic  kernel  function.  The  de¬ 
termination  of  three  relationships  is  necessary  to  derive  a  set  of  aero¬ 
dynamic  influence  coefficients  (AICs)  that  relate  the  control  point  forces 
to  the  deflections.  The  most  bas^  and  difficult  is  the  pres  sure -downwash 
relation  that  is  derived  from  numerical  analysis  of  the  aerodynamic  inte¬ 
gral  equation.  The  simpler  relations  are  the  numerical  integration  of 
thr  pressure  to  obtain  the  force,  and  the  numerical  substantial  differen¬ 
tiation  of  the  deflection  to  obtain  the  downwash.  The  effort  involved  in 
each  step  depends  or  the  planform,  Mach  number  regime,  and  frequency 
range;  a  survey  of  the  status;  of  unsteady  AICs  is  given  in  Ref.  4.  For 
present  purposes,  it  is  sufficient  to  state  the  definition  of  the  AICs  in  the 
oscillatory  case.  We  write  the  aer  odynamic  control  point  forces  in 
terms  of  the  control  point  deflections  as 


{F^}  =  (4ir^/12)p)^bj8(W][Cj^l{hj}  (3) 


where  t  Cj^  ]  is  the  theoretically  derived  dimensionless  (complex)  matrix 
of  oscillatory  AICs,  f  is  the  frequency  of  the  assumed  harmonic  motion, 
p  is  the  atmospheric  density,  b^  is  the  reference  semichord,  s  is  the 
reference  span,  and  (  W  ]  is  an  empirically  derived  weighting  matrix 
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for  modification  of  the  theoretical  AICs.  A  method  for  obtaining  the 
elements  of  the  weighting  matrix  has  been  suggested  in  Ref»  4. 

The  sum  of  the  force  components  may  be  written  now  from  Eqs, 
(2)  and  (3)  for  the  case  of  harmonic  motion. 


{F.}  +  {Fj  =  (4ii-^f^/386)([M]  +  32.174  pb^8(W](Cj^l){iij>  (4a) 

=  (4Tr^£^/386)[M]{hj}  (4b) 

We  next  discuss  the  manner  of  inclusion  of  -  he  rigid-body  degrees 
of  jfreedom  in  Eq.  (1).‘  The  matrix  {h^}  has  been  defini  d  as  the  set  of 
components  of  the  deflections  of  the  control  points  due  ;o  the  rigid -body 
motion  of  the  reference  point.  Each  component  of  the  control  point  de¬ 
flections  h^  is  linearly  related  to  the  rigid-body  translations  and  rota¬ 
tions,  provided  the  rotations  are  small.  Therefore,  we  may  define  a 
rigid -body  modal  matrix  |h|^|  as  the  transformation 

{!>„}  =  [1>r]  Ur}  (5) 

where  is  the  set  of  amplitudes  of  rigid-body  translations  and  ro¬ 

tations  of  the  reference  point.>  As  an  example,  if  we  consider  symme¬ 
trical  vertical  motion,  j^hj^j’is  ci^.nposed  of  two  columns;  the  first  is 
a  unit  column  corresponding  to-the  pliuiging  mode,  the  second  consists 
of , the  x-coordinate  of  each  control  point  corresponding  to  the  pitching 
mode;  |a^|  is  composed  of  two  elements:  the  first  is  the  plunging 
displacement  z^,  the  Second  is  the  pitching  angular  displacement  9. 

Before  proceeding  in  the  derivation,  we  should  review  the  format 
of  the  various  matrices  in  the  case  of  a  multiple  flexible  component 
system.  As  an  example,  we  consider  a  symmetrical  flutter  analysis 
of  an  aircraft  whose  wing,  aft  fuselage,  and  tail  are  flexible,  and  whose 
forward  fuselage  may  be  assumed  to  be  rigid.  We  assume  that  the 
reference  point  (cantilever  point)  can  be  located  in  the  vehicle  such  that 
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,  i 


its  various  components  are  independent*  If  we  choose  a  point  at  the 
intersection  of  the  wing  and  fuselage,  then  the  wing  As  independeht  of 
the  aft  fuselage -tail  combination,  but  the  tail  and  aft  fuselage  niurt  be 
considered  together.  The  motion  of  the  rigid  forward  fuselage  is  deter¬ 
mined  by  the  motion  of  the  reference  point,  and  the  forward  fuselage 
will  not  enter  into  any  flexible  considerations  but  only  into  the  free-  * 
free  boundary  conditions.  From  the  foregoing,  it  is  seen  that'the  vari¬ 
ous  matrices  will  appear  in  partitioned  form.  If  we  denote  the  wing 
and  aft  fuselage -tail  system  by  the  subscripts  1  and  2,  respectively, 
then  the  flexibility  matrix  appears  as  . 


t.i  « 


[ 


Ttv 

t  ^ 


.  46) 


iT 


the  mass  matrix  as 


[M]  = 


'M,  ■ 

0 

U- 

7‘ 

M2 

(7); 


the  weighting  matri*  as 


[wl  = 


w 

11 

0 

(8) 
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the  AICh  uh 


[Chl  = 


Sli  “ 


and  the  rigid >body  modal  matrix  as 


{hpl 


(9) 


(10) 


Two  requirements  should  be  emphasized  with  regard  to  the  AICs.  The 
first  concerns  the  proper  inclusion  of  the  reference  geometry  associated  ' 
with  the  nondime  ns  ional  AICs.  The  dimensional  form  of  Kq.  (9)  may  be 
written 


brSfCj^l 


I 

0 

0 

b?s-C. - 
”2  2  h2 

(11) 


where  b^  and  s  are  the  reference  semichord  and  span  of  the  composite 
system,  b,^  and  s^  are  the  reference  geometry  for  the  first  component, 
and  b2  and  S2  are  the  reference  geometry  for  the  second  component. 
The  second  requirement  is  that  the  AICs  for  each  component  must  be 
determined  for  the  same  "dimensional"  reduced  velocity  V/«,  If  the 
reference  reduced  velocity  is 

1/k^  =  V/b^«  (12) 

the."!  the  reduced  velocity  for  the  first  component  must  be 
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(13) 


1/kj  =  (l/k^)(b^/bj) 
and,  for  the  second  component, 

l/k2  =  (l/k^)(b^/b2)  (14) 

Both  of  these  requirements  can  be  met  in  formulating  the  composite 
AICs  by  choosing  the  same  reference  geometry  in  determining  the  AlCs 
for  each  component. 

The  rigid-body  modal  m.atrix  provides  the  basis  for  a  general 
statement  of  the  boundary  conditions  for  the  free -free  flutter  of  the 
composite  system.  The  boundary  conditions  for  harmonic  motion  may 
be  written  as 


[hj^]^[M]  {hj}+[Am]  {a^^}  =  {0}  (15) 

where  (Am  ]  is  an  incremental  generalized  mass  matrix,  including  aero 
dynamic  effects,  of  any  rigid  component  of  the  system  attached  to  the 
reference  point  (e.  g. ,  the  forward  fuselage  that  was  assumed  to  be 
rigid  in  the  foregoing  example),*  and  is  not  considered  in  the  formula¬ 
tion  of  the  flexible  component  mass  and  aerodynamic  matrices.  The 
form  of  the  matrix  [Am]  may  be  illustrated  by  the  previous  example 
with  e  rigid  forward  fuselage  again  in  symmetrical  motion.  We  may 
wri 


[Am]  =  [Am]  +  [aQ] 


(I- 


N.  R.  :  It  is  assumed  that  no  dynamic  coupling  exists  between  the  rigid 
and  flexible  components.  A  suitable  distinction  can  always  be  made  be¬ 
tween  the  rigid  and  flexible  components  such  that  this  requirement  can 
be  met. 
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where  the  generalized  rigid  component  mass  matrix  of  the  forward 
fuselage  is 


[Am] 


(17) 


in  which  M  ,  S  ,  and  I  are  the  mass,  static  unbalance  about  the  ref- 
o  o  yo 

erence  point,  and  pitching  moment  of  inertia  about  the  reference  point, 
respectively,  of  the  forward  fuselage,  and  the  generalized  aerodynamic 
forces  on  the  forward  fuselage  (if  not  negligible)  are  found  from 


[AQl  =  32.174  pbSjh^jT[C^JIh„J 


Ro' 


'ho^ 


Ro' 


(18) 


where  [hj^^]  is  the  rigid-body  modal  matrix,  is  the  set  of  AICs, 

and  b^  and  are  the  reference  geometry  for  the  forward  fuselage. 
Again  the  AICs  must  be  found  for  the  reduced  velocity  of  the  composite 
system. 

We  are  now  in  a  position  to  eliminate  the  rigid -body  degrees  of 
freedom  and  to  formulate  the  eigenvalue  problem  for  the  flutter  of  the 
flexible  free-free  system.  Substituting  Eqs.  (4b)  and  (5)  into  Eq.  (1), 
and  adding  the  structural  damping  factor  1/(1  +  ig)  to  the  flexibility 
matrix  to  provide  the  artificial  structural  damping  necessary  to  sustain 
the  assumed  harmonic  motion  of  the  flutter  system,  we  obtain 


(I'll  ■  =  (4ir^Kf^/386(l  +  ig)  )[a]  []^  {hj}  (19) 

M{h^}  -  [h^]  {aj^})  =  [a]  [M]  (h^}  (20a) 

=  [Ul{hj}  (20b) 
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where  \  denotes  the  eigenvalue 


\  =  +  i\j  (21a) 

=  386(1  +  ig)/4iT^Kf^  (21b) 

Premultiplying  Eq.  (20b)  by  [m],  and  multiplying  Eq.  (15)  by  X 

and  subtracting,  permits  solution  for  the  amplitudes  of  the  rigid-body 
motion 

\{aj^}  =  -  [U]{hj}  (22) 

where 

=  [hj^]'^[M]thj^]  +[Aml  (23) 

Finally,  substituting  Eq.  (22)  into  Eq.  (20b)  yields  the  generalized 
matrix  equation  for  free -free  flutter 

Xihj}  =  ([I]  -  (hj^][J?il'^(hj^]'^(Ml)(U]  {hj}  (24) 

The  ')lution  of  Eq.  (24)  for  the  complex  eigenvalues  leads  to  the  free- 
fre  frequency  and  the  required  structural  damping.  From  Eq,  (21), 
we  obtain  the  frequency 

f  =  (l/2Tr)y386/K\j^  (25) 

and  the  required  structural  damping 

g  =  (26) 
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Since  the  formulation  of  the  AICs  requires  the  assumption  of  a  reduced 
velocity  l/k^,  the  velocity  follows  from  that  and  the  frequency  obtained 
in  Eq.  (25) 


U  =  0.5921  (2irfbj.)(l/ky)  (27) 

Equation  (24)  is  seen  to  be  completely  general,  being  applicable 
from  the  cantilever  case  (in  which  we  let  [h^^]  vanish)  to  the  case  of 
six  rigid  body  degrees  of  freedom,  and  for  a  vibration  analysis  for 
which  the  aerodynamic  terms  and  required  structural  damping  are  de¬ 
leted.  We  observe  that  Eq.  (24)  is  a  matrix  formulation  of  the  algebraic 
procedures  for  free-free  vibration  analysis  described  in  Ref.  4  (Par. 
11.2). 

From  a  series  of  solutions  of  Eq.  (24)  for  vario  is  reduced  velo¬ 
cities,  the  conventional  required  damping  versus  velocity  stability 
curve  can  be  constructed  for  a  specific  altitude,  and  the  flutter  pointis 
determined  as  the  velocity  for  which  the  required  damping  and  actual 
damping  are  equal.  An  alternative  approach  to  the  flutter  analysis  is 
based  on  a  single  representative  reduced  velocity  and  a  series  of  solu¬ 
tions  of  Eq.  (24),  carried  out  for  various  densities.  The  density  at 
which  the  required  damping  and  actual  damping  are  equal  may  be  used 
to  find  a  stiffness -altitude  similarity  parameter  for  flutter  from  which 
the  flutter  stability  may  be  determined.  However,  at  present,  the  va¬ 
lidity  of  this  latter  approach  requires  further  investigation,  particularly 
the  sensitivity  of  the  results  to  the  choice  of  representative  reduced 
velocity. 

The  generalized  mass  corresponding  to  each  free  vibration  mode 
is  of  interest  in  various  modal  analyses  of  flying  qualities,  stability 
and  control  characteristics,  or  transient  response  of  flexible  vehicles. 

If  Eq,  (24)  is  solved  for  the  free  vibration  modes  (by  deleting  the  aero¬ 
dynamic  terms)  then  the  n'th  generalized  mass  is  given  by 

m<")  =  {h5"^}'^[Ml{h^j"^}  +  {a^"h'^[Am]  {a'jJ'h  (28) 
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where  is  the  n'th  free  vibration  mod  and  the  corresponding  rigid 

component  mode  it  found  from  Eq,  (Z2) 


(.W,  .  .  (29) 

^  386  ^ 
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SECTION  5 


DESCRIPTION  OF  PROGRAM  INPUT 


UN  ns 


The  dimensional  data  required  for  each  component  consist  of  the 
mass  matrix  in  pounds,  the  flexibility  matrix  in  inches  per  pound,  and 
the  reference  semichord  and  semispan  in  feet.  The  aerodynamic  influ¬ 
ence  coefficients  are  dimensionless.  In  U.o  case  of  free-free  analysis, 
the  rigid-body  mass  matrix,  e.g. ,  in  a  symmetrical  analysis  if  and 
ly^  are  given  in  the  foot-pound  system,  the  x-coordinates  which  corre¬ 
spond  to  the  rigid  body  pitching  mode  must  be  measured  in  feet,  whereas 
if  and  1^^  are  given  in  the  inch-pound  system,  the  x-coordinate  must 
be  measured  in  inches.  The  density  is  required  in  slugs  per  cubic  foot. 


CLASSES  OF  DATA  AND  PROBLEMS 

Five  classes  of  data  must  be  provided:  mass,  aerodynamic  influ¬ 
ence  coefficients  (AICs)  and  their  weighting  matrices,  structural  influ¬ 
ence  coefficients,  the  rigid-body  motion  modal  matrix,  and  the  rigid 
component  generalized  mass  characteristics.  The  cantilever  case  does 
not  require  either  the  rigid  body  modal  matrix  or  the  generalized  masses. 
(For  a  vibration  analysis,  the  aerodynamic  input  is  not  required.) 

Several  classifications  of  problems  may  be  analyzed  using  the 
collocation  flutter  analysis  program.  They  are  cantilever  flutter 
analy'^is  —  the  structure  is  restrained  from  plunging  motion.  Free- 
free  .lUtter  analysis  the  structure  is  free  to  pitch,  plunge,  and  roll. 

The  T-free  cases  may  have  rigid  components  and  flexible  components. 
However,  when  flexible  components  are  coupled  together,  the  structural 
atta'  hment  between  component  must  be  statically  determinate.  When 
rigid  body  components  are  used,  any  number  up  to  six  rigid  body  modes 
may  be  used.  Also  vibration  analyses  may  be  performed  when  zero 
aerodynamic  forces  are  used. 
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PROGRAM  RESTRICTIONS  AND  OPTIONS 

1.  The  maximum  niunber  of  control  points  that  can  be  used  for 
all  flexible  components  of  any  system  is  49.  A  maximum  of  49  control 
points  may  also  be  used  for  the  rigid  component  for  the  purpose  of  de¬ 
riving  the  generalized  aerodynamic  force. 

2.  The  maximum  number  of  flexible  components  is  20. 

3.  The  maximum  number  of  values  used  in  the  reference  reduced 
velocity  (1/K|^)  series  is  20. 

4.  The  maximum  number  of  values  used  in  a  density  series  is  20. 

5.  The  program  provides  for  varying  the  densities  with  each 
(1/K|^)  or  for  using  the  same  densities  v/ith  all  (1/Kp)'s. 

6.  The  maximum  number  of  output  modes  is  25. 

7.  The  maximum  number  of  rigid-body  motion  modes  is  6. 

8.  It  is  possible  to  reserve  a  partition  in  the  upper  left-hand 
corner  of  the  flexible  components  AIC  matrices  for  control  points  whose 
aerodynamic  forces  may  be  neglected  or  found  from  an  alternate  theory 
to  that  used  for  the  primary  control  points.  This  partition  is  termed 
the  external  stores  region  since  external  stores  are  an  example  of  a 
source  of  additional  control  points  requiring  such  special  consideration. 
The  maximum  number  of  control  points  that  can  be  reserved  on  each 
flexible  component  for  external  stores  is  48. 

9.  A  weighting  matrix  is  an  optional  input  for  each  flexible  compo 
nent.  The  order  of  this  matrix  must  be  identical  with  the  order  of 

the  AIC  matrices  for  the  particular  component. 

10.  Any  number  of  complete  sets  (decks)  of  input  data  may  be 
stacked  and  run  in  one  ma’':hine  pass. 

DATA  DECK  SETUP 

Loading  Order 

The  data  decks  are  assembled  using  cards  punched  from  keypunch 
forms  and/or  card  that  are  punched-card  output  from  appropriate  com¬ 
puter  programs.  The  data  items  in  each  deck  have  the  following  order, 
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with  the  exception  that  some  of  the  items  may  be  omitted  if  indicators 
used  in  the  control  cards  specify  their  absence. 

I.  Title  card 

Data  deck  general  control  card 

3.  K  card  (flexibility  matrix  scale  factor) 

4.  Data  card  for  change  in  matrix  iteration  tests 

5.  Control  card(8)  for  external  stores  and  weighting  matrices 

6.  Reference  semichord  (b^)  and  reference  reduced  velocities 
(l/k^)'»  • 

7.  Reference  semichord  (b^.)  and  reference  semispan  (S.)  for 
■  rigid  and  flexible  components  (surfaces) 

8.  Density  series  (if  same  densities  are  used  for  all  (l/k^)'s 

9.  Generalized  mass  matrix  ([Am])  for  rigid  component** 

10.  Mass  matrix  ([M]^)  for  each  flexible  surface 

II.  Rig^'-body  motion  modal  matrix  ([hj^^])  for  rigid  component 

12.  Rigid-body  motion  modal  matrix  ([h^.])  for  each  flexible 
surface 

13.  Flexibility  matrix  ([a].) 

14.  Weighting  matrices  ([W]j) 

15.  Aerodynamic  input  repeated  for  each  (1/k^) 

a.  Density  series  cards  (if  densities  vary  with  each  (1/k^) 

b.  AIC  matrix  ([Cj.^])  for  rigid  component  (if  present) 

c.  AIC  matrix  ^or  each  flexible  surface 

Inpu  Data  Description 

/  1.  The  title  card  may  contain  any  aUphanumeric  characters 
desired  in  Columns  2  through  72. 


2.  Data  deck  general  control  card  (FORMAT  1814): 


Col 

1  -4 

5-8 

9-12 

13-16 

17-20 

21-24 

25-28 

29-32 

^33-36 

i— ^ 

Name 

NSUR 

NAERO 

NRIGID 

NFUS 

NDENS 

MODES 

NDELM 

NPUNCH 

NCOM 

•  •m} 

Field 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 
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NSUR: 

NAERO: 

NRIGID: 


NFUS: 


NDENS: 


MODES: 

NDELM: 


NPUNCH: 


NOON: 


Number  of  flexible  components  (surfaces),  £^20. 

Nv;mber  of  reference  reduced  velocities,  s20;  NAERO  = 

0  is  used  for  vibration  analyses. 

Number  of  rigid-body  motion  modes  to  be  input  (=  num¬ 
ber  of  columns  in  ([hj^]):  NRIGID  =  0  for  the  cantilever 
case. 

NFUS  must  =  1  if  AICs  j  =  1,  NAERO  are  input 

for  the  rigid  components;  NFUS  =  0  if  [C^^J  is  not  input. 

If  NDENS  -  0  the  densities  are  to  vary  with  each  1  / 
and  are  input  as  part  of  Item  15;  if  NDENS>0  this  number 
of  densities  must  be  input  as  Item  8,  and  each  density 
value  will  be  used  for  all  (1/k  )'s. 

Number  of  output  modes,  £23. 

NDELM  =  0  if  no  rigid  component  generalized  mass 
matrix  ([Am])is  input;  NDELM  =  1  if  [Am]  is  input. 

This  indicator  is  used  to  obtain  a  printout  of  the  dynamic 
matrix  ([U])  and  to  obtain  the  frequencies  and  modes  in 
punched-card  format;  NPUNCH  s  0  if  no  printout  of  [U] 
or  punched  output  is  desired;  NPUNCH  =  1  if  only  punched 
card  output  is  desired  and  NPUNCH  =  -1  will  provide  the 
printout  of  [U]  and  the  punched  output.  (The  minus  sign 
must  be  in  Column  29  and  the  1  (one)  in  Column  32.) 

This  indicator  provides  for  changing  five  program  test 
numbers  used  in  the  matrix  iteration  subroutine;  NCON  = 
0  if  no  changes  are  desired  and  NCON  0  if  any  of  the 
tests  is  to  be  changed. 


3.  K  card  (FORMAT  6E12.8) 


Col 

Name 

Field 


1-2 

K 

_ Qi _ 

JL 
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Field  1  contains  K,  the  flexibility  matrix  normalizing  constant;  if  the 
[a]  matrix  has  not  been  normalized,  enter  K  =  1.0.  The  flexibility 
matrix  calculated  by  the  program  FLUENC  has  not  been  normalized. 

4.  Datu.  card  for  changes  in  matrix  iteration  tests  (FORMAT 
3E12.8  and  2I4).  Omit  this  card  when  NCON  =  0.  There  are  three 
test  numbers  and  two  control  numbers  that  define  the  convergence  cri¬ 
teria  for  the  flutter  eigenvalue  solution.  A  suggested  set  of  numbers 
are  built  into  the  program;  these,  however,  may  be  changed  by  the 
program  user.  To  alter  any  number,  all  five  numbers  must  be  reentered. 


Col 

1-12 

13-24 

25-36 

37-40 

41-44 

- \ 

_ 

Name 

EPSP 

EPDP 

AITKEN 

NITRSP 

NITRDP 

■ 

Field 

(2) 

(3) 

(4) 

(5) 

■ 

0 

t  ^ 

*  * 

4  ' 

4  4 

*  ^ 

- ( 

* 

EPSP  =  0.5  X  lO"^  or  input  number;  test  for  eigenvector  conver¬ 
gence  when  the  iteration  procedure  is  approaching  a  single  root. 

EPDP  «  0,5  X  lO"^  or  input  number;  test  for  convergence  when  the 
interstlon  procedure  is  approaching  a  pair  of  close  roots. 

AITKEN  =  0.9  or  input  number;  if  this  test  is  met,  the  program 

2 

uses  a  procedure  (known  as  Aitken's  6  method)  to 
accelerate  convergence. 

NITRSP  =  40;  a  maximum  of  40  single  precision  arithmetic  itera¬ 
tions  will  be  performed  for  each  eigenvalue  if  its  eigen- 
vectors  have  not  converged  in  a  lesser  number. 

.  ITRDP  =  100;  if  the  eigenvectors  for  any  one  eigenvalue  have  not 
converged  in  NITRSP  single  precisions  iterations,  the 
program  will  then  use  up  to  a  maximum  of  100  double 
arithmetic  iterations. 

5.  Control  card(s)  for  external  stores  and  weighting  matrices 
(FORMAT  1814);  omit  this  data  when  NAERO  =  0. 
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Col 

Name 

Field 

1-4 

5-8 

9-12 

13-16 

- — < 

ISXTj 

BWj 

ISXT2 

1SW2 

1 

2 

3 

4 

»  m 

t  » 

*  ^ 

4 

Continue  on  successive  cards  to  ISWi  -  ISWnsuR- 

ISXT.  =  Number  of  control  points  reserved  for  the  external 
stores  on  surface  i. 

ISW^  =  0,  no  weighting  matrix  is  to  be  input  for  surface  i. 

=  1,  weighting  matrix  is  to  be  input  for  surface  i. 

Continue  on  next  card,  until  i  =  NSUR. 

6.  Reference  semichord,  b^  and  reduced  velocity  (1/k^)  series 
(FORMAT  6E12.8):  These  reference  parameters  are  used  in  computing 
the  flutter  velocities. 


Col 

1-12 

13-24 

25-36 

37-48 

49-60 

61-72 

Name 

B 

(i/Vi 

(l/kr)2 

(l/kr)3 

(1/V4 

(l/kp5 

Field 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

^  * 

* 

*  0 

0  L 

Continue  {l/k^).'s  on  next  card(s);  i^ZO.  The  b^(feet)  may  be  any  value 
but  it  is  noted  that  the  (1/k^)^  predetermines  the  (1/k^)  used  when  com¬ 
puting  the  AIC  matrices  for  each  surface.  The  1/k.  for  surface  1  is 

found  from  the  relationship  (1/k.)  =  (1/k  )(b  /b  .)  where  b  .  is  the 

1  r  r  n  ri 

reference  semichord  for  surface  i. 

7.  A  reference  semichord  (bj..)  and  semispan  (8|.)  must  be  input 
for  the  rigid  component  if  NFUS  =  1,  and  for  each  flexible  surface  if 
NAERO>0. 
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When  NFUS  =  1  (FORMAl  6E12.8) 


Column 

1-12 

1324 

25-36 

37-48 

49-60 

61-72 

Name 

•ro 

brj 

•ri 

brz 

®r2 

Field 

(1) 

(2) 

(3) 

(41 

(5) 

(6) 

1 

Where  and  are  the  reference  setnichord  and  semispan  fcr  the 
rigid  component.  Continue  on  next  card(s)  until  bf.  = 

•'i  ""  ■>^NSUR’ 


When  NFUS  =  0  (FORMAT  6E12.8) 


Column 

1-12 

13-24 

25-36 

37-48 

49-60 

61-72 

Name 

bri 

•ri 

•rz 

bra 

■r3 

Field 

(1) 

(2) 

(?) 

(4) 

(5) 

(6) 

0 

*  0 

’ 

•  ^ 

0  » 

0  0 

Continue  on  next  card(s)  until  b^j  =  ^rj^sUR  *^1  ~  *^NSUR‘ 

8.  Density  series  (FORMAT  6E12.  5):  Omit  this  input  if  NDENS=0 
in  Item  2.  If  NDENS>0  begin  in  field  1  of  this  card  and  NDENS  densities, 
5^20.  Continue  on  successive  card(8). 

C  -lumn 
Naju.' 

Field 


1-12 

13-24 

25-36 

37-48 

Pi 

P2 

P3 

P4 

1 

2 

3 

4 

9.  Rigid  component  ^generalized  mass  matrix  [Am j.The  {Am]  matrix 
for  the  rigid  component(s)  of  the  system  must  be  compatible  with  the  flex- 
ible  surfaces  product  matrix  given  by  [hj^J  [M][hj^],  i.e. ,  each  element 
in  [Am]  is  based  upon  the  same  rigid-body  motion  generalized  coordin¬ 
ate  as  the  respective  element  in  the  product  matrix.  Input  by  column 
beginning  each  column  on  a  new  card.  Orhit  this  data  when  NDELM  =  0. 


I 

r 


V. 

L 


% 

mr 


m 


iiie 


r* 


20 


Col 


Col 

Name 

Field 

1-12 

13-24 

25-36 

H  i — 

%  % 

^2,1 

^3.1 

•t  r*" 

^(NRIGID-1),  1 

^NRIGBC.  1 

(1) 

(2) 

(3) 

%  % 

(NRIGID-1) 

(NRIGID) 

■T  r— 

1 

Continue  on. successive  card(s)  until 


10.  Mass  matrix  [M]:  The  mass  matrix  is  partitioned  as  shown 
on  page  A-IC,  only  the  nonzero  partitions  are  input;  i.e.,  a  separate 
mass  matrix  ([M.])  is  input  for  each  surface.  The  sequence  for  con¬ 
sidering  the  surfaces  is  the  same  as  that  used  in  Item  5  and  7  if 
NAERO>0.  Repeat  the  following  input  for  each  surface,  from  i  =  1  to 
NSUR. 


Column 

Name 

Field 


Size  Control  Cards 


1-4 

5-8 

9-12 

- i 

NSIZEj 

(1) 

(2) 

(3) 

^  * 

0  <■ 

t  * 

i 

0 

NSIZE.  =  the  order  of  [M]. 

Ofte'  n.Any  of  the  elements  in  a  mass  matrix  are  zero;  the  following 
forma,  has  been  provided  so  that  nriost  of  the  zero  elements  will  not 
need  to  be  entered  into  the  program  as  data. 

Control  Card(s)  for  Omitting  ^eros  (FORMAT  1814) 


Column 

Name 

Field 

t 

1-4 

5-8 

9-12 

13-16 

17-20 

r  n 

LOWj 

LHIGHj 

LOW2 

LHIGH2 

LOW3 

IjlBI 

LHIGH. 

(1) 

(2) 

(3) 

(4) 

(5) 

r 

(2i) 

C 

_ > 

t  » 

*  t 

0  * 

0 

0 
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LOW.  =  The  row  number  in  which  the  first  nonzero  element 
1 

appears  in  Column  i. 

LHIGH^  =  The  row  number  in  which  the  last  nonzero  element 
appears  in  Column  i. 

If  only  one  nonzero  element  appears  in  Column  i  (i.  e. ,  a  diagonal  mass 
matrix)  the  row  number  in  which  it  appears  must  be  used  for  both  LOW 
and  LHIGH. 

Mass  Matrix  Elements  (FORMAT  6E12.8) 

Col 
Name 
Field 


The  elements  are  entered  by  column:  each  column  begins  on  a  new  card. 
Any  zero  elements  in  rows  between  LOW  and  LHIGH  must  be  entered 
or  their  respective  fields  left  blank.  If  external  stores  are  present 
(ISXT>0)  all  store  control  points  must  be  entered  before  the  surface 
control  points:  i.e.i  the  elements  representing  the  external  stcres 
mass  must  occupy  the  upper  left*hand  corner  of  the  mass  matrix. 

11.  Rigid  component  inodal  matrix,  lh|i^]  (see  page.  7  ).  Omit 
this  input  when  NFUS  =  0.  The  number  of  rows  in  must  be  the 

same  as  the  number  of  control  points  considered  when  computing  the 
[Cy^^  1  matrices:  the  number  of  columns  must  agrne  with  NRIGID. 


1-12 

13-24 

25-36 

^l.LOW 

iiml,LOV/+l 

^l,LOW+2 

^1,  LHIGH-l 

^1,  LHIGH 

(1) 

(2) 

(3) 

Size  Control  Card(s)  (FORMAT  1814) 

Column 

1-4 

5-8 

Name 

NROWS 

Field 

(1) 

(2) 

_  .  _  J 

m 

0 

L 


/, 

t 

«tis 


■’  * 


NROWS  =  The  number  of  rows  in  [hn^j 
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I'ASUij  i—H 


t 


i 


I 

I 

I 

I 

I 

1 

I 

I 

I 

I 

I 

I 

J 

I 

I 

I 

I 

I 

J 


,a  — 5-fe-  a.%#»  ' 


Matrix  Elements  (FORMAT  6E12.8) 


Column 

1-12 

13-24 

25-36 

37-48 

49-60 

61-72 

Name 

^^^cl,  1 

^Ro2,  1 

^^03,  I 

•  *  «  s 

s  •  •  • 

^RoNROWS,  1 

Field 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

* 

The  elements  are  entered  by  column,  with  each  column  beginning  on  a 

new  card.  _ 

12.  Rigid-body  modal  matrix,  [hj^]  (see  page  7).  Omit  this  in¬ 
put  if  NRIGID  =  0.  [hj^]  is  to  be  input  by  partitions  [hj^]j,  each  partition 
is  of  order  (NSIZE  x  NRKiID)  for  each  surface.  The  following  data  is 
to  be  repeated  for  i  =  1,NSUR. 


Size  Control  Card  (FORMAT  1814) 


Column 

Marne 

Field 


1-4 

5-8 

-  -  . 

NSIZE 

(1) 

(2) 

^  f  * 

0 

y 

* 

NSIZE  =  Number  of  control  points  on  each  surface 


Matrix  [hj^]j  Elements  (FORMAT  6E12.8) 


Column 

1-12 

13-24 

25-36 

37-48 

49-60 

61-72 

Name 

^Rl.l 

“Rj,! 

«  •  •  • 

•  •  «  • 

^^NSIZE,  1 

Field 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

' 

The  elements  are  entered  by  column,  with  each  column  beginning  on  a 
new  card. 

13.  Flexibility  matrices,  [a]  (see  page  6),  The  flexibility  matrix 
is  partitioned,  only  the  nonzero  partitions  [a]^  corresponding  to  the 
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flexible  surfaces  are  entered.  The  matrix  may  be  formed  by  any  of 
the  well-known  procedures  using  elementary  beam  theory,  force  or 
displacement  methods.  The  program  FLUENC  will  generate  this  matrix 
using  the  displacement  method.  The  punched  output  from  FLUENC  may 
be  used  as  direct  input  into  this  program.  The  following  data  is  re¬ 
peated  for  i  =  1,NSUR. 


Control  Card  (FORMAT  1814): 


Column 

Name 

Field 


1-4 

5-8 

9-12 

13-16 

- i 

m. 

X 

(BLANK) 

IFORM 

IROW 

(1) 

(3) 

(4) 

-  -  m 

»  $ 

S  0 

m^  =  The  number  of  rows  in  [aj^ 

IFC'nM  s  0  if  the  elements  are  to  be  input  using  column 
binary  format. 

s  1  if  the  elements  are  to  be  input  using  FORTRAN 
(FORMAT  6E12.8)  or  FLUENC  output  is  to  be  used 
directly. 

XROW  =  0  if  the  matrix  elements  are  to  be  entered  by 
column. 

=  1  if  the  matrix  elements  are  to  be  entered  by  row. 


Matrix  [a]^  elements  (use  format  specified  above): 

For  IFORM  =  1  and  IROW  =  1  (FORMAT  6E12.8) 


C'  mn 

1-12 

13-24 

25-36 

37-48 

- i 

Name 

^1,1 

^1.2 

^.3 

^1,4  . 

Fi>ld 

(1) 

(2) 

(3) 

(4) 

0 

1 

% _ 

0  0 

*  0 

0  ^ 

0 

Each  row  starts  on  a  new  card. 
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For  IFORM  =  1  and  |ROW  =  0  (FORMAT  6E12.8) 


Column 

1-12 

13-24 

25-36 

37-48 

- -I 

- 1 

Name 

^1.1 

"^2,1 

*3,1 

*4,1 

Field 

(1) 

(2) 

(3) 

(4) 

m 

*  * 

s  * 

s  ' 

[ 

Each  column  starts  on  a  new  card. 


If  IFORM  =  0,  then  IROW  must  =  0:  The  matrix  elements  are  input 
using  column  binary  format;  Column  1  starts  in  Origin  1.  Column  2 
starts  in  location  (1  +  m^);  Column  3  starts  in  location  (1  -I-  2m^):  etc. 

A  TRA’^card  must  end  each  [a^]  deck.  (The  column  binary  format  should 
be  used  only  if  the  data  are  available  as  punched-card  output  from  ap«* 
propriate  computer  programs.)  The  only  advantage  of  the  C-B  format 
is  the  minimum  card  storage  space  required. 

14.  Weighting  matrix,  [W]  (see  page  6),  The  weighting  matrix 
is  partit^v''ned,  only  the  nonzero  partitions  [W]|  corresponding  to  the 
flexible  are  entered.  No  provisions  have  been  made  for  enter¬ 

ing  a  iW,j].  matrix  iu’-  the  rigic'  component;  any  adjustment  to  [Cj,^] 
must  be  made  before  it  is  input  as  data.  If  ISW.  =  0  omit  this  data. 
Repeat  the  following  data  for  i  =  i,NSUR. 

For  {BW).  =  0  *nd  (ISXT)j>0 


Coniroi  card  fcx  external  stores  elements  (FORMAT  1814) 


Column 

Name 

Field 

1-4 

5  -8 

9- '.2 

13-16 

- i 

IQgllll 

1 

, BLANK) 

NFORM 

NROW 

(1) 

(?) 

(3) 

(4) 

I* 

*  ii 

’i^The  TRA  card  has  a  7  and9  punch  in  Column  1,  Column  2  through 
72  are  blank  and  Column  73  through  80  will  contain  the  characters  used 
for  identification  and  sei^uencing  in  the  punched  card  output  deck. 
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NXSTj  =  0  if  no  [W  ].  matrix  is  input  for  the  external  stores 
area  (the  program  will  use  a  unit  matrix,  I ) 

=  n  the  number  of  control  points  reserved  fer  stores* 
NFORM  =  1  if  the  [W]^  matrix  elements  will  be  input  using 
FORTRAN  (FORMAT  6E12.8} 

-  0  if  the  elements  are  to  be  input  using  column 
binary  format 

NROV/  =  0  if  the  [W^]  matrix  elements  are  to  be  input  by 
column 

=  1  if  the  matrix  elements  are  j  input  by  row 


External  stores  elements  [W]..  Format  given  on  control  card  above. 


For  NFORM  =  1  NROW  =  1 


Column 

1-12 

13-24 

25-36 

37-48 

- e 

Name 

*1,1 

"^1.2 

*1.4 

Field 

(1) 

(2) 

(3) 

(4) 

» 

*  0 

t  • 

0 

Each  row  starts  on  a  new  card. 


For  NFORM  *  I  NROW  =  0 


Co  mn 
NUii  T 
l*'ield 


1-12 

13-24 

25-36 

37-48 

- ^ 

*1.1 

'"2.1 

'"3.1 

*4.1 

(1) 

(2) 

(3) 

m 

0  0 

0  0 

0  0 

f 

Each  column  starts  on  a  new  card. 


I 


26 


/ 


I 


U  NFORM  s  0,  then  must  NROW  =  0:  The  matrix  elements  are  input 
using  column  binary  format:  Column  1  starts  in  Origin  1.  Column  2 
starts  in  location  (1  +  IXST^):  Column  3  starts  in  location  (1  4-  ZCCST^): 
etc.  A  TRA  card  must  end  each  [a^]  deck.  The  column  binary  format 
should  be  used  only  if  the  data  are  available  as  pimched-card  output 
from  appropriate  computer  programs.  The  only  advantage  of  C*B 
format  is  the  minimum  card  storage  space  required. 


Control  card  for  flexible  surface  weighting  matrix 


(FORMAT  1814)  The  [W^]  matrix  is  often  sparse,  sometimes 
diagonal  and  may  be  of  large  order,  ^49:  for  this  reason  we  provide 
tor  partitioning  of  the  matrix  and  entering  only  the  nonzero  partitions. 


Column 

Name 

Field 

4 

1-4 

5-8 

9-12 

13-16 

- ^ 

NSIZE, 

NPART 

NFORM 

NROW 

(1) 

(2) 

(3) 

(4) 

e  * 

.  ' 

'  0 

0  0 

NSIZE  =  The  number  of  control  points  used  on  surface  i. 

Do  not  include  control  points  for  external  stores. 
NPART  =  The  number  of  partitions  in  the  [Wjj  surface  matrix 
NFORM  =  1  if  the  [Wjj  will  be  input  using  FORTRAN  (FORMAT 
6E12.8) 

=  0  if  the  elements  are  to  be  input  using  column 
binary  format 

NROW  =  0  if  the  [W^]  matrix  elements  are  to  be  input  by 
column 

=  1  if  the  matrix  elements  are  to  be  input  by  row 


27 


Repeat  the  following  two  data  item  for  each  partition  j  =  l.NPART 


Control  card  for  partition  [W.Jj  FORMAT  (1814) 


Column 

1-4 

5-8 

9-12 

13-16 

17-21 

Name 

Field 

(1) 

(2) 

(3) 

(4) 

(5) 

=  The  order  of  partition  j  of  [Wj^] 


Elements  in  partition  [W.J..  Format  given  on  the  co!iL^‘o^  card  for 

1  J 

flexible  surface  weighting  matrix. 


NFORM  =  1  NRO\.  1 


Column 

1-12 

13-24 

25-3b' 

1 « 

37-48 

- - i 

Name. 

*1.1 

*1,2 

*1,3 

*1,4 

Field 

(i) 

(2) 

(3) 

(4) 

Each  row  starts  on  a  new  card. 


Nf  ORM  =  1  NROW  =  0 


C  u.nn 

1-12 

13-24 

25-36 

37-48 

- e 

Name 

*1,2 

*2.1 

'"3.1 

*4.1 

Field 

(1) 

(2) 

(3) 

(4) 

Each  column  starts  on  a  new  card. 


^  ^«=ri.'  •■"  ’5;“^ 


If  NFORM  --  0,  then  must  NROW  =  0:  The  matrix  elements  are  input 
using  column  binary  format;  Column  1  starts  in  Origin  1.  Column  2 
starts  in  location  (1  4-  DCST^);  Column  3  starts  in  location  (1  f  21XST£) 
etc.  A  TRA  card  must  end  each  [aj  deck.  The  column  binary  format 
should  be  used  only  if  the  data  are  available  as  punched-card  output 
from  appropriate  computer  programs.  The  only  advantage  of  C-B 
format  ie  the  minimum  card  storage  space  required. 

For  (BW).  =  0  and  (BXT).  =  0 

Control  card  for  flexible  surface  weighting  matrix  [W]^  (FORMAT  1814). 
The  (Wj  matrix  is  often  sparse,  sometimes  diagonal  and  may  be  of 
large  order,  ^49i  for  this  reason  we  provide  for  partitioning  of  the 
matrix  and  entering  only  the  nonzero  partitions.  Repeat  the  following 
data  for  i  =  1,NSUR. 


Column 

Name 

Field 


NSIZE  =  The  number  of  control  points  used  on  surface  i 
NPART  -  The  number  of  partitions  in  the  surface  matrix 

NFORM  ^  1  if  the  [w]j  will  be  input  using  FORTRAN  (FORMAT 
"  6E12.8) 


=  0  if  the  elements  are  to  be  input  using  column 
binary  format 

NROW  •  0  if  the  elements  are  to  be  input  by  column 
'  -  1  if  the  elements  are  to  be  input  by  row 


Repeat  the  following  two  data  items  for  each  partition  j  =  1,  NPART. 
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Control  card  lor  partition  [W^jj  FORMAT  (1814) 


Column 

Name 

Field 

4 

1-4 

5-8 

9-12 

13-16 

- \ 

N. 

J 

(1) 

(2) 

(3) 

(4) 

•  4 

4  * 

4  0 

P 

Nj  =  The  order  of  partition  j  of  [W.] 


Elements  in  partition  [W.]j  format  given  on  the  control  card  for  flexible 
surface  weighting  matrix. 


NFORM  =  1  NROW  =  1 


Oolumn 

1 - 

1-12 

13-24 

25-36 

37-48 

Name 

'"i.i 

'"l,2 

'^1,3 

'*'1,4 

Fielvl 

(i) 

(2) 

(3) 

(4) 

Each  row  : 

ctarts  or 

NFORM 

1  a  new  < 

=  1 

:ard. 

NROW  = 

0 

Column  j 

1-12 

13-24 

25-36 

37-48 

Na  Tie 

'"i,. 

^^2,1 

"'4,1 

Field 

(i) 

(2) 

(3) 

(4) 

Each  column  starts  on  a  new  card. 


If  NFORM  -•  0,  then  must  NROW  =  0;  The  matrix  elements  are  input 
using  column  binary  format;  Column  1  st-.rts  in  Origin  1.  Column  2 
starts  in  location  (1  +  m^);  Column  3  starts  in  location  (1  +  2m.):  etc. 

A  TRA  card  must  end  each  [a.]  deck.  The  column  binary  format  should 
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<^ed  only  if  the  data  are  available  as  punched-card  output  from  ap¬ 
propriate  computer  programs.  The  only  advantage  of  C-B  format  is 
the  minimvun  card  storage  space  required. 

15.  Aerodynamic  data:  The  aerodynamic  input  consists  of 
NAERO  sets  of  AIC's.  Each  set  of  AIC's  consists  of  the  AIC's  for 
each  surface  which  have  the  same  reference  1/k^  (see  Item  6).  If 
NDENS  =  0  a  density  series  will  precede  each  set  of  AIC's.  Input  the 
density  series  (if  Item  8  was  omitted)  and  the  AIC  matrices  for  each 
surface  with  the  following  input  order.  Repeat  the  order  for  each 

(1/k  ).  j  =  1.  NAERO. 

J 

Control  card  for  density  series  (FORMAT  1814).  Omit  this  input  if 
NDENS>0. 

Column 
Name 
Field 

NRHO  =  The  number  of  densities  to  be  entered  for 

(1/k  ).s20 
'  r'j 

Density  Series  (FORMAT  6E12.8) 

Column 
Name 
Field 

Rigid  component  AIC  matrix  Omit  this  input  if  NFUS  =  0.  The 

|,ChyJ-  matrix  may  be  sparse;  thus,  provision  has  been  made  to  parti- 
tion  the  matrix  and  enter  only  the  nonzero  partitions. 
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Reference  1/k  for  [Cjj  ].  (FORMAT  6E12.8) 


Column  I  1-12  13-24  25-36  37-48 


1/k 


(1)  (2)  (3)  (4) 


The  reduced  velocity  used  in  computing  the  rigid 
component  [Ch^Jj-  The  AIC's  in  [Ch^]  must  be 
computed  for  a  1/k.  which  properly  recites  them  to 
thejt^'l/k^. 


Field 


Control  card  for  matrix  (FORMAT  1814) 


1-4 

5-8 

9-12 

13-16 

NSIZE 

NPART 

NFORM 

NROW 

(1) 

(2) 

(3) 

(4) 

NSIZE  =  Number  of  control  points  on  the  rigid  component 
NPART  =  Number  of  nonzero  partitions  in 
=  1  for  an  unpartitioned  matrix 

NFORM  =  1  if  the  matrix  is  to  be  input  using  FORTRAN 

(FORMAT  6E12.8) 

=  0  if  the  [ChQ]j  matrix  is  to  be  input  using  column 
binary  format 

NROW  =  1  if  the  matrix  is  to  be  input  by  row 

=  0  if  the  matrix  is  to  be  input  by  column 

Repeat  the  following  data  for  each  partition,  K  =  1,  NPART. 


Partition  Size  Card  (FORMAT  1814) 


Column 

1-4 

5-8 

9-12 

13-16 

Name 

N 

Field 

(1) 

(2) 

(3) 

(4) 

1 

N  =  The  order  of  partition  k 


Elements  in  partition  K  of  Format  is  given  on  the  control  card 

for  matrix.  All  the  elements  in  the  AIC  matrices  are  complex 

numbers,  but  the  complexity  is  considered  in  the  program.  Thus,  each 
partition  may  be  input  as  though  it  is  a  real  matrix  of  size  Nx2N.  The 
real  elenients  form  the  odd  number  columns^,  and  the  imaginary  ele¬ 
ments  in  the  even  numbered  columns. 


For  NFORM  =  1  NROW  =  1 


Column 

1-12 

13-24  1 

1  25-36 

37-48 

- 1 

Name 

a(R.)|,, 

a(Re)j^  2 

\ 

Field 

(1) 

(2) 

(3) 

(4) 

0 

0 

0  * 

*  0 

0 

Each  row  starts  on  a  new  card. 


-For  NFORM  =.  1  and  NROW  =  0 


Coluniti 

1-12 

13-24 

25-36 

37-48 

49-60 

- e 

.  » 

Name 

a(Re)  j 

a(I)i^l 

a(Re)2^  1 

“<®2,  1 

Field 

(1) 

(2) 

(3) 

(4) 

(5) 

0 

L  ' 

0 

0  0 

0  0 

0 

Each  column  starts  on  a  new  card. 
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For  NFORM  =  0,  NROW  must  =  0.  Use  column  binary  format.  Column 
1  starts  in  card  Origin  1,  Column  2  in  location  (1  +  2N),  Column  3  in 
location  (1  +  4N),  etc.  A  TRA  card  must  end  each  deck.  The  column 
binary  format  should  be  used  only  if  the  data  are  available  as  punched- 
card  output  from  appropriate  computer  programs.  The  only  advantage 
of  C-B  format  is  the  minimum  card  storage  space  required. 

Flexible  component  AIC  matrix  The  AIC  matrices  are 

often  sparse:  thus,  a  provision  is  made  partitioning  the  matrix  and 
entering  only  the  nonzero  partitions.  The  following  data  is  repeated 
fori=l,NSUR. 

For  ]SXTj>0. 


Control  card  for  external  stores  partition  of  the  surface  i 


AIC  matrix  (FORMAT  1814) 


Column 

1-4 

5-8 

9-12 

13-16 

- i 

Name 

NXST 

(BLANK) 

NFORM 

NROW 

.  . i 

Field 

(1) 

(2) 

(3) 

(4) 

*  • 

4  * 

* 

NXST^  -  Number  of  control  points  reserved  for  external  stores 
NFORM  =  1  of  the  matrix  elements  are  to  be  input  using 
FORTRAN  (FORMAT  6E12.8) 

=  0  if  the  matrix  elements  are  to  be  entered  using 
column  binary  format 

NROW  =  1  if  the  matrix  elements  are  to  be  entered  by  row 

=  0  if  the  matrix  elements  are  to  be  entered  by  column 


Elements  for  external  stores  partition  of  AIC  matrix 
For  NFORM  =  1  and  NROW  =  1 


Column 

1-12 

13-24 

25-36 

37-48 

r - e 

Name 

a(Re)  j 

.(Re)j_2 

»(«1.2 

"  '  { 

Field 

(1) 

(2) 

(3) 

(4) 

4 

*  # 

f  4 

*  4 

4 

Each  row  begins  on  a  new  card. 
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For  NFORM  =  1  and  NROW  =  0 

Colunm 

Name 

Field 


If  NFORM  =  0  then  NROW  must  =  0  use  column  binary  format. 
Column  1  starts  in  card  origin  1,  Column  2  in  Location  (1+2NXT.), 
Column  3  in  Location  (l-f4NXST^),  etc.  A  TRA  Card  must  end  each 
deck.  The  column  binary  format  should  be  used  only  if  the  data  are 
available  as  punched^card  output  from  appropriate  computer  programs. 
The  only  advantage  of  C-B  format  is  the  minimum  card  storage  space 
requirements. 

Reference  1/k.  card  for  control  point  AIC  matrix 

FORMAT  (6E12.8) 

Column 

Name 

Field 

1  /k  .  -  The  reduced  velocity  used  in  computing  the  flexible 

component  |Cj^  y.  The  AIC's  must  be  computed  for  a 
1  /k^  which  properly  relates  them  to  the  jth  1  / k^. 
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Control  card  for  control  point  AIC  matrix  C,  ..  (FORMAT  (1814) 

n  xj 


Column 

1-4 

5-8 

9-12 

13-16 

- 1 

Name 

NSIZE 

NPART 

NFORM 

NROW 

% 

Field 

(1) 

(2) 

(3) 

(4) 

4 

#  * 

#  * 

.  J 

W 

* 

NSIZE  =  Order  of  control  point  AIC  matrix 
NPART  =  Number  of  partitions  in  |Cjj|ij 
NFORM  =  1  matrix  elements  are  to  be  input  using  FORTRAN 
(FORMAT  6E1Z.3) 


=  0  matrix  elements  are  to  be  input  using  column  binar/ 
format 

NROW  =  1  elements  input  hy  row. 

=  0  elements  input  by  column. 


Repeat  the  following  data  lor  j  s  2,  NPART.  j  =  1  corresponds  to  the 
external  stores  partition 


Control  card  for  partition  size  of  partition  j  FORMAT  (18J4) 


Column 

1-4 

5-8 

9-12 

Name 

N 

Field 

(1) 

(2) 

(3) 

N  =  C 

)rder 

of  parti 

tion  j 

t 


1. 

1 

I 


Elements  of  control  point  AIC  matrix  [Cjjjj  j  partition  j 


For  NFORM  =  1  and  NROW  =  1  (FORMAT  6E12.8) 


Column 

1-12 

13-24 

25-36 

37-48 

- e 

Name 

a(RE)l,  1 

a(I)l,l 

a(Re)l,  2 

a(I)1.2 

Field 

(1) 

(2) 

(3) 

(4) 

*  * 

, 

*  # 

Each  row  starts  on  a  new  card 


For  NFORM  =  1  and  NROW  =  0  (FORMAT  6E12.8) 


Column 

1-12 

13-24 

25-36 

37-48 

- 1 

Name 

a(Re)l,l 

a(I)l.  1 

a(Re)2, 1 

a(I)2, 1 

Field 

(1) 

(2) 

(3) 

(4) 

*  0 

.  ' 

t  » 

0 

Each  column  starts  on  a  new  card 

For  NFORM  =  0  then  NROW  must  =  0  use  column  binary  format. 
Column  1  starts  in  card  origin  1,  Column  2  in  Location  (1  +  2N), 

Column  3  in  Locatit  n  (1  +  4N),  etc.  A  TRA  card  must  end  each  deck. 
The  column  binary  format  should  be  used  only  if  the  data  are  available 
as  punched  card  output  from  appropriate  computer  programs.  The  only 
advantage  of  C-B  format  is  the  minimum  card  storage  space 
requirements 
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For  ISXT.  =  0 
1 


Reference  1/k.  card  for  control  point  AIC  matrix  C.  .. 

»  h  ij 


FORMAT  (6E12.8) 


Column 

1-12 

13-24 

25-36 

Name 

1/kj 

Field 

(1) 

(2) 

(3) 

l/k.  =  The  reduced  velocity  used  in  computing  the  flexible 

component  The  AIC's  must  be  computed  for  a 

1  /k  j  which  properly  relates  them  to  the  jth  1  /  k^. 


Control 


card  for  control  point  AIC  matrix 


FORMAT  (1814) 


Column 

1-4 

5-8 

9-12 

13-16 

- X 

Name 

NSIZE 

NPART 

NFORM 

NROW 

Field 

(1) 

(2) 

(3) 

(4) 

0 

^  0 

0  * 

0  0 

.  ^ 

0 

NSIZE 

NPART 

NFORM 


Order  of  control  point  AIC  matrix 
Numbejr  of  partitions  in  jc^^j 

1  matrix  elements  are  to  be  input  using  FORTRAN 


(FORMAT  6E12.8) 

=  0  mat.’ix  elements  are  to  be  input  using  column  binary 
format. 

NROW  -  1  elements  input  by  row 

=  0  elements  input  by  column 
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Repeat  the  following  data  for  j  =  1>  NPART. 


Control  card  for  partition  size  of  partition  j  FORMAT  (1814) 


Column 

1-4 

5-8 

9-12 

Name 

N 

Field 

(1) 

(2) 

(3) 

N  =  Order  of  partition 

Elements  of  control  point  AlC  matrix  j  partition  j 


For  NFORM  =  1  and  NROW  =  1  (FORMAT  6E12.8) 


Column 

1-12 

13-24 

mm 

37-48 

- i 

Name 

a(Re)l,l 

a(I)l.l 

a(Re)l,2 

»(I)1,2 

Field 

(1) 

(2) 

(3) 

(4)  . 

0 

0 

0i 

»  0 

*  0 

*  0 

Each  row  starts  on  a  new  card 


For  NFORM  =  1  and  NROW  =  0  (FORMAT  6E12.  8) 


Column 

1-12 

13-24 

25-36 

37-48 

- 1 

Name 

a(Re)l.  1 

a(I)l,l 

a(Re)2, 1 

a(I)2, 1 

Field 

(1) 

(2) 

(3) 

(4) 

0 

0  0 

.  ' 

0  0 

f  0 

0 

For  NFORM  =  0  then  NROW  must  =  0  use  column  binary  format. 
Column  1  start  in  card  origin  1,  Column  2  in  Location  (1  +  2N), 


Column  3  in  i^ocation  (1  +  4N),  etc.  A  TRA  card  must  end  each  deck. 
The  column  binary  format  should  be  used  only  if  the  data  are  available 
as  purched  card  output  from  appropriate  computer  programs.  The  only 
advantage  of  C-B  format  is  the  minimum  card  storage  requirements. 
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SECTION  6 


OESCIUPTION  OF  PROGRAM  OUTPUT 

Printed  Output 

1.  All  input  data. 

2.  The  dynamic  matrix  or  flutter  determinant  if  NPUNCH  =  -1 
is  used  in  the  general  control  card. 

3.  ;  For  the  vibration  analysis  or  lor  each  1/k^  in  a  flutter  analysis 

a.  The  eigenvalue  for  each  output  mode  followed  by  the 
number  of  single-  and  double -precision  iterations  and 
the  number  of  Aitken  accelerations. 

b.  The  normalized  eigenvectors  (modes)  followed  by  the 
check  eigenvalues  and  eigenvectors. 

c.  The  frequencies  (omegas)  in  cycles  per  second  followed 
(in  a  flutter  analysis)  by  the  structural  damping  coefficient 
and  the  velocity  (knots)  assoicated  with  each  frequency. 

d.  If  NPUNCH  a  dkl,  the  sequencing  numbers  used  for  identi¬ 
fying  the  punched -card  output  (frequencies  and  modes); 
this  will  conclude  the  printout  for  each  1/k^  used  in  a 
flutter  analysis. 

4.  In  a  vibration  analysis,  the  generalized  mass  corresponding 
to  each  output  (free  vibration)  mode  will  follow  the  above 
printed  output  [see  Eq.  (28),  Section  l). 

5.  A  number  of  different  statements  may  be  printed  which  indicate 
machine  or  program  detected  errors  in  input  data  (wrong  for¬ 
mat  or  incompatibility  in  the  number  of  rows  input  for  a  matrix 
and  the  number  designated). 

6.  If  the  program  or  the  machine  fails  in  the  iteration  portion  of 
the  program,  a  note  will  be  printed  stating  the  type  or  cause 
(jf  failure. 


AO 


I 

I 

I 

I 

I 


7.  If  convergence  is  not  obtained  in  the  allowable  number  of 
iterations,  a  note  will  be  printed  and  the  program  will  con> 
tinue.  (In  this  case  the  eigenvalues  and  eigenvectors  should 
be  compared  with  the  check  eigenvalues  and  eigenvectors  to 
determine  if  the  convergence  is  sufficiently  accurate,) 

6.  The  printed  output  for  the  example  problem  is  shown  on  the 
following  pages. 


! 
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SECTION  7 


EXAMPLE  PROBLEMS 

As  example  problems,  we  choose  a  cantilever  flutter  analysis  and 
a  free -free  symmetric  flutter  analysis  of  the  jet  transport  wing  (and 
rigid  fuselage)  treated  throughout  Ref.  5  and  shown  in  Figure  A-1.  The 
wing  mass,  flexibility  and  aerodynamic  matrices  and  the  fuselage  aero¬ 
dynamic  matrix  (symmetric  case)  can  be  seen  in  the  example  problem 

printed  output  (pages  41  -62)  and  will  not  be  repeated  with  the  list  of  in- 

-7 

put  data  to  follow.  The  flexibility  matrix  normalizing  factor  is  K  =  10  . 

The  wing  weighting  matrix  is  taken  as  unity  ([l])  and,  in  this  case,  re¬ 
quires  no  input.  The  symmetrical  rigid  component  (one -half  of  the 
fuselage)  generalized  mass  matrix  is  given  below  in  the  pound -inch 

system 

17,400  1,370,250 

1,370,250  4,457,907,200 

U 


25%  CHORD 


35%  ELASTIC  AXIS 


T/a» 


PITCH  AXIS- 


IDlSIRiaBagi^i 


\"n' 


75%  CHORD 


ROLL  AXIS 


Figure  A-1,  Jet  transport  wing  geometry. 
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In  the  antisymmetric  case,  the  rigid  component  generalized  mass  nnatrix 
would  be  [Am]  s  <tnd  in  the  composite  longitudinal -lateral  case, 

the  generalized  mass  matrix  would  be 


[Am] 


MS  0 
o  o 


The  symmetrical  case  requires  the  two  rigid-body  degrees  of  freedom 
of  plunging  and  pitching.  The  rigid-body  modal  matrix,  therefore,  con¬ 
sists  of  two  columns:  a  unit  column,  and  a  column  of  the  x-coordinate 
of  each  control  point.  The  rigid -body  modal  matrix  for  the  wing  is 


[hj^l  =  [1  x] 


1  20.25 

1  -81.00 
1  17.85 

1  -71.40 
1  1 5.  80 

1  -63.20 
1  1 3.  30 

1  -53.20 
1  11.05 

1  -44.20 


The  rigid* body  modal  matrix  for  the  fuselage  is' 


[1  x] 


*1  -373.  30' 
1  -248.  30 
1  -123.30 
1  +  1.70 

.1  +126.  70. 


Note  that  the  above  matrix  is  used  in  computing  the  incremental  gen¬ 
eralized  mass  which  results  from  the  aerodynamic  loads  on  one -half  of 
the  fuselage  and,  therefore,  the  x-coordinates  must  be  given  in  the 


43 


proper  order  for  the  control  points  used  in  computing  the  fuselage  AICs. 
[in  this  problem  the  fuselage  AICs  are  hypothetical,  but  the  x-coordi- 
nates  are  given  in  the  order  necessary  for  the  slender-body  theory  AICs.3 

The  reference  geometry  for  the  wing  is  =  5,  468  ft  and  s^ 

=  37.917  ft.  The  reference  geometry  for  the  fuselage  is  b^^  =  5.468  ft 
and  s^  =  18.  9585  ft.  [We  assume  that  the  wing  reference  geometry 
was  used  to  nondimensionalize  the  fuselage  AIC  matrix,  and,  because 
we  require  only  one -half  of  the  fuselage  aerodynamic  force,  it  can  be 
obtained  by  setting  s^  =  (1/2)8  .  ]  Both  example  problems  are  carried 
out  for  the  single  reduced  velocity  l/l'j,  ~  16.  67  with  the  reference 
semichord  for  the  system  b^  =  5.468  ft  and  with  sea -level  density 
p  =  0.  002378  slugs /ft^. 

The  output  modes,  the  flutter  frequencies  and  velocities,  and  the 
required  structural  damping  can  be  seen  in  the  example  problem  printed 
output  (pages  A- 55  through  A-65). 

The  input  keypunch  sheets  are  given,  followed  by  the  computer 
output. 


■BHDHRBEIHDHaHRVDHn 

peuHnHRanaiEiBnHnHnHR 

■BHEXHClHEIHiCBBHIDHaHE] 

■iBainainflinHiBainanHioan 

■RnairHiKjaiBHnHinHnHinHn 

mfa 

PK 

■■in 

■ciHiaHiia 

■mav^am 

■■aacui 

■IMHIMBi 

saHi 

I'M 

EaMEa 

i-IMIt-i 

nHFi 

■■i■i■i■i■i 

BNIiMiHH 

■■■i 

■■iU 

BS 

■■■■■I 

■■ 

■■iia 

■■■■I 

mm 

pH 

■■■■■■ 

■■ 

■■■ 

iHH 

■■■■«■ 

■■■■■■I 

■H 

1BBBC3 

HIHHIBiHI 

Hlii 

BIB^B 

■IHHIIH 

■■1 

■■m 

BBBSB 

SB 

■■iia 

icirii 

intai 

irnm 


ir^i 


ir;i 

iBiMMam 

iqsBmnnEsi 


IWgl— »11 

!13aif3Bini 

IMMEIHKJI 

irtBliilHBI 

nPKIMKI 


leai 


IKII 


IClHKdBI 

inBriBii 


IHBIL 
inrramgliftf! 
ipiaoiBfaiMi 


IHIHni 
inaHni 
^laoaHm 
iinacii^i 

ICJOHIHI 
IldOHMI 
inmMKtii 
lOHam 
■■■■rii 
iitiBhmi 

ICiaiMK^I 

iUHia»i 


lfr1BiF?l 

lUHini 

inHCill 


tsiaiFfli 

MBCII 

_ KIBIMII 

irflUCdHlKI 

inairiHm 

naHiuairfli 

IIViUiElBLJI 

inaiEiBiiKii 

IMHBlHtJI 

IBiHMBim 

MMMaPil 


ItdMnHHI 

IMfiKiMMI 

ir;fHMH»i 

mMPfHRI 


icai 


llrll 


inaiiaiiraaii 

lltiHKIBII-jUBI 

llrJHiBtIHirMrai 

inariM^Mi 

IfdarSBIKHil 

mariHriMi 

IMHRHL'Jni 

iirtHinamw 

inBiMHnHii 

iiiBiciHinm 

luaimHCiiai 

■■■■iaiFnBii 

■■■■■riHi 


lUl 


oaa 

nia 

■lUBII 

■inai 

■IKW 

■uiMKaaiHiaiai] 

■iiaBnaBHHiBii 

■icsaiEaBiHiaiaii 

\nMm 

■mil 

■IBII 

■IHDHIHiHHil 

■loaMiDBiiaBiHiHaii 

naiMMai 

l■i■i■ICI■i 

laariMHi 

l■l■iE»J■i 

IQIHIHI^E] 

■IB 

■■ 

■■ 

■iBi 

BSBBBBir 

■MBIHBIBil 

BBBEBBI 

■■■IBB 

Mfl  BB  ki^  Ah  1 

■■■■iiiai 

laimaM 

Via 

■a 

■in 

■lIDMiaBiBil 

■ICIBiElpHI 

■maiCiHCii 

IHoaiBiaHiiiaiPiaBiiDBi 

■ifiaiHnaMHiaMciai 


ssiBnnRS» 


M-t  'Ann  39ttlt 
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Uitil  uuu  u  u 
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muiBhwwi 


iSSSShrhi 

»Bmis]S 


mi 
^Baom 


IMKtfliMHHtsll 
Wi-1Wfc<Httg]l 
liaianEannaMi 
innBanui 
Kdcai 

llPJkii^l 

_ 

•^■WNKilyDMi^l 
■dUIMIBUir^l 
iKiEicmacmf 


SSEIrai 

■§■11 _ _  _ 

ISB8gggESB^gEl5gg£Eg&. 

' '  ~~'.g8SgSlgBSBSSggBSi 

,^_.i3B:lBlBKIBBa||Hil!^1HHV^MHHHi»BiH! _ _ 

iBgBBgSBBSBBBEKBBBSBBKS 

[aBnHUU3HaHiBnnHII!IBniBBgBBBBIgBH»Hi 


Ht!jaiBiii’iMaigHft^HiHiB:iaaiHggmaiH|g5M 

BBBMnapwMgiEigiBEapriwwBgigigigg 
nMirrwiPiwmnpiMWgiWgf 
mB<niWE3B5BBl:1fliHUnBBHlBHflBHBBH 

HDBiigBiHClHigigUHagignBaHiiigiHiBBiBilMSi 
nag!iBgira»wiwrigiPOEMB»wgiiggg 
laEIElBBiOBgigBBrieaBUPDgBgBBanBgBBMBiHi 

IPQtiUBldaiBBBiPnBgHPBaigBI  ~ 

ocag^naimisifliHivi-iBBBigaDiSgiaii 
:ia|aBitiHDBBgiF4k1PBiUDDBI~ 

.  annPIHmHiBragiBiMBiHiBBBBI 
QBtiMaiOBiaiHigiaMBiHinBii 
panUHIBHIiaaBBiNMBESBI 
PWRIMgiggwBggMiBBgBWPgM 
aanMHiiaHiHiUBBBMBBigHIE3l 

nSN|MSSilSSnSBSSSiaSSlili^S&aii|gii|SSI[B 
qigWMMgggiBWMWPBiWPPPag3BgiBaMB»i 
BPiaMwwtagiwwii^giEEwBawfS 
qfSitwiaigiaaiiiiMiHimL-iBpaErailWHiaiHigigifiS 
iC]aBMgragqBiPHiBiBgi-iBiQDr;ir*'ii  '  ~ 


s  e  R  n  R  s  ic 


1&  U  '.1  1920  21  22  23  2«  25  26  27  28  29  30  31  32  33 


^mumm 

nniilHIHil 
EIKlIigi 


Igffijgtl&ig 

l^fclUkiT 

|{3E3lidiidi 


BiSimsis 


i5i55555£gS5£saS5£555 

nSBBiiiBBSBSRBSBaKaS 

IBSBggSBBSBSBSSBSBB! 


IZIBnHHU 
IBBBEam 


■■■Ol 


BilBtfnHIUESiaFIHMBClBIIBiiMiMBag 

EEgBESBaaBSSaBBRBBBSB 

:in»4iaHii7iaiBKaHtfiiaiaiaiBi>iijBBiiiliii 


iriHtdHHHiHiaBai 

_ IpqaiBIBBHBMBBBHBMBB 

■■aBKlBBiqpBIB— ■■■■■■■!■■■■ 
BBViaiK^BcnMBnaiBHimmHiMiammBB 

■mnBiiMaiHmMPiBBuuBBBPHiBBMBigg 

BBMBrmisnaitiBiBiaaaiBciBiflHBBi 
■nBF4Br3BBiJBBDQBBCnBBBB 

;SSwSSr-SSSSiRBBSSBBBBI 

PBnHSBfliBJBileBBBMBBMBI 
'■■il.iBMBBftlBBBBBBBBMBiBBBI 
BBK^BBBBBBBBBBBBBt^BBBBI 
BBKIBBBBBHBBBBBBBBCIBBBBI 
MBMBBBDriBBBBBaBBMHBaBBBBB. 
BBBlBBBHnBBBBiBQBiBMFBBBBBBBBi 
BBIrlBBEIBIdBBBCIBBBBVhl 
BiBrJBBBBBBBBBDaBBBI 
BBUBBBBBBBDaiaaBBill 
BBBiBBBBBBBBBnilBBm 
BilBBBBBBBBBBBBRi^BI 
BiilBBifliaiaiBHiBBBBBBIIKI 
BBUBBBBBBBBBBBBBFIKJI 
iBBBBBBBBBfliBBBBBCJMiBMBBI 
'■CiMBBBBBBBBBBBBiKiCJicIBBBI 
IrJHBBBBnBBBBBBBnrilUBBBI 
iiftrlBBBBBBBBBBCaEIICIKIOiaBBI 
KBrJBBBBHBBIBBBiaiaitiriBDBBBI 
VilifBiaBBBBnBnnBBfHKJCaBBBBI 
B^mBBBBnBBBBBCItaiJITgrJEII 


|f4l 

■Ml 


E3I 

Kfl 


B! 


IBI 


nHElHtrlBI 


iMFJiiUHmi 
IPlUaWIHBI 


IKI 


caaiMBit!ii:ifli»edieap<iHiHHiHi 

BHItlHIUtEIHiliil - 


r-ii 


UHMHEiiniiilV:iHr«BdHI 
'eHiUHtMBISHMMHiKnni 

UnBWBCIldEiaHiWHilEMHiHHiHH 
,  iBWBKIEIMBBiaHaiilBHIHI 

aiMMwiiiaBSSBSiitaBBB 

IbBSSgHBBSHiaEiBSS 

IBBaiFIB.IElHiBIliBHmBBnH 
ic]Biainn>::aiHiBiaaiCiBHiBi 

SBaiBb  ir-iUBBBBBIdBBB 

IK3MBiCinnaiaiBiElBiK;lHiHai 

IBHMiUMMnaiHiKdHia-lHHiHil 

IMHMMIlMMHiaaBkJMHai 

liaHgUNKiHiUSiZiMMHialHi 

IHHlBMElfiaiDI - 


■I  ~ 


jggmaBgBBBMgggMjBg 


I 


VI  « 

>.  J 


O 
& 

X 

o 

U  COK  oj 
'o  isj 

Ui  vul 

O  .  e  fN-l 

^  X  s  »si 

Cfc  xi\f  s\ 

\  o\ 

•J  e  r«k. 

Oi  “s  tfs 


X  e  c‘ 
X  ^  f\- 
«vl 


CVi  «S.  (V  tv  t\.  « 

a  9, :»  a  ;9  o  i 

uu  a;!  uu  u/  X  oi  Ui  I 

c  C)e  ai  e  a  o  I 

Z  i  s.-.t  -si  j  -s  a  . 

S  ^  elc  s  s  £  e  I 

ST.r|r*.'s£v 
r  c. '  X  <“  I  £  £  I 

>■1  -v  tr'cc  If  I  ©  (V  *0  < 

.  c  -  '■**  ^  I 

x;  —  f-  >6  H  1 


•j  s  e  a  a<  a 
j  C  -  £  C 

'  O  *?  =•  s  r 


UtlRHTtNB  matrix 

soifr*rt  ■  I  .  “  MO  »fEreHTrMB“>»*TRr){ 


I  ^  ^ 

\  'K*.  Cg 

'  •  •  •  • 

•  e  a?*  « 
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PROGRAM  LISTING 
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\  OPTION  rOWTRAN 

T  FORTRAN  LSTOU-npCK 

r  MAIN  FI UTTFR  OVERLAY 

C  JAN  I J. 1067 

nOMM,  .i  |T(?1«  ) 

r  CHANOfc  BCD  TAPES  TO  BINARY  TAPES 
I  ft  CAM  LLINK  (6HPARITT) 

CAl.l  PAPTl 

CAM  LLINK  (^HPAR/'TT) 


CAI t  PART? 

00  TO  in 
END 

i  fortran  lstoii.ofck 

C  HPRINT 

SUBROUTINE  RPRINT  (  A  .  M.  N>  MD»  NTAPE 
n I  HE  NS  I  ON  A(1  ),  IT (6),  C<6) 

EOUIVALFNCE  (IT.C) 

P  FORMAT  MH  ,  4X»  ft(  ^X.  7HC01UHN  IM  )  ///  ) 

FORHAT  (1H  1  M,  X.  («SE  17.8)  ) 

N1  =N 
N2  =  6 
NJ  =  n 
N4  =  l 

4  IF  (N.T-NI)  6,6,«> 

*;  N?sN1-N<'»6 
N.TsNI 

6  K  =  n 

BO  7  1=  N4,N< 

K=K*1 

7  ?TIK)=I 

WRITE  (NTAPF.'’)  (  I  T  (  I  ) .  I  *  N?  ) 

no  0  1  =  1, H 

K  =  n 

I =M0»(N4-1 )♦! 
no  B  J=N4,N3 
K  =  K  +  1 
r(K)=A(L) 

«  l.sUMD 

o  WRITE  (NTAPF.M  1 ,  <  C  (  K  )  ,  K  =  1 ,  N?  ) 

IF  (N3-N1 )  jn.1  I,  M 
111  NJ  =  N346 
N4=N446 
0(r»'  4 

M  RETURN 
FNB 

f  FORTRAN  LSTOil,  DFCK 

r  MPIINCH 

subroutine  HPUNCH(A,N,N, IUUT, ITRA, |ORO,HrU?,HAXH#NTAPF»NCARUS) 
niHI NSION  A( 1 ) 

RETURN 

FNn 

f  LINK  PARI  IT 

f  FORTRAN  LS10U.0FCK 

r  PARTI 

subroutine  parti 

r  nsur=total  niimhfr  oe  surfaces  ailuued. 

C  NOFNSsTOTAL  NIIHBFR  of  nFNSiTIES  allowed. 

C  NRIOIDrTOTAl  NUMBER  OF  RIGID  BODIES  ALLOWED 

C  NSI7F  =  T0TAL  NIIHBFR  CONTROL  POINTS  ON  ANY  ONE  SURFACE  ALLOWED 

r  NHODFS^TOTAi  NUMBER  MODES  INPUT  ON  ANY  ONE  SURFACE. 


HN14007 
HH1 4007 
UM1 4007 
HH1 4008 
HM14008 
HH1 4008 
HMT4008 
HMt40C8 
HM14008 
HHI 4008 
UM1 4008 
HH1 400« 
HHMrtlO 
NM1 4019 
HRI 4809 
HH1 4009 

i 

HM1 4009; 
HM1 4009| 
HM1 4009 i 
HM1  4009  I 
MHI 4010 
HMMOlfl  ; 
HMMOlO 

HM1  401(1 
HM1 4010 
HM1 4010 
HH1 4010 
UHI 401  0 
MHI 4010 


njM>NSIflN  ISXST(/n),  ?SH(?U>.  hR(91),  SC/I).  RHO(?ll>.  UM((S»(S>. 

1  l)HBAR(  b  ,  >  .  RARHRH(b,  I/*  ) .  LOW<mi),  LMIGHC^fi) 

2  .  IT('^IH).  VELCTMXO).  NSlZtSC^fM. 

5  A(^.fi  ,  1  j,(i  )  ,  U(i>n,lU(l),  HP(50»A).  HRT(A»tun). 

4  HTC'S.l'in),  G(A.lUn> 

TMF  rolLOMlNG  STATFMFNT(S)  HAVF  RFFN  MANUFACTURED  BY  THE  TRANSLATOR  TO 

COHPFNSATF  for  the  fact  that  EQUIVALFNCh  DOES  NOT  REORnEP  COMMON - 

COMMON  IT 

FOmVALrNCE  (  I  T  (  I  ) .  I  SXST  ) .  {lT<?n.lSH).  (lT(41),RHn)  , 

1  .<  ?T(b1).NTAPE1  ).  (  IT(Ai;).NTAPF?).  {  1  T  ( 6  t  > ,  NT  APE3 ) , 

?  (1T(.,4)  .NTAPEA).  (  lT(At>l),NTAPE5) ,  (  I  T  (  60  .  nT  APEb  )  .  . 

3  (1T(67),NTAPE7).  (  I7(hB')»NTAPE0),  (  I  T  ( 6f> ) ,  NTAPE9  ) . 

4  (IT(/n).NSUR).  <  IT(71  ),NRI‘M0),  ( 1  T  ( /?  ) .  BREO  , 

'i  (  I  T( /3) -NAERO)*  (  IT(/A),NFUS)»  (  I  T  ( 7*) ) ,  mOENS  ) , 

6  <  1T(/A), MODES)/  <  n( /7).NP01NT).  «  I T  (  78  > ,  NPIJNCH ) » 

7  ( t T( /R)  .MAXR)  .  (  IT(Hu).MAXO),  ( I T ( 81  )  . Ma XS  ) , 

8  (  )T(tr7),NC),  (  IT(P3)  .NSURFS),  (  1 1  <  84  ) .  NPP  )  , 

9  (  ITf^^^I.HR),  (1T(»8A).S),  <  !T(1P7>,VELCtV)  . 
FQUIVALFNCE  (  IT(  |4/),NTAPEU),  (  IT(14H)  ,NRH0)»  (n(U9),NlTRSP)» 

1  (  I  r(  I  All  ),NITRnP).  (IT(181)»FPSP),  (IT(lf.p),FPnP). 

P  <  IT(  lAV), PICON).  (  IT(lt»4),NVEL).  (  I  T  ( 1 ‘jp.  )  .  NCARDS ) , 

3  (  IT(  )86).AITKFN).  <  fT(i87).NG0).  (F.G>,  (F(AU1,n,HT) 

(U.LfiH).  (IJTbl  .1  )  .LHIGH).  (  IT  (  1 88  )  ,  F  K  >  ,  (  |  T  ( 1 1>9  ) .  SQR 
8C0N)  .  ,  (IT(IAO).NSIZFS).  ( 1 T ( 1 80  ) , DM ) .  ( PM » AR , BARMRR ) 

6  Nj  .  (  1T(P16)  .NDELH),  ( ITCifl  /  ) ,  A  ITKED) .  (  IT(?1H)  .tsPARl  ) 

1  FORMAT  <1814) 

'  ’  format  (6E1’.8) 

3  format  ‘(VH  16X,  41H  FLUTTER  ANALYSIS  BY  A  COLLOCATION  METHOD 
1  ^47H  USING  AERODYNAMIC  INFLUENCE  COEFFICIENTS  ///JOH  NSUR  = 


FQUIVALFNCE 

) 

P 

3 

4  . 

bCON) 

ft  Nj 

1  FORMAT  <181 
7  FORMAT  <6E1 
3  format  ‘<10 
1  ^47H 


?  II?.  lOH  NAERO  s  114.  UH  NRIOID  =  11?.  9u  nFUS  - 

3  11?.  HtH  NDENS  s  114.  14H  MODFS  OUT  s  1l9 

4  .  9,H  NDELM  =  lI?,inH  NPUNCH  s  1  I  r-  ) 

4  FORMAT  <lHn.7vX,  UH  B  <REF)  s  lE?n.8.  8X.  4HK  s  lEPr.ft./iHO  P8X> 

1  ZHSURFACE  )8X.  ihB  IVX.  IMS  10 X . POHFXTERNAL  STORES  SIZE  //) 

8  FORMAT  <lHti  mx.  ?l'H  B  RIGIO  COMPONFNT  s  lFl8,8,  br.  bH  S  RIOU) 

1  .i3H  COMPONENT  =  lElb.H  ) 

A  FOPMnT  <IH  OPR.  ?<  5X.  1F7H.8),  lU?  ) 

10  FORMAT  '  <1H1  4HX.  1?H  .MASS  MATRIX  ) 

11  FORMAT  <41Hn  number  OF  CONTROl  POINTS  THIS  MATRIX,  t  114. 

1  48H)  AND  TOTAL  NUMBER  OF  CONTROL  POINTS  EXPpCTED.  <  1l4, 

P  37H)  no  NOT  AOREE.  PROGRAM  CONT I NOEO . . . .  > 

17  FORMAi  (IHl  40X.^4H  RIOIG  BODY  MODAL  MATRIX  ) 

M  forma:  <1H  .vHX.  8HSURFACE  1I7.  -IH.  116.  18H  CO'TROI.  POINTS) 

14  FORMAI  (1H1  41X.  PiiH  FLEXIBILITY  MA.THIX  ) 

l‘i  FORMAT  (IH)  46X.  18H  HEIOHTING  MAtRIX  ) 

16  format  (IHI  ^MX.  Pyp  rigid  COMPONFNT  AERO  MATRIX.  It9.  HH  CONTROL 

1  /H  POINTS  )  i 

17  format  (liil  .MX.  POH  AFRODYNAMIC  MATRIX  8X,  lOH  J./KlR  = 

I  n  '» (1 . 8  ) 

18  format  (IHii  1"X.  P.^H  RIOin  COMPONFNT  MUPES.  1110, 

1  1  /H  CONTROl  points.  )  ’ 

l»»  format  (IHII  ShX.  PvH  RIOin  COMPOKFNT  HASS  MATRIX  ) 

^11  FORMAT  (47111  FHROR  (N  INVERSE  ROUTINE.  PROGRAM  TERMINATED  ) 

n  FORMAT  (  111  38X.  HHSURFACF  IIP.  IH.  ftX,  13H  NO  ..'EIGHTINO 

1  7M  MATRIX  ) 

77  FORMAT  (  3M  7.8.  714  ) 

71  FORMAT  (  8H't|-P5P  =  1E16.8.  iX.  «H  FPHP  =  1E16.8,  iX, 

1  IflH  ATTKKN  r  1E16.H,  IX,  lOH  NITRSP  =  114,  JX, 

7  I  OH  NT  TRDP  =  114  ) 

74  format  (IHII  o;)X.  ?i.H  generalized  MASSES  ////  (l-i  .jdX, 

1  8H  MASS  1 14,  3H  *  IF16.8  )  )  i 
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J'S  roRHAT  (  IMfl  IttX,  23H.  RIGID  GOHPONfMT  MODtS  ///  )  ”79 

FORMAT  (  A<;X,  IA6*  114  )  -  080 

’n  FORMAT  <  IHft  T.1X.  23H  PUNCHFD  CARDS  NUMBERS  1A6#  1l4,  6M  T.HHII  081 

1  1  Aft.  1  14  )  082’ 

>'8  FORMAT  (IHO  ?«X.  3H  K*  lElA.R  )  083. 

DATA  OOn«CT/02n20, 104400(19/ 

Rcnz  aOnnoCT  085i 


%  OPTION  FORTRAN 

\  FORTRAN  LSIOII.DECK 

r  MAIN  FlUTTFR  OVERLAY 

C  JAN  13.1967 

COMMON  IT(21R) 

C  CHANGE  BCD  TAPES  TO  BINARY  TAPES 
in  CALL  LLIJIK  (ftHPARlTT) 

CALL  PARTI 

CALL  LLINK  <ftHPAR2TT> 
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f  OPMON  FORTRAN 

%  FORTRAN  LSintl.UFDK 

C  MAIN  FLUTTfR  OVFRLAY 

C  JAN 

nONMON  lT(?1R) 

r.  CHANOE  BCD  TAPES  TO  BINARY  TAPES 
in  CAll  LLlNK  (iSHPARlTT) 

HAIL  PARTI 

CALI  lI.INK  (ftHPAR.^TT) 

CAM  PARTi? 
on  TO  in 
FND 

t  FORTRAN  LSTOll.DFCK 

C  MPRINT 

SUBROUTINE  NPRINT  ( * . M, N. HD. NTAPE ) 
dimension  A(1  ).  IT(ft).  C(6) 

EQUIVALENCE  (IT.C) 

7  FORMAT  (IH  ,  4X,  M  AX.  THCOLUMN  IH  )  ///  ) 

format  (ih  ih,  x.  (6E  1/.8)  ) 

NlsN 
N2s6 
N3  =  6 
N4s1 

4  IF  (Ni-Ni)  A.A.'S 
D  N?sN1»N3<6  . 

N3sN1 
A  Ksn 

DO  7  I*  N4.N' 

K  =  K*1 

7  I T  <  K  )  a  I 

WRITE  (NT APE,?)  ( I T ( I ) . I =1 . N^  ) 

no  D  1*1. M 

Kan 

L=M0*(N4-n»I 
no  8  JaN4,N.T 
KaK  +  1 
C(K)aA(l  ) 

8  L=L*MO 

9  WRITE  fNTAPF,^)  I  .  ( C  (  K  ) .  K  =  1 .  N^' ) 

IF  (NJ-N1  )  10,11,11 

in  N;SrN3  +  A 

N4aN4+(S 
GOTO  4 
n  RETURN 
END 

t  FORTRAN  LSTOll.DFCK 

n  MPUNCH 

SUBROUTINl  MPIINCHI  A,M,N.  lOUl  .  ITRA.  I  ORO.MCDZ .  MAXM.  NT  A^E  .  NCARUS  ) 
n I  HENS  I  ON  A( 1  ) 

RETURN 

FND 

%  lINK  PAKITT 

%  FORTRAN  LSTOll.DFCK 

C  PARTI 

SUFROUTINE  PARTI 

r  NSURaTOTAL  NIIMF1ER  OF  SURFACES  ALLOWED. 

C  NDF  NS  =  T()TAI.  NIIMRFR  OF  OFNSITIES  ALLOWED. 

C  NRloniaTOTAl  NIIMHER  OF  RIGID  BODIES  ALLOWED 

C  NSI/FaTOTAL  NllMRFR  CONTROL  POINTS  ON  ANY  ONE  SURFACE  ALLOWED 

C  NMOIiF.S=TOTAL  NIIMBFtR  HDDES  INPUT  ON  ANY  ONE  SURFACE. 


I 

I 

I 

I 

I 


HH1  4(1071 
HM1 4U07" 
HMl  4  0(r/R 
HH1  4()0b|[ 
F1M1  4IJ0hI 
FiMI  40  084 
F4M1  4008* 
IIM1  40081 
HMl4n08/ 
HM1  4(II)H(^ 
F4M1  4II0  8T 
MM1  4n09(P* 
HH1 40091 
HM1  4009.'!’ 
HM1  40U9.J, 
HM1 40094 


F1M1  4009(1 
HhV4nii9r 
HM1 4009b 
HM1  4n09‘1 
F1M1  4010(1, 
HM1 40101 
HMl  40]  nil 


HH1 40104 
HM1  40int) 
MM140in(T 
HM1  4010  /* 
HM1 40108 
HM1  40inst 


Oin  UB 
0  0  ?  u| 
no.To 

0  0  4  (I 
0060 
(i  *  II 


niMfNSION  ISXSK/ii),  ISMJVU),  HR<?1),  RHO<?U),  DH(6»ft), 

1  nMBAR(6,l?>),  BARMRR(6.l?),  LOW(‘.0),  LHieH(;>0) 

^  .  IT(9JR).  VELCTY(ifn).  NSIZES{i»«). 

3  A()n,i.in).  Ft50.in«),  U(iO.lun),  HR(i*0.A).  lRT(6ilUO)» 

^  HT('^.lun),  G(A.IOO) 

r.  THF  Foil  OWING  STATFMPNT(S)  HAVE  BEEN  MANUFACTURED  BY  THE  TRANSLATOR  TO 

F  FOHPFNSATF  for  THE  FACT  THAT  EQUIVALENCE’  DOEL'  HOT  REOROER  COMMON - 

CUMHON  I r 

FOUIVALENCE  (  IT(1  ),  ISXST).  {  ITOD#  ISM>,  MT,<41>,RHn) 

1  .(IT('.1  ).NTAPE1>.  (  !T(A2),NTAPE?>,  (  I  T<  6'« ) ,  NTAPE3  ) , 

?  (IK/i^I.NTAPEA).  ( IT(At>).NTAPE^),  ( I  T(6«  )  *  NTAPEA  ) « 

3  ( IT(67),NTAPE7),  (IT(A8),NTAPF8).  ( I T ( 6o ) , NTAPE9 ) , 

4  (IK^ni.NSUR).  (IT(7|  ),hRIOID),  {m7?).BREF), 

^  (IT(/3),NAER0),  <!T(/4),NrUS).  (  I T  ( 7*) ) ,  mOENS  ) . 

f>  (  !T( /A),M0DES)»  ( tT( /7),NPU|NT)>  (  I  T<  78  ) .  NPIINCH  ) . 

7  {IT(/9),MAXR).  (  IT(«0)»MAX0>,  ( I T ( 81  )  * M ^ XS  ) . 

8  ( I  r(H5>),NC).  (IT<R3),NSURFS).  { I  T  ( 84  > ,  Ni>2  ) . 

9  (  I  T(u»;).HR).  (IT<1UA).S)*  (  IT(li’7>»VELCTY) 

FOIJIVAl  ENCF  <IT(  I47),NTAPE0),  ( !  T(  148 ) .  NRHO ) .  (  I  T  { 1>*9 ) ,  N 1  TRSP  ) » 

1  (lT(l'>(i),NlTRnP),  ( |T(151)*FPSP>»  ( I  T  ( 3  R?  ) ,  EPDP )  # 

?  (  fT(  1*1.5), PICON).  ( IT(1t>4),NVEL).  ( 1  T(  1  bR  > ,  NCARDS ) , 

(  IT(l*i^).AITKEN).  (IT(157),NOO)»  <F,0).  ( F  ( AU1»  1 ) ,  HT  ) 
4*  (U.LOU),  <U(i»1  ,1  I.LHIGH) »  (  I  T(  158 )  >  FK  ) ,  (  I  T(  1  &9  ) ,  SQR 

,  ( IT(1A0),NSIZES),  ( IT(1B8),DM),  (DHf»AR,BARMRR) 
ft  ,  <IT(21A),NDELH).(IT(2,17)>AITKED),(It(?1A),KPART) 


FOIJIVAl  ENCF  ( 
1  ( 

?■  ( 

3  ( 

A  . 

*)OON) 

ft 

1  FOPHAT  (18M) 
FORMAT  (6EI’. 
3  FORMAT  (iH  16 
1  4?H  IIS 

‘P  II?.  I 

3  II?,  J 

4  ,  9H  N 


(6EI’.«) 

(iH  iftX.  41H  FLUTTER  ANALYSIS  BY  A  COLLOCATIan  METHOD 
4?H  USING  AERODYNAMIC  INFLUENCE  COEFFICIENTS  ///.I  OH  NSUR 
II?,  iflH  NAERO  a  114.  11H  NRIOID  *  1I2,  Vw  NFUS  * 
II?,  JOH  NOENS  a  114,  14H  MODES  OUT  *  Il9 
,  9H  NPELM  a  lI2,i;iH  NPUNCH  a  il2  ) 


4  FORMAT  dHO  ?vx,  ilH  H  (REF)  a  lE2fl.8,  t»X,  IHK  *  IFZ-^.fl  /IHI)  25X, 

1  7HSURFArE  18X,  1 H8  19X,  IMS  inX»ZnHFXTERNAL  STORES  SIZE  //) 

*.  FORMAS  (lH<i  mx,  ?)H'B  RIGID  COMPONENT  a  IE18.8,  !>v,  8H  S  RIOID 
1  1.3h  COMPONENT  a  1E1H.8  ) 

ft  FORMAT  (IH  1129,  ?(  5X.  lF2fl,A),  111?  ) 

in  FORf.AT  (IHl  48X.  1  ?H  MASS  MATRIX  ) 

11  FORMAT  (41H0  NUMBER  OF  CONTROl  POINTS  THIS  MATRIX,  t  114, 

1  4flH)  AND  TOTAL  NUMBER  OF  CONTROL  POINTS  EXPECTED.  <  1l4, 

?  37H)  DU  NOT  AGREE-  PROGRAM  CONTINUED....  > 

1?  FORMAT  (IHl  4'^X,^4H  RIGID  BODY  MODAL  MATRIX  > 

11  FORMAT  (IH  .18X.  8HSURFACE  H?,  IH,  116,  ISHCOhTROI  POINTS) 

14  FORMAT  (IHl  41X,  PIIH  FLEXIBILITY  MATRIX  ) 

lb  FORMAT  (IHl  4ftX,  1«H  WEIGHTING  MATRIX  ) 

1ft  FORMAT  (IHl  Pn%,  pgw  RIGID  COMPONENT  AERO  MATRIX,  1?9,  «H  CONTROL 
1  /H  POINTS  ) 

17  FORMAT  (IHl  MX.  2nH  AFRODYNAMIC  MATRIX  8X,  IDH  l./K  R  a 
1  II '’n. 8  ) 

Ift  FORMAT  (IHii  mx,  ?.1H  RIGID  COMPONENT  MODES,  1110, 

I  WH  CONTROL  POINTS.  ) 

I'J  FORMAT  (IHh  mx,  P^M  RIGID  COMPONFNf  MASS  MATRIX  ) 

,M1  FORMAT  (4;h*  FRROR  IN  INVERSE  ROUTINE,  PROGRAM  TERMINATFD  ) 

:'1  FORMAT  (  1M  .mx,  RHSURFACF  112.  IH.  6X,  t3H  NO  -EIGHTINO 


'?  format 

yA  FORMAT 
1 
2 


7If  MATRIX  ) 

(  3T12.M.  ?14  ) 

(  HH'ifPSp  =  IFIft.M,  IX.  «H  EPIJP  =  1E16.R,  iX, 

KIM  AITKFN  =  lElft.fl,  IX,  lOH  NITRSP  *•  114,  IX, 
1  OH  MI TROP  a  114  ) 


/4  format  (IM..  4fix.  P(iK  0FNERALI7ED  MASSES  ////  (lu  jnX, 


bH  MASS  114, 


lFlft.8  ) 
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/8 

FORMAT 

( 

1  HU 

5  IX, 

21H 

RIGID  OOHPONENI  MODES 

///  ) 

■1  /9  0 

?6 

FORMAT 

< 

?  I  4  . 

6  2X, 

1  A6, 

114  ) 

•‘boor 

2  7 

FORMAT 

( 

1  Mil 

3  )X, 

?3H 

PUNCHED  CARDS  NUMBERS 

1A6,  Il4.  6H  THRU 

‘•HI  ol. 

1 

1  A6  . 

1  14 

) 

1!  HjMf 

28 

FORMAT 

(IHfl 

211  X 

,  3H 

Ks 

IEI6.0  ) 

••  h  <  0  '  - 

OATA  onuocT/0'->n?n.^n44i)no5/ 

Rcn?  sonnncT 


NTAPno  S  PUNHM  OU 
RT'.Pe?  =  INPUT  TA 
NTAPE3  r  OUTPUT  P 
NTAPF4  r  / 
NTAPF*)  s  /  AN 
NTAPEA  s  / 

NTAPF/  s  / 

NTAPF8  s/  • 
KTAPF9S/ 

IF  (  Nn0-987/,h  ) 
97  ir  (  NAFRO  ) 

)fl  NAFRO  r  NAFHO-1 
IF  (  NAFRO  ) 

99  NTAPFn  s  8 
NTAPFVa^ 

NTAPF.T=^ft 
NTAPF4  ■=  9 
NTAPFN  s  i 
NTAPE6  s  4 
NTAPE7st1 
NTAPFBsl 
NTAPF9  slO 
HAXR  a  •»() 

HAXO  s  6 

NCONail 

MAXSS5A 

991  REMINO  NTAPF4 
REMIND  NTAPER 
REMIND  NTAPF4 
REMIND  NTAPF9 
SORCON  s  SORTT  ^ 


PUNOM  OUTPUT  TAPE 
INPUT  TAPE 
OUTPUT  PRINT  TAPE 


ARE  UTILITY  TAPFS 


90,97,99 

98,170 .9H 

99,99,1/9 


SORT<  186.0)  /  (7.n«<.14H»9  ) 


EPSP  =  .6E-i)6 
FPDP  =  .''E-(l7 
AITKFN  =  .9 
AITKFD  =  .9 
PirON  =  .*,9;M( 
NCARDS  =  )l 
NVFI  =  II 
NCn 

NITRSP^-lli 
NITROPalH U 
RHO(  1  )  =  II.  II 


‘,9;m  •o.n-l.l  4189 


r.  READ  IN  TITLF,  FONTROLS  ANO  CONSTANTS  AND  PRINT. 
lUfl  CAM  RDI  N  <NTAPFv,ntAPF.T,1  ) 

READ  (NTAPi  7, 1  )NSUR.  NAFRO,  NRIOIO,  NFUS»  NDENS,  HQi'ES 
tNDEIN.  NPUNFH,  NfON 
READ  (NTAPI  •>,?  )FK 

MRIIF  (NTAPI  I.DNSUR. NAFRO,  NRT6T0.NFUS,  NDENS, NODES 
1  ,  NIlFIM,  NPUNCH 

NR7rNRI«l  t)*Nf! 

NSIIPFSsNSUR*  Nf  US 

If  (  NOON  )  in2,10'«,»07 

III?  REAli  (NTAPF.?,??)EPSP, FPDP, AITKFN,  Nl  TRSP,  N I TRDP 


OES,, 

I!  8h  0 
•18  b  0  j 

0BOu-‘ 
090-' 
•.K'l!.  r 

"9?li]. 

09Jii 

0  94^, 

0  9bi)  } 
0  961)“'* 
0  97!) 

0  9HI1  'J 
ii  9  9 1)  ^ 
1  0  n  tl 

I 

’•  0  ■«  ■; 

1  4  r, 

1  0  *3  'J  ;{ 


1080" 
1  090^ 

n  no 


111  0 ' 

1130 

n4tiT 

n  8  ij  ± 

'1  16  0 
117  0^ 

n8o| 

1190 
i?no_ 
i?i  0 ; 
1  i;?o-^« 
1  ?3u 
1  ?4li  ** 
1  ?80^ 
1260 
5  27  ll;- 

:  j  8  0  ! 

129J  '• 
130  0 
1 .51  0 
1  v5?0  !. 
1  33  0 
1  34  0  ■ 
138  0  i 
1  360  ■" 
1370 
1  38  0  i 
1  390  ^ 


I 

I 

I 

I 

I 

! 

1 

I 

I 

Cf 

I 

1 

I 

I 

I 

I 

1 

1 

1 


RtAf)  (NTAPf'>.??)Af  TKtn 

MRITP  (NTAPn,^.i)tPSP.  EPOP»  AITKEH#  NITRSP,  NITROP 
104  IF  (  NAPRO  ) 

103  WRfTF  (WTAPfi3,'>H)FK 
nOTO  111 

10^  READ  (NTAPfcP.1 )(ISXSTri)*ISU(f ).I«1*NSUR) 

READ  <NTAPFi»,5>)BREr.  ( VELCTV(  I  ) » I«1  ,NAERO > 

NCs? 

RR?r  NR16ID  «RC 

READ  (NTAPE?.P)(BR( I).S(I )<|:1*NSURrS) 

IF  (  NDENB  )  106,1U7*1UA 

10^  READ  (NTAPEP.P)(RHO(l).I«1«NOFNS) 

107  IF  (  NFUS  )  10ft. 100,106 

1')«  WRITE  (NTAPE^.'ilRHd  ),  S(l) 

109  MRITF  (NTAPFI.oirrPF  .FK 
DO  no  IsI.NSUR 

J=I*NFUS 

110  WPITF  (NTAPFT.A)  I,  RR(J).  S(J)  .iSXST(i) 


READ  IN  FIISF'  OF  HASS  CHARACTERISTICS 
111  IF  (  NRfjlll  )  11?, 117, 11? 

11?  IF  <  NDELH  )  1l5.1i:i,11S 

1  13  DO  1  14  1=1  ,  NR  I  OH) 

no  114  J=1 .NRIOID 
114  DMII.JIsO.O 

OOTf)  117 

1  I*?  DO  lift  1  =  1  .NRIOID 

lift  READ  (NTAPE?,?)(DN(.I,I),J=1. NRIOID) 

WRITE  (NTAPE1,19) 

CALI  HPRINT  (DH, NRIOID, NRIOID, HAX0,NTAPE3) 

READ  HASS  HATRIX  FOR  EACH  SURFACE,  STORE  SYSTEH  HASS  HATo)X  ON  NTAPE4 
1  17  Kl  =  i| 

WRITE  (NTAPFT.I'I) 

DO  1?n  Isl.NSMR 

RTAD  (NTAPE^d  INSIZF 

READ  (NTAPE'^,1  )  ( I  OW(  J ) , LHIOH( .1) ,  J=1 ,  NSl ZE  > 
nn  119  Jsl.NSlZF 
DO  lift  Ksl.NSiZE 
IIP  A(K.J)  =  ii.)| 

Nl=|  OW(  J) 

I H I HH ( J ) 

119  RtAD  (NTAPE?,’)(A(N, J),NsHl,N?) 

WRITE  (NTAP(4)NSI/E.  NSIZE,  C ( A( N. J) ,N=1 , NSI ZE > , J«1 . NS  I ZE ) 

WRITE  (NTAPn.l  ni.  NSIZE 

CALL  HPRINT  ( 4 .  NS  I  ZE .  NS  I  ZE ,  HAXR,  NT  APE.i  > 

NSI/FS( I IsNSIZF 
1?n  Kl=KI*NSIZF 


NPOINT  =  TOTAL  NlIHRER  OF  CONTROL  POINTS  ON  ALL  SURFACES. 

NP01NT=K1 

READ  IN  RiOin  BODY  HODAL  HAIRIX  FOR  EACH  SYSTEH. 

IF  (  NFIIS  )  1?t,121,121 

1?1  READ  (NTAPF'>.  1  )NSIZF 
DO  1??  Jsi. NRIOID 

W?  RIAI)  (NTAPf  9,'^  )(HR(  I,  J).  i  =  l  .NSIZE) 

WRITE  (NTAPI 1, iftINSiZF 

CAM  HPRINT  (MR. NSIZE. NRIOID. HAXS.NTAPF.i) 

WRITE  (NTAPI 4 )nSI/E.  NRIOID,  ( ( HR( I , J ) . I =1  . NS  I ZE ) . J= • . NR  I GI 0 ) 
1^3  IF  (  NRIOID  )  1?4,1i?B,1?4 

1/4  Klsl 


14flU 
141U 
14?0 
1430 
1440 
1450 
1460 
1470 
1  460 
1490 
1500 
1510 
15?« 
1530 
1541) 
1551. 
15f  li 
1570 
1560 
1590 
1600 
1610 
16?U 
1630 
164  0 
1650 
1660 
1670 
16BU 
1690 
1  700 
1710 
17?U 
1730 
174U 
1750 
1  760 
1770 
176* 
1  790 
1  fcl'O 
1610 
16?0 
1  630 
1650 
1860 
1670 
IflftO 
1690 
1  900 
iVtU 
19?0 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
?010 
Z0?0 


— 
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I 

I 


MRIIF  (NTAP(^.i^) 

no  1?^  isuRrt.Nsnp  '/’o-iij.- 

PPAn  (NTAPFV.l  )NS1ZF  VWbwl 

K/  =  K1  ♦NSI/F-1  POftU*' 

no  ip*i  j^i.NRiGin  '^rj\] 

RFAI)  (NTAPEP.P)  (MR<  I  .  J)  .  1=K1  .KiO  PfiRii*' 

MRirF  (NTAPF^.I  OISUR.NSIZF  Vo*?.)*. 

HALL  MPRINT  (  HR(  K1  , 1  ) .  NS  I  ZP  .  MRIG I  D  »  HAXS  .  NTAP63  )  ?lf.u 

ivft  Kt=KltNSI7F  ''nil-- 

NRITF  (NTAPFS)NPOINI  ,  NRIQIO.  <  {MR(  H.  J ) ,  N  =  1  ,  NPO I  NT  ) »  i  =  1  » 

1  NRlGln  )  '/'Uit*' 

IF  (  K1-NP01NT-1)  1?7, 128,127 

127  MRTTF  (NTAPK1.1  1  )K1  ,  NPOINT  ^1^-; 

•  •••••#*«»*«•««»•«•»»«•««»•••«•»•**«»«•••*••«»•••«•••••«• 

RFAO  IN  FLFXlBn.1TY  MATRIX  FOR  tACH  SURFACE.  STORE  -N  NTAPfc5  23^0 

128  Kls<>  219!" 

MRITF  (NTAPF-i.M)  2?ii(l;, 

no  1*>3  In.NSUR 

RFAH  (NTAPI ’.I INSIZE.  J.  IFORM,  IRONS  2??e- 

CATI  HRFAO  (  A.,  NSI  ZF .  NSI ZE .  I F  ORH,  IRONS.  II .  1  »  F  .  HAXR,  NTA>'b? .  NT  APED 
NRITE  <  NT  API:  <,11)1  ,  NS  I  ZE 
CAI,l  MPRINT  (A.NSIZF.NSIZE.MAXR.NTAPFi) 

NRITF  (NTAP.F8)NSIZE,  NSIZE.  ( ( A  (  J,  K  ) ,  J=1  ,  NS  I  ZE  ) ,  K^I  .  mSI  ZE  )  2?6IF 

1S3  Kl  =  K1  ♦NSIZE  ?28ui. 

IF  (  KI-NPOINT  )  ISO, 16-1, 1»>8  229'l 

1‘>8  WRITE  .(NTAPET.I  l)K1  .  NPOINT  vani.** 

164  fr  (  NAERO  )  166, 168, 166  2310 

168  NRHOsl  73?o' 

fib  6110  lal.NRIOin  Z3.H- 

no  500  Jal.NR-.'  ?34U 

6110  OQ(|.J)a'1.n 

WRITE  (NTAPF6)((nH(I,J),|a1,NRIGtD).Jal.NRlGlO)  236U 

NrXsNC  23«ir 

OOTO  PI  6  ?39l|. 

?4n« 

C  READ  IN  MEIGHTINO  MATRIX  FOR  EACH  SURFACE . STORE  ON  NtaPE6,  241 

166  .  WRITE  (NTAPE3,V,)  24?(- 

no  1/8  Isl.NSUR  243i;* 

IF  (  ISW(I)  )  167,177,167  2440. 

167  NIsISVSTd)  ?45( 

If  <  N1  )  168,172,168  246U. 

1  68  no  ^  J: 1  . N1  2470 

no  I'.R  l.  =  l.MAXR  24Rr- 

A(J,l.  )"n,ii  249l« 

169  A(l,J)  =  ll.(l  2800*' 

I/O  A  (  J  .  J  )  s  1  .  (I  2  8 1  II, . 

RIAII  (NTAPF'>,1  )NXS|.  J,  IFORM.  IRON  ^62i 

If  (  NXST  )  171,172,171  .'830- • 

1/1  FAIL  MREAP  ( A , NXST , NXST , I F OHM, I  ROM, 0 , 1 , F . MAXR , NT AP' 2 , NT APb3 )  284o 

1/2  K=N1*1  ?88(.f* 

RIAD  (NTAPI-2, 1  )NS1ZF,  NPART.  IFORM,  IRON  ?86t;I. 

If  (  IFORM  )  174,1/3,^74  28/0 

1/3  CAILMNfcAII  IA(K,K),NSIZF,NSIZF,n,n#0,l,F,MAXR,NTAf’E?,NTAPE3)  28RO.- 

K^K»NS1ZF  289u; 

GOTO  1  /6  2600"' 

1/4  nn  1/8  JrKNPARI  2610 

RFAn  (NTAPF2, 1  )NSIZF  26?nr 

N^K.NC  263Ui'- 

llf)  1/41  l.rK,N  264  0 

nil  1/41  Ms1,K  2680'' 

1741  A(M,l)  =  li.(i 


I 

f  I 


t 

t  I 


I 

I 

I 

I 

I 

■« 

t. 

I 

I 

I 

I 

1 

I 

1 

I 

1 

I 


I 


1  /h 
1  /<S 


71 


CAIL  HPKAi/  (A(K.K),NSIZe.HSI7e.  I  .  lROW.n,0.f  ,  HAXR»  NTAPtV .  HI  APfc  J  ) 
KsK4NSI7F 
K  =  K-1 


WRITE  (NTAPF^.l  Ol,  K 

TAIL  HPRINT  (A.K.K.RAXR*NTAPF3) 

MRITF  (NTAPEA)K,  K.  ( { A ( J» L ) . Jsl • K ) »Ls1 »K > 
GOTO  17B 

1?7  WRITE  (NTAPE’^.'/'DI 
17R  CONTINUE 
I /O  REWIND  N^APE^ 

REWIND  NTAPF7 
NCXsNC 

NVFl  s  NVEL+1 


IF  SERIFS  OF  DENSITIES  FOR  EACH  V/o  OMECA»  READ  IN  THAT  SERIES. 

IF  (  NIlFNS  )  1Bn,181.lBft 

ion  NRHOrNDENS 
OOTO  IBP 

1H1  READ  (NTAPE'^.  1  )NRHO 

READ  (NTAPF'’.?)(RH0<  1),  Ia1  .NRHOi 


1B3, 184,181 
184,188,184 

189,188,188 


READ  IN  COHPLEX  AFRODYNAMIC  MATRIX  FOR  EACH  SURFACE 
1HP  no  Plfl  Isl.NSURFS 
Ksl 
K?sl 

IF  (  NFUS  ) 

181  IF  <  I-l  ) 

184  LsI-NFUS 

IF  (  iSXSTd.)  ) 

188  K=ISXST(L)H 

KV  s  K+ISXST(L) 
no  188  Jsl.K 

00  188  L=1.KP 

1H8  A(J,L):*l.O 

RFAD  (NTAPE^»,1  )NXST.  J.  IFORM.  IRON 
II  (  NXST  )  187,188,187 

147  N~  NXST  *NC 

CAM  HREAD  < A.NXST.N, IFORH, IR0H,n,1 ,U,HAXR,NTAPE2,NT»PE3) 

148  READ  (NTAPF7,P)VFLn 

RFAD  <NTAPEP,1  INSIZF,  NPART.  IFORM,;  IROH 
IF  {  IFORM  )  l9n,18R,1Vn 

180  N=  NSIZF  •NC 

CAM  HREAD  ( A ( K . X  . NS  I ZE , N  , IFORM . I  ROW . 8 , 1 , U , HAXR , NT APE? , NTAPE3 ) 
NSI7FsNSlZE«K-l 
OliTO  19.1 

ion  no  19P  .1=1.  NPART 

READ  (NTAPFP.1 INSIZE 
N=K/»Nn 

no  191  MrK’.N 
no  191  I T I .K 

101  A(l,M)sn.« 

Ns  NSI/r  -NC 

CAM  MREAD  (  A  ( K ,  K/* ) .  NS  I  ZE ,  N,  1  , 1  ROM ,  0 , 0 ,  II,  MAXR,  NTAPEP  ,  NT  ARE1 ) 
KsK4NSI?F 
107  KPsKP*N 

NSI 7FsK-1 
NsKP-l 

10.1  IF  (  1-1  )  199,194,100 

104  IF  (  NFUS  )  199,19/,  t9*i 

19»>  WHITF  (NTAPE  1,1  i)NSI7E 


V?6  7U 
P680 
P690 
?700 
'.^710 
?7?a 
2730 

2746 

2790 
27.80 
2770 
2780 
2790 
2800 
2810 
?8?U 
28.30 
2840 
’890 
2880 
2870 
2880 

?91i0 
?V1« 

2920 
2930 
2940 
2990  ;  f 

2960  I 
2970  I 
2980 
2990 
3000 
.101  0 
3020 
3030 
3040 
3090 
.3080 
3070 
3080 
.3090 
3100 
3110 
3120 
3130 
3140 
3190 
.3160 
3170 
31 BU 
3190 
320  0 
3210 
.3220 
3230 
3240 
.129  0 
3260 
3270 


■  .  'a#  . ' 
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HAM  MPHINT  (  a  .  KS  I  7F  .  M.  MAXR.  NT  APR.n 


n  COMPIITF  (I)Q)  =  (HR)T  *  (CH)  *  (HR) 

Kir) 

KV-; 

M  =  NC  I 

READ  (NTAPt^lNSlZE.  L.  (  ( U(  L .  J  )  .  L  =  1 .  NS  I  Zh ) .  Js1  ,  NRI G 1 1)  > 
nAll  HHIILTO  (A,N.U,n.r.NSIZE.NSIZF>NRIbIU.HAXR*HAXR.HAXR) 
no  196  L=].NSjZE 
on  196  j=i.NRioin 
1  A(  J#L)tU<L. J) 

CAl.t  HNULTn  (  A.O  .r.N.  !)0.NR|0|0.NSt7E,NRlQlD.HAXR*NAXo.MAXU) 
nOTO  ?10 

1  97  no  1  9«  1  =  1  . NRIOin 

no  19R  Jsl.NRv 
19R  nO(L.J)  =  ''.A 

1 99  I =|-NFnS 

IT  <  isw(L>  )  9n?.?on,vii'> 

WRITE  (NTAPI-/)NSIZE.  N,  (.(  A  <  J  .  M  ) »  J=1  .  NS  1  ZH  ) ,  H«  1 ,  N  \ 

GOTO  ?nft 

READ  (NTAPE<^  )|  .  L.  (  ( F (  J . M ) .  J=l  . L > . Hsi  . L  ) 

IE  (  NSiZE-l.  )  9n4.?flA,, 

Pn4  write  (NTAPE  1.1  DNSlZE.L 

9ii6  CAI.I  HHULTO  (  E  .  •) .  A  ,  NC-I  .  U.  1. .  L  .  1. .  HAXR.  HAXR,  HAXR  ) 

WRITE  (NTAPF/)NSTZE.  N.  ( ( U( J > N ) . J=1 » NS  1 ZE ) # N  =  1 , N ) 

WRITE  (NTAPE.1, 1  7)VELC 
LsI-NEIIS. 

WRITE  (NTAPEI.i  nL.  NSIZE 
CAI.5  NPRINT  (  A.NSIZE.N.MAXR.NTAPED 
91.0  CONTINUE 

C»»«»»»**««**»«**«**»»» 

C  CARRY  ON  ERON  HERE  TO  END  ONCE  FOR  EACH  DENSITY. 

?1S  no  3(10  IRHOsi  .NRHO 
KsNC*NP01NT 
no  916  Isl.NPOINT 
no  916  J=1,K 

916  IK  I  . 

REWIND  NTAPFi 
REWIND  NTAPES 
REWIND  n^apf; 

CONsRIiOl  lRlin)*nR(  )  )#»?»S(1  )  *’Z.174 

IE  (  Nl’MlIin  •>1R,99?,?1R 
9)R  no  9^0  l=l.NRIGin 

no  99(1  J=l,NR9.NrX 

nMBAR(  I ,  .iM  )=noN»no(  i .  j+» ) 

Kr  J/NONC-I 

99(1  DMHAPI  I  ,.11  =  DM  (  I  ,  K  )»CON*DO(  I .  J) 


C  READ  ENTIRE  HR  MATRIX. 

READ  (NTaPPOI.I.  (  ( HR  (  I  >  J ) .  I  =  1  .  NPO  I  NT  ) .  J  =  1 .  NR  I  G  I  D  ) 

99?  Ki=n 

no  916  ISOR-I.NSUR 
,  K  =  ISUR»NEIIS 

CONsRHIK  lRHn)»BR(K)»»9»S(K)  74 

C  RF  AD  ;M)  I 

RIAl)  (►  TAPE4  )NSIZF,  NSIZE.  ( ( E  (  I .  J  ) ,  I  =  1 .  NS  I  ZE  ) .  1  =  1  ,  NS  1  ZF  ) 
IF  (  NAERO  )  991  , '991  , 

991  1)0  996  1  =  )  .  NS  I  ZE 

no  996  J=).NSI7F 
996  A( I , J)=F  { 1  .  J) 


I 

I 

I 

r 

I 


.19  Mil 
.3990 
.13(10 
..131  0 
'3911 
.1330 
3340 
.13^0 
.M6I1 
.137  0 
330(1 

339  0 

340  1) 
<410 
'490 
M430 
144„ 
‘4^H 

'  4  f.  ta,, 

I4f*0  I 

.14  9  0 
.^bOll  . 

3^ioir 

3»>9  0'*^ 
3^3  0 

3»,60| 
3S7l! 
3580_ 
3»)9oI 
1 6  (1  0  ^ 
.161  II 

*69ur 

36  3  0'. 
3  64  0 
361>1) 
366  0  ■ 

36  7  0”- 
1680 

37  0((.*» 
.1/10 
1/?()«*‘ 
3730  > 
3  74  0^’ 
>’710,., 
3  7  6  0 

3  7  7  tl->‘ 


3/80 

3790r 

36011' 

«» *. 

38)  u 
■1 R  9  0 

383  0  ■ 

384  u" 
3  880 
3861)“' 
3H/0i! 
3bH0 
3890 
390  01 


n  n 


1 

i 


/O 


I=NSIZE 
ROTO  9^■\ 

n»» 

n  RFAD  (W)*(nH>  I  ''  . 

9-^  R£An  (NT4PE7>NSI7F,  L.  (  ( A(  I .  J)  *  I  « NSl  ZE  ) «  Js]/.  L :) 

n  nOMPUTE  (M  BAR)  a( (H)*RHO«RR*»P«S*(W)»(CH)  )  I 

no  ?H8  IsJ.NRIZF 

no  ??8  jsi.L.N':x  , 

K=J/?+1 

A<  I,  J)s(:ON«A(  I,  J)4F<  I.K) 

A( r, JH ) -noN*A(  I .  j*n  ;  s 


RFAD  riFXIRILlTY  MATRIX  FOR  SURFACE 

?n  RFAfi  (NTAPfS)NSIZF»^-NSl[2E>  « (F  <  I .  J)  .  I  *1 ,  NS  I ZE  )  /  J*!  » I ZE ) 
K-Nr*K1*1 

CALL  HNULTD  ( F  .  0 .  A .  NC-r>  U( Kl«l  .K )  *NSI  ZE>  NS)  ZE *  NSI  ZE>  ^'AXR.  HAXR. 

1  MAXR)  . 

IF  (  NRIGID  )  P3?,?46,i>3P 

. . . 

C  FtNn.(  IITTLE  H  bar  >  s  (OFLTA  M  )  ♦  (H  R)TRANSP0SED  •  .(‘«  BAR)  •  (N  R) 
no  ?4n  iri.NsiZF  . 

XsKlH 

no  j=i.NRiGin 

?4n  0( J# I )=HR(K. J) 

KsNf;«K1*l 

CALI  HHIIl  TO  (G.O.A.NC-1  >HRT(1*K).NRI6li)*NSIZE*NS!ZE*‘'*AXU»HAXR» 

1  MAXO)  < 

CAlt  HHIU.T0  <HHT<1  ,K).NC-1  »HR(K1«m  )>ll»A,NRI6)D>NSfZE*NPI6ip»  ' 

1  HAXU.HAXS.HAXK) 

no  ?44  in.NRiGin 
no  944  Jsl.NR,^ 

94  4  RARMRRd  ..J):RARNRR(I.J)«A(),J) 

74^  K1 =K1 ♦NSI ZE 

L=NC*NP0INT  . 

IF  (  NRiGin  )  947,P6A,947 

947  ROTO  (?52,?4«),NC 
94R  no  P*>0  Isl  ,  NRIGID 

no  ?Bn  J=1,NR9.NCX 

K=J/?*1 

n( I ,K )rHARMRR( I , J+1  ) 

9*,n  BARMRR(  1  .K)=RARMRR(1»  J) 


c 

(; 


THFN  (IITILL  M  HARIINVFRSF  ANO  )INAL  U  MATRIX  STORED  ON  mTAPE4. 
CAM  MNVRSX  (  R A RMRR ,  G .  A  ( 1  . 1  )  .  A  (  1  *  HAXO  )  #  NK  I  G I  D »  I R *  NC- 1  ) 

IF  (  iR  )  3in,Pt»4,nn 

Rom  (  2MI ,  V'W. ) ,  NCX 
7',<,  on  p*>«  i  =  i.NRinin 

no  9Sft  JtI .NR/.NCX 


K=NR?-J 
M=K/?« 1 


9iR 
96  0 


RARHRHI I . K )3RARHRK( I .H) 

PARMRRI  I  .  KM  )=G(  I  »M) 

CAM  MMIlLTn  (RARMRR,Nr-1  ,HRT»NC-l.  F » NR  I G  f  0»  NR  I G I D ,  POINT, 

MAXO. MAXO, MAXR) 

CAl I  MMOl  Tn  (F,NC-1 ,u,nc-),mrt,nrioid,npoint,npoint.maxr,maxr. 

KAXO) 

CAl  I  MMULTR  (MR.  <.f .NC-I  ,A,NP01NT.NRTGIO.NPOINT,MAXR.MAXR.NAXR) 
no  ?<S4  1  =  1.NP0INT 

nn  9^?  .(=1.1 

F(  I  ,  J)  =  ll(  1  .  J) 


3V10 
3920 
3930 
3940 
3950 
3960 
397u 
3980 
3990 
4000 
4f,io 
4020 
4030 
4040 
4.050: 
4060 
4070 
4  080 
4  090 
410  0 
4110 
4120 
413  ft 
414i 
4l5o 

4160 

4170 
4180 
4190 
4200 
421  0 
4220 
4?30 
4240 
4250 
4260 
4270 
4280 
42s  0 
4300 
431  0 
4  32!) 
4  3  •. 

4340 

4  35  0 
4360 
4370 
4  38  0 
4  39  0 
4400 
4410 
4420 
4430 
444  0 
4  45  0 

4451 

4452 
4460 
44/0 
4480 
4490 


74 


A(  j ,  J)  =  -A(  r ,  j> 

K=  ?*NC-NC< ' 

A{|,K)  =  1.I!*A(1,K) 

CAl  I  MMtiLTD  <  A.NC-1  .  F  ,  NC- I  .  U .  NPO I  MT »  NPO  I  NT  ►  NP(I  I  NT »  HA  vR ,  M  A .  MAXK  > 
IF  (  NAERO  ) 

REAP  (NTAPf  *>)A(1  .1  ) 

NRI.TF  (NTAPFSX  (HRT(  I.  J),  1  =  1 .  NR  IS  1 1) )  ,  J  =  1 .  NPO 1  NT  ) 

MRIIE  (NTAPFR)NPOINT.  L.  ( (  U(  I .  J  > ,  I  =1  ,  NPO I  NT  ) .  J  =  1  *  L  ) 

I'ln  CONTINUE 

rewind  NTAPF‘) 

REWIND  NTAPFA 
rewind  NTAPF/ 
rewind  NTAPFH 
REWIND  NTAPF9 
KPART  =  ? 

RETURN 

Tin  MRTTF  (NTAPI  l,V(i  ) 

ROTO  9R1 

T.>n  REWIND  NTAPFi 
REWIND  NTAPF‘. 

IF  (  NRIGII)  )  140 

T  in  REAIi  (NTAPF*!)  I ,  I ,  Ad.l) 

TT?  Kls)  . 

DO  Isl.NPOINT 

no  .T34  Jsl.NPOiNT 
TT4  F(I,J)sti.(1 

DO  T34  ISURsl.NSUR 
K>‘sK1  ♦NSIZFS(  ISUR)-I 
READ  (NTAPFD)!,!,  Ad.l) 

READ  (NTAPEA  )NSIZE.  NSIZE.  ( (F(  I .  J ) .  IsKl ,  K;/ ) ,  J=K1  K?) 

T.1A  KJsKI+NSIZF 

READ  (NTAPEA  )NSIZE,  MODE,  < < A< I . J ) . I *1 . NS  I ZE ) . Jrl , MnDF ) 

MsMODE/NC 

CAM  MMULTD  (F,n,A,NC-l,U.NP0fNT.NP01NT,H.MAXR.HAXR.MAXR) 

IF  (  NDELM  )  .T;i7,  T4R/<37  , 

T17  IF  (  NRIOID  )  .T3fl,14fl,T3R 

TIP  READ  (NTAPEH)((DH<1,J),Is1jNRI&IU).Js:1.NKIGID) 

READ  (NTAP>  H ) ( (n( I , J)  , l  =  l .NRIGID) . J=1 .NPOINI ) 

DO  T.10  1  =  )  I  Gin 

DO  '.'R  Jri.NPOlNT 
TTR  0( I , J>=-G( I , J) 

CAM  ID  (R.n,A,n,HRT.NRI6ID.NP0INT.M,MAX0,MAXR,‘'AX0) 

READ  (NTAPF4)  ( HR ( I . 1  ) . I  si . H ) 
no  MO  l  =  l.M 

DD  1411  JtMNRIGID 
dtl  G(  J,  I  )  :  URT(  J.  I  )/HR(  Id) 

WRMF  (NTAPr  1, '-»*>) 

CAM  HPRINT  (f;.NRIOin,M,MAXM.NTAPE.T) 

If  (  NPUNCM  )  3-1.1,  14S,  14^ 

M3  WRMF  (NTAPI  'i,:«MNRIOin.  M.  DCDZ.  NCARDS 
NCARD  =  NCARDS* 1 

CAM  MPUNCH  (  O.NRinin,M,!l,1,1,ncnZ,NAXO,NTApFU.NCAoO) 

WRMF  (NTAPf  1,?/ )HCD7,  NCARDS,  BCDZ.  NCARD 
NCAHDS=NCARD* 1 
M*>  CONTINIIf 

CAI  I  HMUl  ID  (DH..i.G,(1,HRT,NR{GID,NRIGID.M.MAX0,MAX0.MAXQ) 
no  MA  l  =  1.NR1fiFD 

Dd  MA  Jr)  .mode 

MA  HR(J,  I  )  =  n(  I  .  J) 

MU  no  3‘.(i  1  =  1  .npoint 

no  Mil  iri.MnnF 


4 0 II  _ 
'1  b  1  !| 

4  b? 


1b6!i 


4  60  1 
4610* 
*‘6?0 
^63  7 

464 
4  6‘>  i| 

4  6  A  ^ 
46  /  j 
4680 


4  6  9 IL 
4  7  II  T 
4  7  1  4* 
47?o 


-17  IT 
47bu 


4  76tj. 
477  ; 
4786“ 
4790 
4807 

4  8?0 


4  851*- 
484' 

4  851)* 

486U 

487;' 


4H8t» 
4  89  0 
4  90'^ 
490  I 
490.1 
4900. 

-‘'>1  I 

4  9 

493  0 
494"* 
‘>0  0-.. 
*>  0 1  II 

b  n 
-XKIf 
•>  II  4“r 
b  0  A  I) 
•)  0  6  S 
oil 


bll8u 
b  OP'*'* 

bin' 
bl  Ul 
b  1  ?  y. 

bl.1  I 

b  1  4  n  • 


75 


t 

c 

i; 

r 

I 


.l‘»n  F<J,  I  )sA(  I  ..n 

CAI  L  HHIILTD  ( F .  n  .  U.  n  .  A ,  H.  NPO  t NT  . H.  HAXR>  HAXR>  HAXR  ) 

IF  (  NnfLM  )  .^5?.3!>R,.’<b? 

T>?  IF  (  NRIGIR  )  •^54,.^5f|,:454 

■574  CALI  HHULTD  ( HR .  •) . HRT.  (I .  F . H» NRIOI D* H« HAXR .  HAXO* HAXR  ^ 

no  5'jft  1*1. hoof 
no  556  Jsl.M 

^56  A(J, 1  )  *  A( J. I )  ♦  F( J. I ) 

55A  URITF  (NTAPF5,94)  (  I,  A( I . I )» 1*1 »M) 

ROTO  991 

FND 

FORTRAN  LSTOII.DFCK 

HRrAH 

MHPAO 

MATRIX  RFAO  SUBROUTINE 

CAII  MRFAr)  (A.M.N.  IFORH.  IRON.  ITRA*  tOR6,T*HD,NTAPF2.NTAPE? 


) 


A  *  MATRIX  TO  REAR  IN 
M  *  NUMBER  OF  ROUS 
H  *  NUMBER  OF  COLUMNS 
IFORM  =  -1,  F0RHAT<12A6) 

*  ">  COLUMN  BINARY 


ITRA  *  n,  TRA  CARO  oFTtR  MATRIX 

..'A  CARO  *FTER  EACH  RON 
/OR  COLUMN'  ) 

lORG  *  ORIGIN  OF  FI:>ST  C.R.  CaRO 
T  «MDXN  temporary  GrLLS 


I 

I 


1 

I 

I 

I 


I 

J 

1 

I 


*  *1,  F0RHAT(6E12.tt)  MO*  DIMENStONEO  NIIMoER  OF  ROUS 
IRON  *  .0.  matrix  by  COLUMNS  IN  A 

so,  matrix  by  ROUS  NTAPE?  *  INPUT  TAPE 

NTAPE3  *  OUTPUT  TAPir 

SUBROUTINE  MRFAD  ( A . M. N* I  FORM. t ROU* I TRA , 1 0R6, T. HO* NtaPE? »NT APt3  ) 
OIMFNSION  A(l),  T<1) 

1  FORMAT  <6F1'.>,B) 

?  FORMAT  (IPAM 

5  roRMAT  (  / /  ^6H  THAIS  ALL  YOUR  DATA.  ) 

4  FORMAT  (4eiA.>0 

MN=Mn»N 
no  5  1*1, MN 

KDsOifl 
«>  A(  I  )  =  «.« 

If  (  IFORM  )  59,15,6 

6  IF  (  !ROU  )  «,/,M 

7  K5*1 
K4  =  N 
K4rMH 
K6=M-1 

K?s1 
GOTO  9 
n  XV=MN 

K.lsMD  ^ 

K4  =  H 

<•>=! 

K  6  =  li 

0  no  1 1  1  =  1 ,K4 

K1  *|»KS-K'>*  1 
If  (K6)  111.  11  1.10 


1  n 

K2SK1 4X6 

1  1  1 

If  (  If  ORM  -  1  ) 

54.110,109 

1  119 

REAli  INTAPF2.4) 

( A(L)  .L  *  K1 ,K2.K5) 

GO  TO  11 

1  1  n 

PfAl)  INTAPf'’ 

.  1  )  (  A(L),L«K1  ,K'J',K3) 

1 1 

CONI  INIIE 

OOTO  36 

)'» 

K  1  :tN 

K?  =  M 


51SU 

5160 

5170 

5180 

5190 

5200 

5210 

5220 

5230 

5240 

5260 


MOOlO 
MOOPO 
Mfl(j50 
M0040 
Mr.  050 
MOfifii 
M0078 
MOOBii 
MM.90 

Momu 
Mono 
M012U 
M0130 
M0140 
M0150 
M0170 
M0180 
H0190 
M0195 
MOpnu 
MO  21  0 
Mopj'a 
MO  2  5(1 
M0240 
M025m 
M0260 
MO  2  7  0 
M(/2Ru 
Mn29U 
MOSOU 
MO  31  i; 
M0320 
M0330 
Mn34li 
M0350 
MO  360 
M0370 
M038<i 
M039U 
MOAOI) 
M04n2 
M0404 
MOUIIb 
MOAIO 
MO  41 5 
M042U 
Mn430 
H0440 


\ 
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K.Ui 

IF  (  IOHO-1  )  16,1/. »7 

16  K  I  r  V 

17  If  r  IPOW  )  IH.IO.IH 
!>• 

Kl  =H 

lO  IF  (  ITRA  ) 

Klrl 
:79  Ksft 

HO  {>;5  I=1.K1 

K<-K+K3 
Kv  =  l 

HAll  BINRD  (T(K4).  KS,  1..  NTAPFl?.  NTAPE4  ) 

onto 

K  =  K4K/! 

IF  (  IRO«  ) 

?4  L  =  n 

IF  (  lORR-ll 
L=IORG-I 

,/6  no  ?/  1=1, R 
j=i*Mn-Hi) 
no  PV  K=1 . H 

j=j*i 

i.=i 

n  A( j»  =  T(i ) 

GOTO  16 

/f\  l.sn 

IF  (iORO~i)  r«n,.5n,?9 
i'Q  Ls|ORO-1 
.10  no  31  Ksl.N 
Js  K*MD-MI) 
no  .11  IsK.MN.N 
JsJ+1 
KlsL+1 
A(J)  =  T(K1) 

.16  RFTORN 

.10  MRI1F  (NrAPF1,.1) 

55TOP 

.19  RFAII  (NTAPr-^,9)(A<  I  ),  Is1  ,M) 

OOTO  16 
FNn 

%  FORTRAi  '  Simi.  DFCK 

r  niNRo 

SIlRROllTlNf  IMNRII  (  T  .  K  ,  I  ,  NTAPFl  .  NTAPb?  ) 
niHI  NSI ON  T( 1  ) 

RI  TURN 
FNn 

%  FORTRAN  I  S  mil.  IIFCK 

r  Rill  N 

SlIRROllTiNt  RIH.N  (NTAPF?,  NTAPb  1,  I  ) 

1  FORHAKRIIH 

I  ) 

5»  FORMA  I(  IHI  ) 

1  FORMA  I  ^  HIM  1 
Rf  Alt  (  NTAPI  •> ,  1  > 

GOTO  (4,0  ).  I 
4  WRMf  (NiAPt  1  ,  ■'  ) 

00  Til  6 

•>  WRl  IF  (NTAPf  1.  M 
6  MRMF  (  NTAPI  1 .  1  ) 

RFTORN 


I 


MM  46  I! 
MIMAl) 
Ml'47l) 
MMAHlI 
M(>4RII 
MObnii 
MM61  U 
MUbPll 
Mn63(i 
Mn64ll 
MOb^O 
Mn66i; 
MU  6 7  II 
MM^Ai) 
Hi>6Qij 
Mo'-  ^  , 

M  fi  ^  I. 
M  M  ^  ?  (I 

MM  Mu 
M0ft4(l 
MM65II 
Mnft6M 
M II 6  7 1) 
MIlAMil 
MIIMOU 
MO  711  u 
MO/l  u 


I 

I 

I 

I 

I 

I 

I 

I 


Mil/?i) 
Mil  7. ’HI 
MO  7  41) 
M076I) 
MO  7  6 1) 
Mn77lj 
MO/R(l 
MM79U 
MIIHIIU 
MI181  u 
MMM.li! 


I 

I 

I 


MO,04  IJ 
MII86U 
MM  06 11 


011  1  U  <<■ 
niipu 
on.iu  T 
(I  (1411  X 


II  M  I  •! 
II  M  V  (I 
(I  I)  .1 1) 
0  II  4  I- 
II II'.., 
0  It  /.  (I 
I)  (I  /  M 
OHM  II 
(I  null 
li  1  ll  II 
I1 11 II 
in  ?  n 


i 

I 

I 
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FNn 

FORTRAN  LSTOU.HFCK 
MMULTn 


I 

SURROIITINF  HHIII  rn  CA.Nl.B. 
niHFNSION  A(1  ).  B(1 )>  C(t> 

IC  =  1 

1 

I AsHCvK 

IH:><A*N 

X 

T 

ID?HA 

IHs-HC 

1  JsMC 

£ 

IF  (  N1  )  4.. 5, 4 

5 

IF  (  N2  )  /.O./ 

4 

IB 

1 

If»  =  /»ID 

IF  (  N2  )  5.6,5 

5 

\z=y 

I 

OOTO  7 

:! 

•  a 

6 

INsz-^lH 

■  > 

7 

IC  =  3 

I  As2«|A 

'e 

IJ  =  '>*1  J 

R 

no  16  1*1, M 

1 

INrsA 

no  11  Jsl.IA.IH 

ii 

ri  j)s(i.u 

1 

INxlNC 

DO  10  L*I.IR,ID 
lNsIN*l 

aiM 

1  n 

r( J)sC( J)« A(L )*R( IN) 

T 

4 

1 1 

INCsINC^HH 

INCsO 

-  k 

«» 

OOTO  (16,12.1‘'),IC 

1? 

no  14  J:I . 1  A. I J 

1 

IFsl »HA 

IF  sJ+HC 

-I* 

IN* INC 

no  1.5  L*IF.IB,IP 

i 

INslN^I 

I0sIN*HH 

“H 

1  3 

n(IF)*R(IF)*A(L)»H(IM) 
C(J)si;(J>-A(L)»R(  10) 

1  4 

INRsINC*H» 

'•I* 

1*^- 

omo  IH 

IF  =  |4Mn 

\ 

IF  *|4MA 

** 

00  17  JsIf.IA.IJ 

••  » 

IN*  ini: 
ri j)*o.o 

•  » 

nn  16  LriF.iR.in 

-# 

1 

1  6 

ins|N4| 

ri  J)sC< J) • AIL )*R( IN) 

A 

1  / 

iNn3iNn»HH 

1  H 

CONI  INIIF 

HI  TURN 

i. 

FNn 

FORTRAN  LSinil.llFCK 

r  MNVNSX 

r  SliRROUTINI  HNVRSX  (  AH ,  A  I  .  R « 0 » KS7 »  I  BOOFi) .  NOP ) 

|,  DIMfNSION  At  (A.  6).  R(A.^).  C(A,A) 

l(;nnFi»=n 


UldO 


601U 

A0?0 

0040 

oo?u 

AOAA 

0A7tl 

OURO 

0A90 

AJOO 

Ulltt 

OlPO 

ni3tt 

n)4u 

01A9 

0170 

OIBO 

0190 

0200 

0210 

0220 

0230 

0240 

02&0 

0260 

1)270 

02n0 

0290 
0300 
031  0 
0320 
03.50 
0340 
0350 
0360 
«370 
0  3R0 
0  39  0 
0400 
041  0 
0420 
0430 
0440 
0450 
0460 
0470 
0460 
0490 
0500 
0510 
0520 
0530 
0540 


ooto 

0020 

0040 
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I  f  (  NOP  )  j  I'  • .  I  01  . 1  np 

HM  -^Al  1  INVFRS  {AR.KS7.  inOvIFD) 

an  TO 

10/  nOMTlNIlP 

f)0  I  K=1  .KSZ 

no  1  i.si  .Ksz 

I  R(K,l  )  =  AR(K.I  ) 

N0=0 

CALI  lMVFRS<tt.KS7.Nn) 

IF  (NO)  ?.3,P 

r  RFAi.  matrix  not  SINOULAR 

C  MOLT  R#Al  STO.  C 

;5  no  4  K=i  ,Ksz 

no  4  L=1 .KSZ 
C(K,L  )5n.(t 

no  4  iri.KSZ  ■ 

4  C{K,L)=C(K,l)*R(K. 1 )»AI ( I ,L) 

C  MULT.  AI*C  ♦  AR  STO.  R 
no  ^  Ks1 .KSZ 
no  'J  L=1  .KSZ 
R(K.L)  =  AR(K,I.) 
no  H  in  ,KSZ 

^  H<K.l  )  =  B(K.I  ).A I (K. I  )*C(  I  ,L  ) 

NOsi.  ‘ 

CAM  INVFRSIR.KSZ.NO) 

IF  (NO) 

C  SFCONn  MATRIX  NOT  SINGULAR 

C  MUL-T.  -C*R  STO.  A I  ALSO  SET  ARsB 
7  no  H  kal.KSZ 

no  0  LSI .KS7 
AKK.DsO.n 
AR(k.L)sB(K.I.  ) 

no  H  isi ,Ksz 

n  A|(K,L)sAI(K,l.)-C(K,I)»B(I.L) 

GO  TO  5»n 

C  RFAL  MATRIX  OR  SFCONI)  MATRIX  SINGULAR 
no  9  K=1.KSZ 
no  y  L=i .KSZ 
9  R(K,L)=AI (K.L) 

NOs(. 

CAM  INVIRSIH.KSZ.NO) 

IF  (NO)  1  ti,  1  1  .  1  II 
r  IMAO.  NOT  SINOIII  AP 
C  Mill  T.  lUAi:  Sin.  C 

II  no  I  9  K  =  1 , KS7  ■ 
no  :P  I • 1 ,KS7 

C  (  K  .  I.  )  =  II .  (1 
1)0  IP  I=1.KS7 

IP  C(K.I  )=C(K.I  )»R(K. I )*AR(I,L) 

C  Mill  T.  AR»r.  A I  STO  rt 
no  1 ^  K=1 ,KS7 
no  I  .A  L  =  1  .KS7 
n(K.l  )  =  AI (K.l  ) 
no  M  1=1 ,KS7 

11  n(K,|  )  =  «(K.|  MARIK.  I  )*r(  ML) 

NO=i. 

CAM  INVfRSin.KSZ.NO) 

II  (NO)  1il.M>.lil 

r.  IMIRO  MATRIX  NOT  SINOUIAR 
C  Mill  r  -(  •U  STO  AR  ALSO  SET  AI=-R 
!•>  1)0  I A  K  =  1  ,KS7 


THY  IMAG.  ROUT'' 


II  0  >5  U 
niiAiil 
(I  n  /  0 1 
0  0  8  0 

li 

01  (Ml  I 

0 1  1 1I 

'>1.1  II I 

III  4U* 
0  1  5  I) 
ni  Aiil 
0  1  7  U I 
n  1  ftu 
OIQIM 
npf.ii 

?  1 0 ' 
It  ’/  ?  ii 
^  1  ■! 

<1  /  4  ,)  ■ 

'  \ 
IM.  !•'' 

0  P  h  IJ 
OPOO. 
n  1  (I  II 
(1.51  ir 

n.ipo 
UA.IU’ 
(1.54  0, 
(i3'>ir 
(I  .^All, 
0  370 
0  38  ir 
0  3911 
0  Ml  U ' 
1141  I), 
0  4  ?  n 
f|43u 

0  4  4  -1 
4  S  ' 

('  4  7  .1 
(MHO 
0  40(1 
0  ►)  II  1) 
OHI  U 
fM^Pl) 
OH  311 
I)  4  II 
(iHhil 
IIHAO 
0  S  7  II 
II  b  B  U 
ItbOO 
0  60(1 
0  M  (I 
06P0 
0  6311 
064(1 


•'9 


no  16  t  si .KS7 
AR(K ,L)s« 

AI(K,L)s-B(K.|.) 
no  16  1=1, KS7 

6  AK(K,L)sAR(K,i )»C(K, I  )«R( I .L) 
no  TO  ?» 

0  1  BOOF  0*1 

n  RETURN 
ENO 

FORTRAN  LSTOtl.OECK 
INVFRS 

nUNROIITiNE  INVFRS  (A,N.|  GOOF  R) 
niMFNSION  A(6.6).  1.(6),  H(6) 

CAM  OVERFI  (KiHt'iFX) 

no  Toisnii.hniD.Kiim.Fx 

00  CAI  I  OVFRFL(Klin«FX) 

GO  TO(6fll,‘io  I  ),K(iniiFX 
01  CAI.I  lIVCHK  (KHfUlFX) 

GO  TO('»n^,‘>'i;') .K'iniiFX 
•  0?  IGOOFOsO 

.SEARCH  FOR  LARGEST  ELEHENT 

no  HO  Ksl ,N 

i.(K)sK 

H(K )sK 

niGAsA(X*K) 

no  vn  fsK.N 
no  vn  jsK.N 

IFT  ARS(RIGA)-ARS(A<  I  .  J ) ) )  10 »  ?ii .  ^0 
lO  niGAsA(I*J) 

L(K)sl 
H{K  )*J 

>0  CONTINIIF 

INTERCHANGE  ROUS 
JRORs|.{K  ) 

IFd  (K)-K)3‘j.  I*),?*? 

?»>  no  .10  1=1  ,N 
HOI  n  =  -A(K,  1  ) 

A(K.  1  )sA( JROU. I ) 

TO  AC.IROW,  I  )*HOU) 

interchange  COI.UHNS 
.IS  ICO|sN(K) 

IF(M(K)-K)A*>,4»;..1/ 

.17  no  4  0  J  =  1  ,N 
MOI  l»s-A(  J,K  ) 

AI.I.K  )sA(  J,  ICOI  ! 

4  0  A(.l.  lent  )sHUI  11 

niVlOE  nOlllHN  RY  MINUS  PIVOT 

45  no  • S  ic=l ,N 

46  If  (  IC-K  )SM,‘>r, ,»,« 

SO  A( lr.K)sA( IC.K)/(>A(K.K)) 

55  CONIINUE 
RfOllCF  MATRIX 
no  ».*>  1=1  ,N 
no  /.'*  J=i  ,N 

56  If  (  I-K  )s;,6‘>.*>7 

57  IF(.)-K)M>,6‘’.6II 

60  A( I , J)=A( 1,K)*A(K, J)^A( I, J> 

6S  CONTINUE 

niVIOF.  ROW  OY  PIVOT 
no  /s  JRS1,N 
6S  If (JR-K)70,/S,7n 


0660 
U67U 
06R0 
0690 
0700 
0710 
0  720 
0730 
0740 


0010 
II 020 
0040 
0090 
0060 
0070 
0060 
0090 
0100 
‘J 1 1  'I 
f  1 '  II 

0130 
0140 
0190 
0160 
0170 
niGO 
0190 
0200 
0210 
0220 
0250 
0240 
«?50 
0260 
0270 
02H0 
n?90 
".TOO 
f‘31  (• 
0320 
0330 
0340 
0390 
0360 
.  0370 
0380 
0390 
0400 
0410 
0420 
0430 
0440 
0490 
0460 
H470 
0  480 
0490 
0900 
0910 
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\ 


t 


i 


i 


I 

« 

I 

V' 

5 


I 


I 

B 

I 


R 


n  A'X, JR)sA(K,JP>/A{K,K) 

75  nONHNtlF 

r  roNTiNiJFT)  pkomiCT  or  Pfvnrs 
r.  PFPIACF  PIVOT  HT  RFriPROCAL 
A(H,K  )sl.«/A(K.K) 
r.  COHTINUF.  OOMPLFTF  OPFRATION 
HU  CONTINUE 

CALL  nVCHK  (K'lOOrX) 

00  T0<5tft,5n.>),K'in(iFX 
5n.T  CAI I  OVERFLIKiinoFX) 

GO  T0(5l0.bft4  )  .K'ln'iFX 

504  CALI.  OVFRFLiKlinOFX) 

00  TO<510.,5n5).K''nitFX 
C  final  row  ANO  rOLIlHN  INTERCHANGE 

505  KsN 


100  Ks(K-l) 

IF(K)15n.l5(i.1!M 

Is|  <K) 

|F(  |-K)l?0.15»ii.lnT 

no  110  Jal#N 
HOinaA(J.K) 

A(.I.K  )s-A(  J.  I  ) 

A(.i.  I  )sHOLD 
.JsHIK) 

IF(J-K)inu.lrt'i,l?5 

no  130  1*1.  N 

HOLIlsAIK.  I ) 

A(K,nx-A(J,n 
AU.DxHOLn 
00  TO  inn 
RETURN 

I  doiiF  n=i  . 

.00  TO  i»>n 
ENn 

T  LINK  PAR?TT.PAR1TT 

f  FORTRAN  LSTOU.DFCK 

C  PART? 

OUPROUTINE  PART? 

C  PARti.,..?  VIBRATION  AND  FI UTTR  ANALYSIS  BY  A  COLLOCAllOw  METHOD, 
n  I  mens  I  Of  iT(?1  S),VFLCTY(?0).NSIZES(?0)»nM(6.6).RHO('»H) 

COMMON/I  I/  U(4V,196)»0UFSS(49,2).H(4V»5n),EIGf5n), 

1  |eMP<?/3/);.  NAKSR(^5),  NAK0R(?5),  NlTER(/5), 

?>  OMFGAI^b),  nAMP<25),  VELC(?,5) 

C  THF  FOIIOWINO  SIATFMFNT(S)  HAVF  BEEN  HANUFACTUREO  BY  THE  TRANSLATOR  TO 


10  1 


105 


11  U 
1?ll 

1-?5 


1.1  II 

150 
51  n 


nONPFNSATr  FOR 
COMMON  I T 
d'OUf  VAl  FNCE 
I 
? 

.1 

4 

5 
5 
/ 

A 

V 

Foil  I  VALENCE 

1 

? 

.1 

4 


THE  FACT  THAT  EQUIVALENCE.  DOES  NOT  REORnER  COMMON— 


( IT( I ). ISXST). 

( I r<5l ) .NTAPE1 ) 


(  IT(54) 
<  IT(fi7) 


NTAPF.4  ) 
NTAPE/) 


(IT</ii).NSUR). 

( |T( /1),NAFR0)» 
(|T(/A),NOnES)/ 
<  n(/9).MAXR). 


( IT(?J )» ISM) ,  ( lT(4l),RHn) . 

.  ( IT<A2),NTAPF?),  <IT<6«).NTAPE1) 
.  ( IT(A!»),NTAPF5),  ( I  T(6*  ) »  NTAPE6  ) 
.  ( IT(5B),NTAPFH),  (|T(6o),NTAPE9) 
( IT(7i ),NRIOID),  (1T(/?).HREF). 

< IT(/4),NFUS).  ( IT(7b)»MUFNS)» 

( IT</7),NP0|NT).  ( IT<7«i.NPUNCH), 
(  IT(Ali),MAXO>,  (IT(81)*MaXS). 


(  IT(t«?).Nr).  (  IT(Hd),NSURFS).  (IT(84)*Np?). 

(  I  T(H*»)  ,HR>,  (IT(1II5),S) 

{ T T( l?/),VFLCTY).  ( I T ( 14 / ) , N  I A  PE  0  ) ,  ( I T M 4h  ) » NRHO ) 
<?T(i4«n.NlTRSP).  {  iT(15U).NlTKl)p),  (  I  T  / 1 5  1  ) ,  EPSP  ) 
(  I  T(  i5;o,FPnP).  ( IT(  l»iJ)»PICON) .  (  ITdbO.NVEL). 

(  I  T(  I55).NCAR0S),  (IT(15C»).AiTKEN).  (  1 1  ab  /  ) .  N60 )  . 
(IT(15H),FK).  (  1 1(1*19  ),S0RU0N),  (  IT(IAIm,NSIZES)» 


Ob?  II 
Ub.lU 
0  5411 
05511 
«5<S!) 
f-5  7  0 
ObHO 
(1590 
(1601) 
061  0 
(I6?0 
06  1') 
0  64  0 
l)65u 
06f  u 
ft  6  7  0 
n6H0 
11690 
(I  7  0  I) 
071  0 

(17  10 
0  74  0 
(I  /  5u 
(17611 
0/70 

(I7AIJ 

(1790 
(1800 
(181  II 
(18?  0 
OHIO 
1)840 
0  85  0 


I 

I 

I 

I 

I 

I 

I 

I 

I 

I 


onio 

I)  0  ?  II 


r 

L 


0  0  4  0 

0  0  /  (I  I 
0080 
(1119(1  *»■ 
0150  1 
0  1 6  (I 
017  0-^ 
018  0  ,i 
0190 
(I  ?  0  0 

n?j  0  T 

o??o  i 

()?.1U 
fl?40 
fi?5  0  J 
ll?60 
l)?70 

(i?8o  r 

(l?On  I. 
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6  (lT(l«n),nM).(IT<'^17).AITKtO).  (IT(218).KPa8T)  031 

900  FORMAT  (IHI  ^HX,  IVH  OUTPUT  DATA  tit  IH  11X»  fLUTlPR  03?« 

1  ANALYSIS  RY  A  COLLOCATION  HEiHOD*  USING  AERODYNAHIC  033(i 

t  ^4H  INFLUENCE  COEFFICIENTS  //  IHO  14X*  11m  DENSITY  s  0340 

)  lE?n.o.  Sx.  20H  REDUCED  VELOCITY  =  IE20.««,  //  IHO  33X*  0380 

4  116.  39H  RIGID  RODY  DEGREES  OF  FREEDOM  /////  >  03Att 

9111  FORMAT  (1H1  VOX.  8H  NODE  7X.  1PH  ONEGA  ICPS)  lUX.  HH  DAHPINO  037D 

1  4X.  17H  VELOCITY  (KNOTS)  ///  (  IH  2-X.  114,  «3AI 

2  3F'^ii.«)  )  03R0 

9)19  FORMAT  (  214.  A9X.  1AA,  1  14  )  0400 

903  FORMAT  (  IHf  3.1X.  V3H  PUNCHED  CARDS  NUHRERS  1A6,  If 4,  AH  THRU  0410 

1  1AA.  lit  )  II4?0 

9i|4  format  (1H1  4f)X.  1AH  DYNAHIC  MATRIX  )  *1430 

9iiS  format  (.51M  MEXIRLE  MODE  SHAPES.  SURFACE  1  14  )  0440 

9i)A  FORMAT  (8M  MODE  llA.  3?H.  GIVES  AN  IMAGINARY  rREOtlL*»CY.  )  0480 

OATA  QOlinCT/O90.^fl.1044i|fl08/ 

RCD7  aOnOliCT  0470 

NVElaNVEL  0480 

HAXP*49  1)490 

KaNC»NPOINT  0500 

no  290  Isi.NPOINT  0510 

DO  290  Jsl.K  0520 

9Vn  U(I,J)a(i.O  ’0530 

REWIND  NTAPE9  454# 

no  314  IRHOsl.NRHO  '(55» 

MODEsNOOES  »  •J550 

READ  (NTAPE9)N8I2F,  NSIZE2,  (  (U(  1 .  J )  *  I  ,  NSl  ZE  ) ,  f'570 

1  Ja1.NSlZE2)  0580 

IF  (  NPUNCH  )  302,104,(04  0590 

ll.i  WRITE  (NTAPE(,904  )  0600 

CAI I  MPRINT  (U.NPOINT.K.MAXP.NTAPEI)  0610 

304  WRITE  (NTAPE(.2ijn)RH0(IRH0).  VELCTYI MVEL  ) ,  NRIGID  0620 

no  305  1*1, MODE  0630 

NAKDR(I)sO  0640 

305  NAKSR(I)«li  0650 

CAI  I  miters  (U. guess  .0  .NPOiNTrMODE  ,MAXP  *NC  .EPSP  *  0660 

1  EPDP  .NAKSR  ,NAKUR  .  N1 TRSP*  MI  TRDp,  A I  TKE>‘ ,  A ITKEO,  0670 

2  IR  .TEMP  .H  .EIO  .NITER  . NTAPEn , NT APE3 )  0680 

IF  (  NAERO  )  3055.3054,3055  0690 

3054  WRITE  (NTAPE4)  NPQINT.  MODE.  ( (M( I , J) . I *) . NPOI NT ) , J*) . MODE >  0700 

WRITE  (NTAPE4)  ( E IG( I ) , I *1 . MODE ) . NPOIHT . MODE ,NPOI NT . HODE.MODE  0710 

30', 5  CONTINUE  0720 

MOm  S2  s  NCtMOOF  0730 

no  310  1*1 .HODESP.NC  0740 

Ks|/NC  INF  -1  0750 

If  (hlOd))  3051 ,3052,3052  0  761) 

3061  WRITE  (NTAPf  <,2i)6)K  0770 

OMFOAlK  1*0.0  07An 

ftOTfl  3053  0790 

10.2  nMlGAlKl*  SORCON  /  (  SORT(  FK^EIGd))  )  OROU 

3II‘»3  tiriTO  (306,3110)  .NO  081 U 

1(16  DAMRIKlaO.n  0820 

VHC(K)aO.O  0830 

GOTO  110  0P40 

(HR  DAMPIK)*  FlOdd)  /  FlOd)  0850 

VdfMK)*  PICnN*nMEOAIK)*RRFr  •  VEI.CTY(NVEL  1  0860 

on  CONTINUE  08  70 

WRITE  (NTAPI  (.901  )(K,  OMEGAIK).  OAMP(K).  VELC(K)  ,  0880 

1  K*T . MODE  )  0H9U 

IF  (  NPUNCH  )  312,114,31?  0901) 
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Hz'  K1=1 

MRfTF  (NTAPf  n,9o?  )NS1ZFS{  n,  HOOFS’#  BCI)7,  NCARDS 
NCROS  =  NCAKOSM 

CAM  MPliNCH  »nMEftA.MOl)FS7.1»n,»,t,flCDZ.HA)(P.NTAPEo.N'’HnS) 

MRMF  (NTAPH,:»ll3)BrDZ,  NCARHS.  RCOZ.  NCROS 
NCARDSsNCRDS^  t 
no  ISUR  =t,NSllR 

IF  I  ISUR-1  )  Tl?l. U’l 

.U?1  WRMF  (NTAPFn,?(|7)NslZFS(  ISUR),  MOOFSi:  >  hCOZ  -  MCARDS 
NCROSsNCARDS* 1 
00  10  .5123 

51;'?  MCROSsMCARDS 

51?.1  CAM  MPUNCH  IH(K1  .1  ).MSIZeS<  ISUR'.MOnSSSU.l  .1  >UCDZ»*‘AXP,NTAPbfi# 
1  NCROS  ) 

MRMF  (NIAPI;  5, 9(1*4 )  ISUR 
WRITF  (NTAPFI.’U.I  )flCOZ .  NCARDS  .  RCOZ  .  Rr 
NCAHOSsNCROSH 
11.5  KIsKI  fNSiZFSI  ISUR) 

114  CONTINUE 

REWIND  NTAPFO 
Non  *  9H765 
KPAPT  s  1 
RETURN 
FNn 

F  f,..'TRAN  LSTnii.DFCK 

r  N0RM7 

SURPOUTINE  NORM  < A . H . N» C . I NOFX . MAXR, NC, NP ) 

OIMINSfoN  AM),  H(l),  C(l).  T(?) 


IF  (  INDEX  1 
nt<'  00  TO  (110»14»>).NP 
11«  00  TO  (5'«<1.1?«),NC 
i?n  oio  CM  )**?«r(? )«#? 

00  TO  (SOO 

OAM  DNORH  -  A . R, N . C . HAXR. NC* T ) 

liO  TO  7«fl 
’:'0  !  NOr  X  =  1 

NSTAR^^NP^t 

NSTnPsN»NP 

KsNp*NAXR 

IF  (  NSTARI-NSTOP  ) 

’n*>  00  TO  (  la ,  4'in  ) .  NC 
71  n  OICrAHS  A(1  )  ) 

1)0  ?.51I  -NSTART  NSTOP.NP 

II  (  11  |0-AUS{  A  (  I  )  )  ) 

?-'n  IND|-X=I 

n!n  =  AiiS(A(  I )  1 
?.5ft  OONTINUF 

00  in  /loii 

IIII!  M|n  =  AM  )«•'»  «  A<K*  1  )*•? 

no  :mo  i=kstart.nstop,np 

Mf*  I 

■f  (  B?0-(  A(  I  )**.'*A(  J)«*V  )  ) 
1711  mo  s  AM)**'^  ♦  A(J)»#? 

INDFX  =  I 
550  CONTINIU 
4iin  .)  =  NP»NC 

f  =  IM)f  X*(NC-i  )»K*  (NP-i  ) 

C(  .))  =  A(  I  ) 

.»  =  .)-  <  NP-1  ) 

I  =  I-(NP-1  ) 

C  (  J  )  =  A  (  I  ) 


1  0  R  ,  ?  II  (1 »  ?  I)  II 


?M'>.2li'>»40ll 


?20,?i;n»?.5li 


>5?n,32n,  130 


II V 1  fl 

It9?|| 

UP  5|| 

1194(1 

nviil 

n9M|[ 

0  y  7  ii 

II9RIU 

II99|| 

1  (I  n  if 

1  0  1  II 

1 11  ?  (I 

1  11.5  J 
1040 
1  (>b(l 

1  (I  f)  J 

10  7  0 
)  (I  ft  u 

1  npul 
n  nu! 
1  n  u 
1  1701 

1  1 .1  o] 

1  1  4i< 


MTR4fiOO^* 

MTR4001JI4 

HI  0400061 

HTR400t)/i 

HTR400nf) 

HlR4onn9i 

HTR4flniu 

HTR4nnl2^ 

HTR4nni 3 

HTR40ni^] 

HTR40016J 

HTR4noi ; 

HT'/40()18- 

HTR4fl|i  1  9  ^ 

HTR'i-ifi’ii 

H 1  R  4  *  0  7  e 

H  i  9  4  0  (1 7  „ 
HTR4>;.'7!>' 
MTR4(1 07u 
HTRIOO?'/ 
HTR4(I  ll?a 
HTR4!i  ()?y 
HTR4U Oil 
MTR4(l  0:57 
MTR4''  0.1.5 
HTR4n()34 
HTR4  0  0.1b 
MTR40n3fi 
HTR4(in3H 

HTR40  (I4II 
HTR4(Ut41 
HTR4n()42 
HTR4no4 J 
MTR4()  044 


U  T  n  4  r* 


4  r 


^.yvirrinr  *Wi,’>’'M(rren»E«?*S!! 


I 

I 

I 

I 

I 

I 

I 

I 

I 
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T; 

ri 

li 


l  =  INni.-X*<NP  1  ) 
r(NP)sA(I) 
n(1 )SA( INDEX) 

DO  10  (41(1, MO), NP 

4in  00  TO  (5l)fl,Aiin),NC 
5.)n  no  '»in  i»i.n 
sto  P( I )sA( I  )/c;i ) 

00  TO  VOO 

AflO  olfi  s  c(l  )*»’*n(P>*»? 
no  Ain  lal.N 
JrMMAXR 

T=  A( I )*C(1 ) »A< J)*C(P) 
niJ)  5  (  A{  J)*C(1  )-A(  1  )*C(?*)  )  /  816 

AiO  R(|)s  T/Bin 
700  iwnpx  *  INDFX/NP  ♦NP-l 
RETURN 
FMn 

FORTRAN  LSIOM.DECK 
DNOPNZ 

OUOROIITINP  liNDRH  ( A  » H.  N. C .  NAXR* NC . T  ) 
noUHIE  PREOIOION  A(1).  8(1 ) .C(1  )*T(1  ) 

410  00  10  (Aoil.Aon  )  ,NC 
SOO  no  OVO  Irl.N 
•>10  R(  I  )sA(  1  )/n(  1  ) 

00  TO  701) 

ADO  Riosr(i)*C(n  ♦  o(P)»c(?) 
no  Atn  isi.N 
JsMHAXR 

Ts  A( J  >*C< 1 )«A( J)«R<P) 

R(J)  a  (  A( J)*C(1  )-A( I  )«C(9)  )  /  big 

AtO  R( I  )a  T/8i0 
700  RETURN 
FND 

FORTRAN  LSTOII.DFCK 

PQHS 

SUBROUTINE  PUN  ( LHORN, LMRDt » I.MRUP*HN> HNl  > Hi  , H2* N. Nr  ) 
liMENSION  |.MRnN(t),  LMROKDr  IH8D2(1)»  HN(1),  HN1  ( ■  ) » 
1  H?(i),  A(?) 

noURLF  PRECISION  LHRON*  (.NBRi^  I.NBO?,  HN*  HNl,  HI*  u?, 
on  10  (RilO.  1  nr  ),NC 

1 1  0  I  a  J  ,  U 
K::UN 

A<1)  a  LMUn^l 1 )*HN( I )-LHRDN{?)»HN(K) 

A(9)  a  LNHnN(n*HN(K)«LHBDN(V)«HN()) 

MKI)  a  LMRIPd  )»HNl(|)-lNRO'MP)«MNt(K)-A(l  ) 

HKK)  «  IMODMi  )»MNM  )*LNRI)?(9)»HN1CI)-A(R) 

H,(l)  =  A(  I  )-LHHni  ( I  /<HN1  ( I  )aHBni(2)»HN1  (K) 
nu  H>'(K)  s  A('')>LMHni(t  )*HN1(K)-LHBD1(?)*HN1(1) 

RE  TURN 

9  nil  nu  RIO  lal.N 

AM  )a  IHHDNd  )»HN(t) 

HI  (  I )  s  LHRnvd  )»HN1  { I  )-Ad  ) 

RIO  HRd)  s  A(M>LHHni*HNl(|) 

RETURN 

END 

lORTRAN  LSinu.DECK 
A  I TKNS 

SUBROUTINE  MIKNS  ( HN ,  f!N1  ,  HN^ ,  HNFN  *  R*  N  *  MAXR ,  NC  *  NP*  I R  > 
niMINSION  MNd).  HNId),  HNRd),  HNFH(I)*  A(2).  B(4«, 

IRall  ! 

l.aNP*N 


HTR4n04A 

HTR4004/ 

NTR40n4B 

NtR400*»l) 

HtMl)85l 

HtR40095 

I4TR40854 

HfR400***> 

NTR40057 

HTR40098 

HTR400Aa 

HtR40061 

NtR40062 

NTR4006i 

HTR40flA4 

NTRROOftiJ 

NTR40067 


HTRAOOAO 
HTR40071 
HTR4nn7? 
81840074 
HtR4flM7i 
MTR40076 
NTR40070 
HTR4ft079 
HTR40080 
>^TR40081  I 
Hf84«082  1 
.NtR4n083 
NfR48084 
HTR4n086 


HTR4noR9 
Hid).  HTR4nu9l 

NTR4nn92 
A  NTR400R.j 

HTRAOORi 
NTR40097 
HTR4009V 
HTR4«J01 
PTR40102 
NTRAOIOA 
HTKAOlOtJ 
MTR4nin7 
NTR4ntOb 
NTR4niio 
HTR4ftll2 
NTR40114 
HTR4ail6 
HTR4011/ 
HTR4nil9 
HTR4ni21 


HTR4ni2i 
U(2).  0(2)  NTR4i)12b 
HTR4ni27 
HTR4U12y 


\ 
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I 


I 


f 

■| 


) 

< 

? 

i 

t 


{ 


i 


i  f 
■  h 

i 


CC 1  )  =  R**? 

•  no  1 1 .1  I  =  1 .  t  . 

p  ( 1 )  i  II 

n<  ? )  =  0 

no  1  no  Jr  t ,  Mr 

Ks( J-l ).HAXR« I 

Hd)  s  tld)  ♦  (HNHK  )-HN7(K  >  >#«i' 

1"**  H(?)  s  IH  V)  »  (HN(K)-HNl(K)  )»»9 

If  (  n( 1 )  )  1  U5» 1 1 n, 1 

rONTfNUE 

If  (ncd/njt )  -  cd))  iin.un,pii' 

110  CONTINUE 

00  TO  (.Ton  ,‘^rifi )  ,K|P 

?iin  CAI  I  AIT.KND  (HN.HNJ  .HN2,HNEH,N,MAXR.Nn.  A,B,C,D) 

00  10  7011 

ion  no  6(10  1*1.  N 

00  TO  (4'in  ,5(10  ),NC 
4(in  0(1)  s  HN(  n-?,»HNl(  I  )*HN9(  (  ) 

* ^  (  0(1  )  )  .  410 ,611(1 , 4 1  II 

4in  HNFH(I)  a  HN2(|)  -  (  ( HNl  (  1  )*HN2(  I  ) )  ••/!  /(Ml)  ) 

00  TO  600 
Sun  A ( 1  )  a  0 . 
rd  )  a  2. 

0(1).*  0 . 

no  510  .1*1.2 

.  N  =  !♦( J-1  )»MAXR 
B(J)  a  HN(K)-.>.#HN1  (K)*MN2(K) 

C  a  C*(HN.l  (K)-HN2(K)  ) 

*  a  -<HNMK)  -  HN2(K))*«9  -A 
SIO  n  a  R(J)*»t>  4.0 

■  ■  -W ..  (  n  )  *>20,600,520 

5/0  MNFW(I)  a  HN/(()  -  (  0  ( I  )  •  A  ( 1  )  ♦»(  2  )  aCC  1  )  )  /  11(1) 

,HNFM(K)  a  HN2(K)  -  ( R  d  )  *0  d  )  oRl  2  )*A  (  1  )  )  /  Dd  ) 

5iifl  CONTINUE 
7  10  IP  a  1 
O.jO  RETURN 
FNO 

f  FORTRAN  LSTOll.nECK 

C  AITKND 

‘JURPOUTlf'"  AITKNO  (  HN,  HN1  ,  HN'/ ,  HNFN,  N.,  MAXR ,  NC ,  A ,  H ,  C  .  0 ) 
niHFNSin'.  UN(i).  HNId),  HN^(  I  ).HNEMd  ),  Ad),  B(2)  C(') 

noUHl  E  1  'ECISION  HN,  HN1  ,  HN/.  HNFN,  A,  H,  C,  U 

inn  no  6(10  1 :  I .  N 

on  TO  ( '•i  "II  ,  50 -I  ) ,  fn(j 

0(1)  r  HN(  I  )-2,sHN)  (  I  )4HN2(  I  ) 

IF  (  AOSfCd))  -.iKiiiniinioii  )  6uo,6iin,4tii 
4l«  HNFWd  )aHN'’(  f  )  -  (  (HN1  (  I  )-HN2(  I  )  )«(HN1  (  |  )-HN2(  I  )  )  /C(l) 
on  TO  6110 
50  0  A(1  )  a  (I, 

■  Od  )  a  2, 

0(1)  s  0 , 
no  510  J=|.2 

K  =  I ♦ ( J-1 )*MaXR 
0(J)  =  HN(K  )->',«HN1  (K  )*HN/(K) 

0  =  C*(MN1 (K)-HN2(K) ) 

A  a  -(  IHNI (K)-MN/(K) )»(HN1 (K)-HN2(ld  )  )  -  A 

510  ii  =  B(J)»H(J)  4  0 

(  D  )  S2n,6iin,*./ii 

5/0  HNFW(I)  =  HN-  ( I )  -  ( R (1  )  # A (1  )  aR ( 2  )  aC ( J  )  )  /  0(1) 

HNFM(K)  a  HN/(K)  -  (  R  ( 1  )  •€  (1  ) -R  (  2  )  a  A  (1  )  )  /  I)(1) 

600  rONTINUE 


HTR4-;l  J(j 
HTR4n) 1  " 
HTR4ni.5.1«, 
MTR4ni  .TA 
HTR10136'r 
>•10401 
HTR4i)14ij 
HTR4fH4l_ 

^104(11  4i  I 

MTR4()144^ 
HTP4n J  45 
HTP4n) 46  T 
MTR401 4a  4 
HTR4II 5  49 
MTP4015(1  y 
HTR4015/  I 
HTR4ni54 
HTR4CI155 
HTR40156  T 
HTR4(j15»<  * 
NTR4'H59 
1  *.  1 

KTP4‘i3*  i  J 
HTR4lllftj 
NTR4mft4  ^ 
HTR4(I166  ! 
HTR4ni6/  ■*> 
HTR4ni A8 
HTR4ni6y  r 
MTR4IM  70  i; 
HTR4I1]  72 
MTR4ni7j 
HTR4lt174  ; 
MTR4itl76 
M1R4ni 76 
MTR4naaii  T 
HTR4ni8/  ^ 


HTR40lfl4  ; 
0(1)  **jp4n ] Pft  “* 
HTR4018/  _ 
MTR4r,i  ;^g 
NT R4 01 91  « 

MTR4U102 

H7R4II19.)  *• 

NTR4nn)5 
MTR4(1 1  96 
HiR4ni9a  -• 

MTR4ni09  * 

HTR4fl2(>o  " 
MTR4()2ii  1 
HTR40/n^  7 
NTR4n2li4  1 
HTR4  0  2115 
MTR402n6  »** 
HTR4(l2n/  ■ 
MTR402fl9  " 
NTR4'l21ll  _ 
NTR4l)?li 
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( 111  n ,  9  lift ) ,  nr 

nABS(  I  AHRHAd  )-LAHBOAr>)  )  )/(1 

r)ARS(Pos<7)/POSrn)-i.)  /  fc? 

IIAHS(PaS(4)/P0S(  d)-1  .  )  /  k? 
(Rt  .  -*.1 
(R2A,V,1 ,?,A) 


/lift  IR  s  1 

Rtin  RtTIlRN  ■ 

CNR 

rORTRAN  LSTOU.OFCK 
l.fRISS 

<;(IPROlJTINE  IFRIS  (LAHHnA.POS.Ei  .P2.NC»  IK) 

DIHFNSION  LAHRnAd).  P0S<1) 
noilRLE  PRECIRinN  LAHHnA.  POS*  Ri  »  RPA*  R2b 
fRsli 

no  TO 

l.jfl  R1  s  { 

R9A  s  RAHSI 
RVA  s  OAHSC 
CAII  HPRINT 
CAII  HPRINT 

CAM  HPRINT  (R2R.R.1 .2.6) 
d"  rr  (  Rl-RPA  ) 
dn  IF  (  R1-RPB  ) 

Mn  IRsRC 
14n  RFTIIRN 

2iin  R1  =  nSQRTI  (LAHRHAd  )*LAHRDA(  A))»<LAH'BUA(1  )<LAHBDA( 
1  (?)-LAHRnA(4))»(i.AHBDA<2)-i.ANBDA(4>)  )  /  FI 

P2A  s  0ABS(  nsORT<  < ( POST  5 )*POSd )*POS( ? )*POS( 4 ) )•{ P 
1  P0S<'»)*P0S(4) )  ♦  (POST  d*P0S(?)-P0$(4)«P0Sd 

V  P0S<’)-P0S(4)*P0S(1 ))  )  /  TPOSd )»PUS(1 )4POS 

3  -1.  )  /  F2 

R2R  s  RARST  nSQRTT  (  ( POST  7 )«POS(^ IMPOST  6 )»P0ST8 > )  / 

1  (5)  4  POS<A)»POS(6))  )•(  TPOS(7)»PO$(V)4POS(6 

2  (PUS(H)«POS<H)  ♦  P0ST6)«P0S(6))  )  «  (  (POS(7) 

A  H)*POS(^))  /  (PQST6)«P0S(5)  ♦  P0S(6)*P0S(6)) 

4  tPUSd.)  -  POS(H)*PQST‘i))  !  TPOS(t>)*PuST^)  ♦  P 

»>  )  )  -1  .  )  /  F? 

dft 


t2A.l4n.l4i) 
liSA.I  40.1411 


HTR4021‘> 

HTP40217 

HTR4AP19 


MTR40272 
HTR40P74 
HTR4n?2b 
HtR40??7 
HTR40779 
HTR40731 
HTR4ft732 
HTR4073J 
HTR407.<4 
HTR4ll?3‘> 
HTR4n?36 
MTR40238 
HTR4n?40 
NTR4n?42 
NTR4A?44 
1))  (LAHRRA  PTR4n?4.) 

P|ft4<>?4 ; 

♦  HTR4A?4»j 

HTR40?49 
HTR4II750 
HtR4n?5l 
MTR407»).$ 
HTR40?54 


>ST3)*P0S(i  ) 
d»TP0ST3)* 
)2)*POS(?))  ) 


no  TO 
FNR 

FORTRAN 


{P0S(6)«P08 
)»PUSTH>)  / 

»P0ST6)  -  PUS(MTR4n26)> 
I  •  (  (POST/)  HTR402'>6 
^S(6>*P0S(6>)  htraapn/ 
HTR4n?6o 
HTR4n?6u 
MTR40762 


LSTOll.  OFCK 


C  MAO? 


SIIRKOUTINE  HAAS  ( A.B.C.NStZF.NC) 
nOOBlE  PRECISION  Ad).  B(1).  0(1) 

KsNn*NSlZE 
no  ion  1*1, K 
nin  r(  I  )sAT  I  )-Rd ) 

RETURN 

FNR 

%  FORTRAN  LSinil.OFCK 

r  HIM  IS 

SUBROUTINE  HU|  T  (  A ,  R .  C .  (.  I  Z .  N I Z .  H I Z .  MAX  A .  HA)(B .  HAXC .  Nr* .  NP  ) 
•i  IN,  NS  I  ON  Ad  ).  lid  ).  Cd  ) 

KA^KCtHAXA 
KBrNCtMAXH 
KCsNOHAXC 
no  10  (201).  iiliD.NP 

tun  CAM  OHULT  ( A . H . C . I  I Z . NI Z . M I Z. KA. KB. KC. NC ) 
no  10  700 

?iifl  no  Ann  i.si  .nz 

nil  6fl0  M=).HI7 
K=  (H-1)»Kn  »l 
1=  (M-1  )*Kr  ♦!. 
rt I  )  =  o. 

no  TO  <  1110,4  m ),Nr 
.tun  no  310  Nsi.Niz 

Js(N- , )*KA«t 


HTR4flV64  1 
HTRAO^OS  1 
NT  041)26  7  i 
MTR40269  I 

••TR4427.t  I 

PTR4427?  I 
HTR4ft773  I 


HTR40275  I 
HTR4fl?77  I 
HTR4n279  I 
HTR402BII  i 
MTR402B1  I 
HTR4A2A3  I 
MTR402Bb  I 
MTR4n2B6 
MTR402Bb  I 
HTR4n290  I 
HTR40292  I 
HTR40293  I 
HTR4n295  I 

MTR4U297  4 
HTR4n299  I 
MTR40301  i 


JB  :  (N-1  )*K 

•llrt  C(  I  )5C(  I  )♦*(  J)*B(  JR) 

no  TO  ^1:0 

41)0  IC  =  I*MAXC 

CdC)  *tt. 

DO  410  N=l.NIZ 
Js  (N-1  )*KA«L 


MTR4()J0J  p 

MTR40304 

MTR4n3no 

HtR4(!;i(i/ 

HTR4njn9 

HTR4IKn  1 


JH*  (N->  ) 

JC*  J»MAXA 
JRCs  JR^RAXO 

n( I )an( I )»A( J)«H< JR)-A( jn)»H( JRC) 

410  C(  IC)sn(  to  »A(  J)«R(  JRC)«A(  JC)*R{  JH) 

StlO.  CONTINDF 
AUn  CnHTINUH 

700  RFTORN 
f-ND 

%  FORTRAN  LSTntl.DFCK 

C  nHlH.TS 

SURROOTINE  DHOLT  < A , R . C . L 1 2 . N! 2 . HI Z , KA , Kb . Kf , NC ) 

DOUBLE  precision  A(1).  B(1),  C(1) 

MAXAsKA/? 

HAXPsKR/? 

?00  no  ADO  Lsl.LIZ 

no  *>00  Msi  .M|Z 
Ks  (M-1 )»KR  ♦  I 
Is  (H-1 )*KC 
C< I  )««. 

00  Ti)  <  ^il.(>«4iin).NC 

Ton  no  J)o  NsiiNiz 

Js(N-n«KA4l 
JH  S  (N-l  )*K 

Tin  .C(  I  )SC(  I  )fA(  J)#B(  JB) 

00  TO  soil 
40n  ICsI*KC/7 

C(IC)  =n. 

DO  410  Nsl.NIZ 
Js  (N-1  )*KA*L 
JRs  (N-1  )  ♦K 
JCs  J4HAXA 
JH'  *  .yR*MAXH 

n(  I  )=n(  I  )«A(  J)*H(  JR)-A(JO»H(  JBC) 

4  10  C-,  '  ')sn(  IC)  »A(  J)«B(  JRC)4A(  jn)*R(  Jb) 

son  CON  I  IN. 'I 

Alio  CONTINUh 
700  RETURN 
END 

%  FORTRAN  LSIOM. DECK 

n  POLHS 

SURROUTINf-  POI  M  ( PN ,  PN1  ,  ON .  ON  i  ,  E  I  HBDl  » LHRI)''* ,  NC,  I  P»  f  GO  ) 
niMINSION  PN<1),  PNKI),  ONO),  0N1(1)»  E?(1),  lHPI)i(1), 
DOUMI  I  PRICISION  PN.  PN1  ,  ON.  ONI.  LHBD1  .  LMbDl'.  A(?) 


HTR4(i  M/ 
PTR4fi:)l  3  T' 
HTRAOjU  Jo 
HTR4031 6 
PTR-i'l.l)  /  -- 
PTR4n3l-.  .I‘ 
HTR-i  -'3?  I 
PVR4  s')>  '  ^ 
PTRAf-^?-? 


HTRAO.i^/  “7 
HTR'IOJPV  £ 
HTR4  03.111 
HTR4  0.1.11 
MtRAHJT?  f: 
NTR4n.i.l4  *“ 
HTR40;<  16 
PTR40.1.1/ 
hTR4n:i.)9  U 
HTR4n34i 
MTR4«34J7f 
HTR40;14>>  Ij 
MTR4n.146  *' 
UTR40347  „ 
HTR4n34tt  If 
MTR40.1'>ti  ii 
hTR40.i6l 
HTRoOo'o" 
HTR4035^  I, 
MTR4n3‘l6 
HIR4(i.l6;  .. 
MTRAO.I'jn  I' 
PTR-inj^l) ' ' 
HIR40;iA) 
HTR40.16,1  Y- 
MIR-MMAb  i;. 
HTR4n.1A/ 
HTR4II369  T' 

I 

!  • 

f1TR4ll.i7^  . 

LHBD'.'d  )HTR4fl3;4  f 
MTR-iii.lZ'j  ■'  • 


IR  =  n 

DO  10  (ion.  Ill  <1.1  HO),  I  no 

lofi  DO  in  (1  )n,'>iiii  ),NC 

110  If  (  nAHS(  (PN-PNl  )/(flS0RT(  nARS(ON)  ))  )  -FX  )  1i»  .lim.ll? 

11?  IF  (  DARS(  PN-PNl  )  -  )  1  iri. ,  1 PO ,  1 /O 

1?o  IF  t  IIARS(  (0N/()N1  )-1  .  )  -  F2  )  1  3  ■  ,  1  30 , 1 /« 

1111  I  HH|i?  =  (PN«PN  -'l.*ON) 

I.MPli?  =  nSORT  (DAHSd  MRD?)  ) 

I.MPni  =(-PN  ♦  I  MHO?)  /?. 

I.MRD?  =(-PN  -  LMBII?)  /?. 


MTR4n.176 
NTR4  0.17  /  j  j 
MTR4n379  I  I 
H1R4I13HI 
HTR403H2  : 
HTR4n3P4  I 
NTR4 11306  •  • 
HTR40387 
HTR4n3R3  | 
HTR4fi;iB9  1  . 


IF  (  ARSILhRDi)  -  ABS{LMHD2)  )  140,1^0,16(1 

140  4=  I MBD1 

LMRni =l MBO? 

I.MRI)?sA 
RO  TO  160 

isn  IF  (  LMflDl  )  140,160,160 

160  IRrNC 
I/O  RETURN 

IHO  RO  TO  (l3n,J'9iO,NC 

?.tO  A  =  nSORT(  ON(1)*QN(.1)  ♦  0N(R)»0N(7)  > 

IF  <  DSORTI  (  (PN-PN1 )#(PN-PN1 )  ♦  ( PN( 2 )-PN1 ( ? ) )*( Pn( 2  1-PNl < 2 ) > 
1  /  A  )  ,-F2)  210,213,170 


140,160, 160 


/  A  )  ,-F2)  210,213,170 

210  A  =  ON1»ON1  ♦  ONI  (.2)*IJN1<?) 

IF  (  (  <  (0N»0N1-1}N(2)»0N1  (2))/A  )•{  ( bN*ONl -ON(  2  )«'’N1  ( 2  )  )/ A  )  ♦ 

1  (  (UN»ON1 (2)40N1*0N(2))/A  )•(  ( ON»ONl ( ? ) ♦QNl »ON( V ) )/ A  ) 

2  -1.  >  -  F?)  220,220,170 

220  A(1)  =  PN#PN  -  PN(2)«PN(?)  -4.*0M 

A(?)  =  2.»(PN#PN(2)-?.»qN(2) ) 

IF  (  A<1  )  )  230,2»>n,?'>6 

210  Ad)  =  nSORTI  (-Ad  )  +  DSORT(  A  ( 1  )*A<  I )  tA  ( 2  )  ^A  ( 2  )  )  )  /?.  ) 

A(2)  =  A(2)  /  (2.#A(1  )  ) 

212  LMBIll  (1  )  =  (-PN*A(2)  )/?. 

LMROK?)  =  (-PN(2)*A(l))/2. 

LMBn2(1 )  s  (-PN-A(2) )/2, 
l.HRD?(2)  a  (-PN(2)-A(1  ))/2. 

Ih  (  DSORT(  lMR01*LMBni*LMBDl (2)*LHBD1(2)  )  -  DSORT(  LMHn2#t«8n2 
1  LMBD2(2)»LHBD2(?))  )  240,160,160 

240  Ar|MBni{2) 

LMBOl (2)rLHBn2(2) 

I.MRI)2(2)a  A 
00  TO  140 

2>>n  Ad)  a  nSORT(  (A(l>40S0RT(  A  (1  )#  A  (1  )♦  A<  2  )#A  ( 2 )  ))  /2  .  ) 

A(2  )  a  A(2)  /  (2.*Ad  ) ) 
l.HRIil  aAd  ) 

Ad  )rA(2) 

A(1>  )aLHBJ)) 

00  TO  2.T2 
FNO 

FORTRAN  LSTOll.OFCK 


MTR4n39l 

MTR40392 

MTR4n393 

NTR4n304 

MtR40395 

HTR40396 

MTR40398 

HTR4040U 

HTR40401 

HTR40403 

)HTR404n6 

NTR404«6 

HTI<404n7 

NTk‘f404nv 

NTR4n41ft 

MTR40411 

HTR40412 

HTR40413 

flTR4041!> 

HtR40417 

MTR40418 

,RTR40420 

MTR4(I4?1 

HTR40422 

NTR40423 

♦MTR4n42t> 

HTR40426 

HTR40427 

HTR40428 

MTR40429 

HTR40430 

MTR40432 

HTR4n433 

MTR4n434 

HTR40436 

NTR40436 

HTR4fl437 

NTR40439 


BURROUTINE  PO  (FL,  HN,  HNl ,  HN2,  HN3,  P,  0,  NC,  HAyR  ) 
niMFNSION  Fid),  HNd>,  HNl  ( 1  ) ,  HN2(1),  HN3(1),  P(li,  0(1), 


NTR40441 

HTR40442 


I  A  ( .1 ) .  R  ( ,1 ) 

nmtMI  E  PRtClSION  FI,  MN,  HNl,  HN2,  HN3,  P,  0,  A,  R 
MA)(H  a  MAXR/2 
00  10  (2110,1110).  NO 

100  Ad)  a  rL(l)#HN  -  FL(2)#HN(HAXR41) 

A(2)  s  rL(3)*MNl  -  f L(4)*HN1(MAXR41 ) 

Ad)  =  Fl-{‘>)»MNV»  -  FL(  6)#HN2(MAXR4l) 

Rd  )  =  d  d  )*HN(MAXR4|  )♦  FL(2)#HN 
0(2)  =  FL(.1)«HNl(HAXR4l  )4FL(4)«HNI 
R(  i>  a  FI  (6)«HN2(HAXR*1)4FL(  6)«MN2 

Pd)  a  A(3)*mn'i  -  B(3)»HN1  (HAXR41  )-A(l)*HN3  ♦  B  ( 1  )#H‘i3  ( KAXR+l ) 
P(2)  =  A(3)<.HN1  (HAXR41  )4B(3)«HN1  -  A ( 1 ) »HN3 ('HAXR*1 ) -P ( 1  )#HN3 
Ud  )  =  A(l)*HN'’  -  Rd  )«HN2(MAXR41)-A(2)*HN1  ♦  B  ( 2  )#Hm1  ( MA  XR*  1  ) 
0(2)  =  A(l)*HN2(MAXRM  )4Bd  )»HN2  -  A  (?)  •HNl  ( HAXR4 1 )-'!(?)  #hNl 
Ad)  =  A(2)«HN«  -  B(2)#HN.1(MAXP41  )-A(3)*HN2  ♦  B  ( 3  )*Hf2  (  MA  XR4 1  ) 
«(1)  a  A(2)«HNMMAXR41  )4B(2)#hN3  -  A  (  3  )  ♦HN2  ( MAXR*  1 ) -M  3  )  tHN? 

A(2 )  =  A( I  )*P( 1  )  ♦  h( 1  )«P(2) 

R(2)  =  A(1)«P('>)  -  H(l)«Pd) 

A(^)  a  Ad)«0d)  ♦  H(1)«0(2) 


MTR40443 

HTR4fl44‘> 

HTR4n447 

HTR4fl448 

HTR404*»0 

HTR40451 

MTR4045? 

HTR40464 

HTR4046‘> 

MTR40456 

HTR40458 

HTR4n459 

HTR40461 

HTR4n462 

MTR40464 

HTR4n46t> 

NTR4U46/ 

NTR40468 

MTR4n469 


88 


B(.l)  s  A(1  )»0(?)  -  8(1  )*Q<1  ) 

A(1 )  =  A<l)*A(1  )  ♦  8(1  ,*R(1  ) 

A(?)  =  A(ii)/A(1  ) 

8(?)  a  B{i?)/A(1  ) 
hlS)  s  A{i)/A(1 ) 

R(.l)  a  B(3)/A(1  ) 

P(1)  a  FL(3)*A(?>  -  FL(A)*B(?) 

P(?)  a  FL(3)»B(2)  ♦  FL(4)»A('>) 

A(?)  a  FL(3)*FL(*i)  -  Fl(4)*FL(  6) 

R(P)  a  Fl(4)*Fl(%)  ♦  FL<3)»FL(  6) 

0(1  )  a  A(2)«A(.n-H(2)«R(3) 

0(?)  a  A(2)»R(3)4A(.1)*R(?) 

'4AXP  a  P«HAXR 
RCTURM 

Ptift  A(1)  a  (FKJ'iaHNI  *883  -  FL  ( 3  )*HNP*HN7  ) 

P(1)  a  (FL(3)«HN‘>(1  )*HN1(1)  -  FL ( 1  )*HN(  J  ) *883 ( 1  ) )  /  a(1) 
0(1)  a  (  FL(1  )*H8(1)*HN2(l)  -  FL(2)»H81( I )*HNl ( I )  )  /  A(l) 
P(1 >  a  P(1 )*IL(’) 

0(1  )  a  0(1 )*Fl(2)*FL(3) 

HAXR  a  ?«HAXR 

RETnR8 

P8n 

%  F0RTRA8  LSTOU.DFCK 


HTR40A7fii 

MTR4047 

8184047^. 

MTR40474 

8184047  - 

818404?^^ 

81840476 

8TR40479p 

8TR4046  I 

8TR404Rf 

8TR404B4 

8TR404ft'T 

8184046  4 

H1R404B6 

HlR4049r^ 

81R4049  J 

81R4049,f 

81840494 

8184049'! 

8184049# 

81R40499 

818405)01 


.4 

c 

c 

n 

c 

c 

c 

c 

c 

c 

c 

n 

r 

c 

c 

c 

n 

r 


niQSBS 

S(JBR0Uli8E  (iLOSF 


COMPUIES  2  CLOSE  ROOIS. 


IJ  a  8AlRfX.  niHENStOHED  (HAXR.2«NC*HAXR) 

H  a  S1AR1I80  OUESS)  0 1 HE8S10NEO( (HAXR* ?»NC*4 )«?*mC*8 ) 
NSI7E  a  SI7E  OF  HAIRIX 

HAXR  a  Ul8E8SI08En  8U8RER  OF  ROWS  OF  U  a8D  H 
8AX1RY  a  8AX(8UH  8U8BER  OF  DOUBLE  PRECISION  I1ERa1I08S. 
EPSl  a  SI8BLE  ROOl  C08VER0ENCE  CRIIERIA 
EPS?  a  double  ROOl  C08VERGENCE  CRIIERIA 
Ra  AllKENS  C08VER0E8CE  CRIIERIA 
IRR  a  error  RFl(iR8  INOICAIER.  «1 

a2 

a3»  ROIH  OVERFLOW  AND 


overflow 

DIVIDE  CHECK 
ROIH  OVERFLOW 


1 


DIVIDE 
CHECK. 

IlERSa  NlIHRER  OF  IIERAIIONS  PERFORHEO.  •  FOR  DOU'LE  ROOl 

♦  FOR  SiNnLE  ROOl 

2,  COHPLFX 

(U.H>NSIZE.HAXR#R#EPS1,EPS2,8C# IRK. 8AX1RY* I lERS 
NAITKN. INDEX1 . INDEX? • VALUE* 8SI ZE ) 

VALUFd  ) 


NC  a  1 

siibrouiine 


•  REAL 
Cl  OSES 


DIM) NS  I  ON 
CAM  OVFRFL 
CAI I  DVCHK  ’ 
IRRrO 

NXa;^»NSI7E 
N2Cs2#NC 
CAI I  CHANOE 
CAI I  CHANGE 
l6aN?u«NSIZi 
Ilal 

I  ?  a  1 1  ♦  I  A 


11(1 ) 
( 


H(  I  ) 


inVFLW  ) 
inVDCl  ) 


(«I,HSI7F.NC*8SIZE.NA 

(H.NSI7E.NC*NSIZe.l) 


4  # 


1  ) 


81R40'i0,^ 
HlR40t>0  j 
81R405>05' 
81R40»>06 
818405(0  J 
8lR405>oX 
81R40509 
81R405>1  V 
8184051 j 
81R405lf 
81R40SU 
8184051  f 
8184051  rJ' 
81840516 
8184051  7* 
81840511^ 
81840519 
,81R4fl5?V 
818405? 
81840524* 
81840526 
8184052”: 
MTR4052., 
81840529 
8184(153' 
8184053, 
81R4053.r 
8784053^ 
8184053(1 


I3al2*16 

I4al3»16 

I5a|4»I6 

K1r|1 

K4a  12 

K.1aI3 

K2aI4 


81840538 

81840  53 '-f 

8184fi54(j, 

81840541 

81840542- 

81840547! 

81840544*’ 


89 


ITPRS  s  0 
Ksn 

IGO  8  1 
fTf-RSsITERS^l 

IF  (  !TERS-M4XTRY  )  14(»,140»130 

IF  {  K-ft  )  13?«140«13? 

fTFRS  s  ITER9-1 
GO  70  RflQ 
KsK«1 
IsKI 
K1  sK4 
K4SK.1 

•(2sl  ! 

OAll  HULT  (II  .H(K?).H(Kl)«NSlZF*NStZe,l,HAXR»N<;| 

1  NC.?) 

INOFXsfl 

iKsr*5*<4-K)«N;*r. 

CAM  norm  (H(Kl)«H(Kl)«NSIZE*H(tK)*|N0EX*NSlZE«NC*2t 
CAII  OVERFL  (  lOVFLW  ) 


.H(K?).H(Kl)«NSlZF/NStZe,l  ,HAXR»N<;iZE.NSlZE. 


Ufl  CAI  L 
GOTO 


?nn»9ft«»44ii 


?20.?2n«30fl 


CAII  OVERFL  (  lOVFLW  ) 

GOTO  (190.1411)  .lOVFLW 
ISO  |RRsIRR4l 

Ufl  CAI L  DVCHK  (  IDVUCT  ) 

GOTO  (I7«,l«(i)  .lOVOCT 
1/0  |RRstRR«7 

190  IF  (  IRR  )  ?nn»7ft«.440 

TE9T  FOR  CONVERGENCE  TO  A  SINGLE  ROOT 
7M«  GO  720  I«1  ,16,? 

jy  X  K?*I-1 
J.i  s  K1*l-1 

IF  (  ABS(H(  J7)-H(  J3)  )-eP5J1  )  ?20.?2».3o(l 

770  CONTINUE 
GOTO  790 

30  0  GOTO  (j7u.1’tl.l7H,.1?()),K 
.170  JI  =  1G4(  |G04o)*M?C 
J2=I94( IOO*9)tN2C 
.i  =  7#INI»FX-7 
J.)sKl4j 
J‘>aK74j 
J/rKl4j 
JVaK4*J 

CGNPilTE  P  N  AND  0  N. 

CAII  PO  (H(l<)).  H(J3).  H(J9).  H(J7}.  H(J9).  H(J1).  ii(J7).  NC. 
GOTO  (.190.340  ) .  IGO 
MO  .Ms  I9«4»N7C 
.)7r  JJ4N7C 
.Ms  .I74N7C 
.14=  J34N7C 
.19=  I94N7C 
KSs  J94N7C 

TEST  FOR  DOUBLE  ROOT  CONVERGENCE  AND  IF  SO*  COMPUTE  LAMBnA  1  AND  2. 
CAI  I  POLM  (H(.I7).H(J1  ).H(  J4).H(J3).ePS8*H(J9).H(J6).»>C*  I R*  1001 
GOTO  (,144.344.4011).  IGO 


GOTO  (,144.344.4011). 
IF  (  IR  ) 

11  (  ITFRS-MAXTRT  ) 

no  340  Is|.N'»C 
H(  .11  )sH(  J7  ) 

M(  J.1)=H(  J4  ) 

.ns.Jl4l 

.M  S.)?4l 

J.1S  J3  4l 


344.  14^.400 

.147.800.801' 


MTRGO'iAS 

MTR4G54? 

HTR40948 

HtR40549 

lltR40550 

|«fR4055l 

MTR40552 

HTR40983 

HtR4p‘>5'> 

HtR40996 

HTR40997 

MTR40998 

MTR40959 

HTR4Ub80 

HTR40962 

HTR40963 

HTR409A4 

MIR40945 

MTR40946 

MTR40948 

MTR40969 

MTR40970 

MTR40971 

MTR40972 

MTR40973 

MTR40974 

HTR40976 

MTR40978 

MTR40979 

MTR40980 

MTR40981 

MTR40982 

MTt409»3 

MTPAt^A'i 

MTR4098? 

MTR4(t988 

MTR40b90 

MTR40*)91 

MTR40‘>9? 

MTR40993 

M1R4fl994 

MTR40996 

NX)MTR40997 

MTR40600 

HTR406fll 

MTR40602 

HTR40603 

MtR406R4 

MTR4060b 

MTR40606 

MTR4n6n8 

MTR40610 

MTR40611 

MTR40612 

MTR40613 

MTR4n614 

MTR4061b 

MTR40616 

MTR40617 

MTR40616 

MTR40919 
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J4  =  J4  +  1 

.T>n  K=3 
IG0=? 

00  354  Jsl.< 

l.sl5*(3-J)»N^n 
Llsl5  +  (4-J>*N,>C 
on  354  1=1, NC 

LsL*( Irl )#NC 
LlsLl ♦{ 1-1 )»NC 
Hd.l 

^94  H<U  ♦l)sH(L*1  ) 

CAll  AITKNS  (H(iKn,  H(K?).  H(K3),  H(K4),  R,  '  ,,6,  tsSIZF. 
1  IR  ) 

<  IR  )  3<)f),  ij»n,  ^611 

36«  Ir|t1 
•(1=r4 
K4rK3 
K  J  =  K? 

K2sl 

NAITKN  s  NA1TKN«1 
GOTO  inn 

400  CAI.I  POH  (H<1'>1,  H<J5),  H(J6),  H(K1),  H(KZ),  H(K3),  u<K4), 
1  NO 

GOTO  (404,40J»1,NC 
40?  VALUE(9)*H< J5*o) 

VAU)E<4)sH(  JA  +  O) 

404  IHOFXrO 

VAl UF<1 )sH( J5) 

VAlUE(NC*l)sH(jn) 

GAll.  NORM  (H(K3),H(K3).NSIZE.H(J1). INDEX. NSlZB,NC»y) 
INOFXl  s  INDFX 
ITFRS  s  -ITFRS 
isK442*|NDEX-;^ 

H(J1 )sH( I  ) 

H(J1+1)*H(I*1 ) 

GOTO  (4?n,4tn),NC 

410  I  s  ffNX 

H(J1*?)sH< 1 ) 

H(.I1  ♦3)=H(  l4l  ) 

420  INOFXs-d'OFX 

hail  Nn.  H  (H<K4).H{K4),NSIZ6.H(J1).1N0FX,NSIZE,NC,'  ) 

I  ’  440.500,450 

440  J.j=13 
44?  .14  =  14 
470  K=1 

460  GO  46?  J=1 .  |6 

J1 =J3+ J-1 
.»?- J44  J-1 
46?  H(JnsH(j2) 

GOTO  (5?i»,5Hi.51'i.f,nii),K 
4/0  .I.5M1 

GOTO  44? 

400  (F  (  K.1-  13  )  490,490,4711 

490  K=3 

13=  M 
.14=  1.3 
GOTO  460 
500  Kso 
.M=|3 
.14  =  11 
GOTO  46(1 


HTR406PU 
MTR40622  ;* 
MVR40623 
HTR4n624 
MTR4  06?5  - 
HTR4.)6?6  j 
HTR4n6?7  '■ 

HTR406?a 

HTR4n6?9 

HTR406.3a 

HTR406.31 

NC.  2,  HTR406Ti3~' 

HTR4('f;  »  4  . 

PTR4063/  - 

HT»4fi6^H  ' 

HTR4»h.ly 

HTR4064U 

HTR4n64i 

HTR40642  « 

HTR4064J 

NSIZfc.  MTR4n645  - 
MTR40646  * 
MTR40648 
HTR40649 
MTR40650  T 
MrR4n651 
HTR4065? 
HTR40653  ■“* 

MTR40654 

NTR40656  " 

HTR40657  .. 

HTR40658 

HTR40659  -* 

HTR4n66(l 

HTR40661  ■' 

MTR40662 

rTR4!)663 

HTR4n664  • 

MTR4n666 

MTR4n66/ 

MTR4066fl 

MTR40669 

MTR40670  - 

HTR40671 

MTR4n67i; 

MTR4I1673 

MTR4n674 

HTR4n677 

hTR40676 

MTR4067/  ■ 

MTR4067H 

MTR40679 

MTR4n6H0  , 

MTR4n6fll 

HTR406P2  •' 

HTR406«v3 

HTR4n6fl4  '• 

HTR  •'6115 

HTR4il»..  • 

MTR4nAR/ 
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S  1  ft 

.11=11 

MTR406RM 

SI? 

14=  1? 

MTR40AR9 

GOTO  4‘)ll 

MTR4n69U 

S  '  4 

H  =  I.T 

MTR4IIA9i 

GOTO  *>1? 

MTR4n692 

•i  xr> 

CaU  mads  (H( 11 ).H( |3)«H( I2).NSIZE*NC) 

MrR40O9J 

.JsNC»NSIZe*l 

INHf-Xsn 

HAll  NORM  (H( J).H(J).NSlZE*H(l&)*tNDeX.NSt7E,NC»l) 
J=J*lNnFXl-l 
Ji  =JMNn-i)«Nsize 
M  (  J  >  n . 

M(.ji  )sn. 
mn^x?  =  iNDFx 
CAM.  OVERFL  (  lOVFI.H  ) 
aOTO  (ftl  )  .lOVFLW 
<Sin  IRR=IRR*l 

(s/n  CAM  ovcHK  {  invncT  ) 

ROTO  (63U.64IO  .IDVOCT 
IRRsIRR*? 

<S4n  CAM.  CHANGF  <  II.  HSI  ?F .  NC*HS  1  Zb .  HAXR*? ) 

700  RETURN 
7bn  tKsl5*(4-K)«NVn 
VALIlFd  )  s  M(  IK) 
lKslK*?*(NC-1 ) 

VAI  UE(Ni:)sH(  IK) 

INIir  XIsINDEX 

Ks4 

JJell 

HsK1 

IF  (  Kl-!1  )  46n,ftU(i,-l6n 

inn  IGOr.l 

INDEXlslNIlFX 

Jsib 

.11  =  IS*4*f^'.- 

nilA  CALI  LEGIS  (  H(J).H(J1)>FPS1.EPS?.NC.1R  ) 

IF  (  IP  >  34n.Mn.7^>n 

FNH 

FORTRAN  LSTr>:.DFCK 
MITFRS 

A  IF  FTORFD  IN  CORE  AT  A.  (HAXR  X  NC«NP»NSIZE  > 

NTAPIIT  IS  A  UTHITY  TAPF,  FOR  CHECK  VECTORS  IF  OESIREn. 
hPSP  =  FPSItON  ONF  =  SINGIE  PRECISION  CONVERGENCE  TEST  NUMBER 
f  PDP  =  FPSiioN  run  =  nniiHLE  precision 

NC  =  1 .  IF  RFAl  NP  S  1.  IF  SINGLE  PRECIS 

7,  IF  COMPIFX  =  IF  DOUBLE 

NGIIFSS  =  n.  IF  FIRST  GUESS  IS  TO  RE  A  COLUMN  OF  ONES. 

MOlinilT  s  NO..  OF  MODES  TO  RE  COMPUTER. 

NAKSH  -  NO.  TIMES  AIFKNS  ACCELERATION  NAS  USED  IN  SINGL*  PREC 


H(J).H(J1)>EPS1.EPSP.NC.IR  ) 
340. ^40.7^0 

DECK 


=,  <?.  IF 
COLUMN  OF 


SINGLE  PRECISION 

DOUBLE 

ONES. 


NAKSH 
NAKOR 
MAXSR 
MAXDR 
IRH  = 


NSi  ;r 
RSP  = 
HDP  s 
MAXR 


-  NO.  TIMES  AIIKNS  ACCELERATION  NAS  USED  IN  SINGL*  PRECISION. 
-..  .•  ..  .•  ..  ....  DOURL^'  . .  . 

'  MAXIMUM  ITERATIONS  ALLOUED  IN  SINGLE  PRECISION. 

^  ..  ..  ..  ..  Double  . .  . 

I RROR  HE1URN  ^  t.  FOR  OVERFLOW 

FOR  DIVIDE  CHECK 

.W  FOR  ROTH  overflow  AND  DIVIDE  C.ECK 
t  NO.  OF  ROUS  AND  COLUMNS  OF  A 

H.  AITKENS  ACCFll RATION  CONVERGENCE  CONTROL  FOR  S'NGLF  PRECIS. 

R  AIIKENS  ACCFI FHATION  CONVERGENCE  CONTROL  FOR  lOUbLE  PRECIS 
:  DIMENSIONED  NIIMHER  OF  ROWS  OF  A  AND  GUESS 


SURHOIITINE  MITFRSI  A.OUESS.NGUESS.NSiZE.MODOUT.MAXR.Nr.EPSP.EPDP. 


MTRAOdR*) 

HTR40697 

MTR4l»69b 

MTR4n7l»ij 

MTR417PI 

MTR48;K3 

MTR4ll7fl4 

MTR4n7ni> 

MTR407n7 

MTR4fl7IIB 

MTR4n709 

MTR407n 

MTR4n71? 

MTR4(»713 

MTR4071‘> 

MTR4n71 / 

HTR40/19 

MTR4U7?» 

MTR4fl7?l 

MTR4n7?2 

MTR407P3 

t*TR4n724 

MtR40  7?!> 

MTR4n7?ft 

MTR4n7?/ 

MTR407?y 

MTR4n730 

MTR40731 

MIR4II732 

MTR4073J 

MTR40734 

MTR4n;.36 


NTR40738 
MTR407.19 
MTR4n74(i 
MTR40741 
HTR4n742 
MTR4074 j 
HTR40744 
MTR4n74‘> 
MTR4n746 
MTR40747 
MTR4n748 
NTR4n749 
MTR4II78II 
MTR407^l 
MTR407‘>if 
MTR40783 
MTR4n7^4 
.MTR4«7'»!> 
MTR4a796 

MTF'4«758 
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1 

NAKSR 

.NAKPR.MAXSR.MAXUR. 

RSP,ROP, 

IR> 

,H 

,  VECTOR. 

MTR4n7»,V 

?. 

VALWF 

•NITER, NTAPUT, NTAPOT) 

MTR4ii7Ai(  -p 

OIMINSION  A(1).  GlIFSSd).  H(1>.  NITERd) 

,  NAKSRd  )  . 

NAKDRd  )  . 

MTR4I'  76  J  j 

1 

ATI TLF(N)  . 

VFCTORCi ),  VALUE (1 ) 

MTR4n7A4 

OATA  ( ATITLFI I  ) .  l  =  1 .N) 

/  UiH  OVERFLOW. 

.  1  OH 

• 

MTR4n/At)  * 

1 

1  IIH  PI  VIDE  .  1  DM  CHECK 

.  lUH 

IN  , 

MTR40/A/  [ 

2 

IliH  CLOSF 

.  l»H  SHEEP 

,  lOHSUBRO" 

TINE  , 

MTR4n7Atf*' 

3 

1  IIH 

/ 

MTR407AN 

in 

FORMAT  (IHfl  •iX.  AH  MOOE  l^X.  IIH  EIGENVALUE 

loX. 

MTR4077]  '{ 

1 

24H?TFRAriONS 

s.p.  n.p.  AX, 

?0H  AITKFNS 

S 

.P.  D.P. 

MTR4II77?  J 

? 

/  U  ) 

MTR4n77.l 

1? 

format  dH  llli.  ?F19.ft#  119?,  11/, 

1119, 

1  1  ' 

) 

MTR4n7  74  -1 

\  4 

FORMAT  (IH  1111.  IF2R.B#  iIU.  IW, 

1119, 

1  1  » 

) 

MTR4n  /7‘>  J 

! 

t 

FORMAT  dH'i  1  U.  43H 

MOOES  ARE  CORRECT, 

COMpUTATIO* 

TERMINATED. 

)MTR4{I776 

FORMAT  (IbHliHOniFlFII 

MOOF  1|A,  ??M  IS 

CORRECT, 

Ti>UF 

MODE 

MTR4n7d 

1 

?liH  CANNOT  OF 

COMPUTED.  ) 

MTR40  7  7f<  1 

<'(1 

FORMAT  (1HO  ///  1 HH 

4AX,  t4H  EIGENVECTORS 

/  ' 

/ 

) 

MTR4fl7  /  V  J 

FORMAT  dH1  mX.  .'AH 

CHECK  EIGENVALUES  AND  EIGENVE" 

TORS  n 

MTR4I!  7ft.| 

1 

(AHA. 8)  ) 

MTR4F7* ;  1 

r 

FORMAT  (5H  MODE  114 

.  ?0H  HAS  NOT  converged  IN 

MA/IHUM 

MTR4(i7f.i.  J 

1 

i 

l?H  ITERATIONS 

) 

HTR4II783 

FORMAT  (37H1*»«**  ERROR  IN  ITERATION  SUBROUTINE 

PAID  > 

MTR4n784  . 

CAIL  OVERFL  (  lOVFLH  ) 

MTR4n7Bt> 

1 

CAll  OVCHK  (  IDVnCT  ) 

MTR407ft/'^ 

1 

HSI7F  =  NSI2F 

MTR4078H 

ISIZE  s  NCtMAXR 

MTR4II7B9  j 

r 

IRRsO 

HTR40790  1 

1 

riNO  A  AND  STORF  ON  TAPF  IF 

NECESSARY. 

MTR4(I79? 

IF  <  NTAPUr  ) 

Id?, in?, 100 

MTR4n79vi  • 

r 

1 1I  n 

REHIND  NTAPIIT 

MTR40/94 

MsNSIZE«ISIZE 

MTPnnzRb  * 

0 

MRIfE  (  NTAPUT  ) 

(A(  I  )dsl  ,H) 

MTR4n79A 

REWINO  NTAPUT 

MTR4079/ ' 

r 

1 1I? 

IF  <  FPSP  ) 

1l|A,1  114, 1UA 

MTR4()799  .. 

i 

1114 

FPSP  s  .^F-liH 

MTR408(ilJ 

1  IIA 

IF  (  FPDP  ) 

1in, 1118,110 

MTR408n? 

>1^ 

1  (i« 

EPFP  =  ,^E-i'n 

MTR408nj, 

i 

rl 

1 1  n 

.IS17E  =  NSIZE*(NC-1  ) 

MTR4n8nt,  ' 

OFFINF  FIRST  GUFSS.  IF  NOT  GIVEN. 

MTR4n«n  7 

IF  (  NGUESS  ) 

111,114,111 

MTR40HI)9 

Ill 

IF  (  MA^  -NGI7F  ) 

1  1 A,1 JA,1  1? 

MTR4n81 0 

*P 

,  1? 

00  11.1  Isi.nSIZE 

MTR4081 1 

Jr  NSl/1  ♦  1 

MTR4n81 2 

Jl r  MAXH* 1 

MTR4«H1  .1 

1  1  3 

miFSS(  J)=utirss(.ii  > 

MTR4IIH1  4 

GOTO  1  1  A 

MTR40HH 

- 

1  14 

no  IIA  |sl.NSI7l- 

MTR4n81 7 

JrJSIZFM 

MTR4()H1  « 

nuFSS(J)  r.  II. 

MTR4r)H1  9 

1  1»» 

niiFssi  ?  )  =  1 . 

MTR40rt?u 

OFFINF  PROGRAM  I.ONSTANTS. 

MTR4(IH?? 

1  )  A 

MOOF  =  11 

MTR408P4 

n  d 

MTR4n8?t, 

I?  =  11 ♦ IS! 7F 

MTR40H?6 

1.1  r  1?*  ISIZE 

MTR408?/ 

S>  T 

FOR  MOfif  COHN  I 

MTR4n8?9 

1  1fi 

Moni  =8001  ♦  1 

MTR4n8Jl 

»fsi/f  r  nchsi/e 

MTR4np,l? 

II  (  MoiiF-  Honoin  ) 

1 4  n ,  1 4 11 , 

,  ^  il  II 

MTRnilH.I.I 

1  -1 II 

NSRAK  =  (l 

MTR40M^V 

«  • 
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NORAKsO 
ITFPSRsn 
ITFPnRsn 
K1  =  I1 

Ki=I? 

no  ISO  Ksl.NC 
J1s(K-t )*NSI7E 
Js  (K-1  )*M«;iZE*K1-1 
no  ISO  Fsl.HSIZF 
JsJ*1 
JlsJI ♦! 

Ion  H(.l)xGui-SS(  J1  ) 

1  >9  N*r;(l 

l-'.fi  RArrNAK*1 

IS’  ?TFRSR=lTbPSR*1 

IF  (  ITFRSR-HAXSR  )  1/0,1/O»?*>« 

1/n  isKi 
K1  =k.^ 

K3sK? 

K?s| 

C  SFT  ....  NOW  MAKF  ONF  ITERATION. 

r 

''At  L  HOLT  (A.M(K9).H(K1  ) .  HS 1 ZE  •  NSI ZE » 1 .  HAXR. HSI ZE*  HSl  ZB » NG .  1 ) 
INnr  xso 

IAS  M*KSlZE*Nr»(  WNAK) 

CAil  NORH  (H(Kl).H(Kt ).H8iZE.H(IK). INDEX, HSIZE»NC»1) 

CALI  OVFRFL  (  lOVFLU  ) 

CAM  nVCHK  (  IDVilCT  ) 

OOTO  .lOVFLW 

IHO  IRRsIRR+1 

1B9  ROTO  (104. J«6)  , IRVOCT 
184  IRRsIRRt? 

1H/S  IF  (  IRR  )  Ipn.IPn.OOfl 

I'iO  GOTO  (l6il,?<in,?n>i),NAK 

9110  no  ?in  lal.KSiZE 
.ll=K1*I-l 
.l9sK?+I-1 

II  (  ABS(H{.ii)-H(j?))-EPSP  )  >ifn 

’in  rnNTINUE 
goto  400 

990  GOTO  (IftO,  J  ),NAK 

910  CAIl  AITKNS  (  H(KI).  H(K?).  H(K3).  H(K3).  RSP,  NSIZE.  HSIZE#  NO, 
I  1 .  IR  ) 

II  (  IR  )  940,9,19,940 

040  |.s9«Nf: 

.l=IK  ♦!. 

Ilf)  9  14  Isl  .1  ,Nr 
Jl * J*Nr-| 

.»9s.)-I*1 

.1  ^^.)1  -NC 
M  r  .)9-NC 
H(JI  )sH(  j.n 
9  14  H(.J9)=H(J4> 

GOTO  1f>9 
9  4  0  I  - K  I 
K1  -  A  1 
K  1  K  9 
I 

nsrak=n<;rak«' 

GOTO  1S9 


'll 

RTR4(ie?6  } 
flTR406.i;  [ 
HTR40«3n  1; 
NTR40a40 
HTR4n841  I 
HTR40D42  j 
NTR4a844  t 
HTR4084'>  5 
••TR40R46  3 
mtP4!)i»4> 
nTR4%fe«!y  I 
l•|R4AA>4?  ’j 
NlR4»JhSH  I 
NTR40HS2  ;] 
HTR40bSJ 
HTR4nflSb  I 
NTR408S/ 
MTR4n85V  ;i 
HTR40860  J 
MTR40b61 
MTR4nAfS2  I 
HTR408A4  I 
NTR4086t>  J 
HrR4nR66  ^ 

NTR408AIJ  i 
HTR408A9  : 
HTR40870 
HTR40872  I 
HTR40873  ! 
htr4087*>  *; 
HTR40876  i 
NTR40H77  1 
MTR40878  : 
HTR40HR0  i 
NTR4n8fl?  j 
f<TR408l14  1 
NTR408ftS  ■ 
HTR40886  j 
NTR408fl/  ■ 
HTR4n«H8  = 
MTR40809  j 
HTR40891  j 
HTR40893  i 
MTR40H94  ■ 
MTR4U89‘> 
HTR40896 
MTR4089/ 
HTR4n89H  1 
MTR40899  : 
NTR4n90U  !i 
NTR4n9ni 
HTR40902 
?iTR40903 
HTR409fi4 
HTR409n^ 
NTR40907  ■ 
HTR40908 
HTR409rt9 
HTR40910  i 
HTR4n9n  ; 
HTR40912  ■ 
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'»'>n  HAI  I  nvf-Kri  (  invU.M  ) 
r.Ai  i  nvCMK  (  invncT  ) 
nOTfi  .lOVFLW 

JRp5|RR»1 

9*’  90Tn  .invncT 

9f.4  ?RRsIRR4? 

IF  (  IRR  )  ?7fl.‘>7rt,/.un 

9/n  IF  (Kl-Il  ) 

•>!•?  no  71^  I  =  J  .KSIZfc 
J=Kl+l-l 
J1=I1 ♦!-! 

?/4  H( J1 >  s  H( J) 

980  Js(H0DF-1 )*Nr*1 

ITFRSR=|TERSR-1 

CALI  CLOSES! A.H( n  ) . HSl ZE • MAXR. ROR. EPSP . EPOP* NC, IRR 
1  NRRAK.  INDEX1  *lNnEX2.VALtlE(J).  NSI7E 

IF  (  IRR  )  7HP»'»«9,A10 


IF  (  IRR  ) 

P89  IF  (  ITERDR  J 


p«*:»98R,9bft 


IF  (  KSIZE-TSIZF  )  P«*:.98R,9bft 

9H4  II  =  I1+2*KSI7F 
12  =  IP^KSIZF 
no  284  Isl.KSIZE 
Jl=  Jl-l 
J2=  J?-l 

984  H( J2)aH( Jl  )  ^  ^ 

2»R  INREX  a  INOEXI 
iNniX  s  INDEX 
29(1  Ml  s  NSlZE-1 
JlsINOEX 

IP  <  J1-H1  1  29?,?V?,29H 

292  no  294  Ks1 . ISIZE 
t=(K-l )»MAXR 
LisL+INDEX 
HnLlla  A(L1  ) 
nn  294  j5ji,Mi. 

I=L  +  J 

294  A(  I  )  =  A(  1*1) 

294  A(l+1)sH0Ln 
298  l.iNSI  7E-M0nE4  I 

.1  =•  (MODF  -1  )  *  NC  •  MAXR  *1 

CALI  SW, (PX  (VECTOR! J), A,  H.  A(L),  VALUE.  MODE.  MSi 
1  INIIIX.  i:PSP.  NSIZE.  NITER(H0D0UT4  n  . 

CAI  I  OVt  Ri  I  (  lOVFLM  ) 

CAM  nvCHK  (  IDVnCT  ) 

OOTO  (  .lOVFI.  W 

Tim  IRR=IRR*1 

.Tii2  noTO  (  (u4,(itA)  .iiivncT 

.Tii4  IRR=IRR42 

TII4  IF  (  IRR  )  .11  n,  H  fl,62(» 

no  is(NC-i  )»NSI  ;r- 

no  .T|2  JS  INDEX.  NSIZE 

I  =  I  4j 

nilFSSd  )amiESS(L*  1  ) 

119  OUT SS( J)  =  nUESS(j4l  ) 

HSI/E  r  MSI/F-1 
N||»R(MOnE)  -  I  TFKSR4  ITEROR 
NAKSR(Mnl)|  )=  NSRAK 
NAKDRlHOnE )=  NORAK 

IF  (  ITERDR  )  .120/ 1611,. 160 

IPO  HonE=MnnF4i 

ITFROR.-I TFWnR 


.Tin,  n(i,42(» 


.”120  ,  1611 ,  .160 


MTR4"91 4 
HTR4091  t> 
MTR4II91  6 
MTR4n91 7 
MTR4II918 
MTR4()919 
HTR40920 
MTR4n922 
MTR4«924 
MTR40996 
MTR400?6 

HTR4099  / 
MTR40S9  ^ 
MTR409? V 
.  HAXHR,  |TERbR.MTR4u9-^^ 

)  HTR4(i9T, 

MTR4n94t) 

MTR4n946 

MTR4n94rt 

MTR4098() 

HTR409SI 

HTR41I982 

MTR4n98.J 

MTR4II9S4 

MTR4l)9h'j 

MTR40O*;/ 

HTR4r)98b 

HTR4n9‘>9 

)f  T  P  4 'm  6 1, 
PTR4  *  V*  \ 
»'TR4n94? 

HTR4I194.J 
MTR4n964 
MTR4II946 
PTR4n946 
HTR4n94/ 
HTR4n96H 
HTR40949 
MTR4U970 
MTR4U972 
7E.  MAXR.  NC,  HTR4IIV7/ 
(RR)  MTR4(l97b 

MTR40979 
MTR4098IJ 
MTR4II9H2 
MTR4()98v1 
M|R4n984 
Ml R4n9H6 
MTH4n9R4 
M'fR4ll9Hrt 
MTR409R9 
MTR40991) 
MTR40991 
MTR4n992 
MTR40994 
MTR4099t> 
MTR4fl996 
MTR4no97 
MTR4(I99H 
MTR4n999 
MTR41 nno 
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MiTf  R(Honr.-i  )=  (I. 

INfir  XsINDFX:' 

IF  (  INDFX-NSIZF  > 

IF  (  INlIfcXI-INnFX?  ) 
Tip  iNPIXs-JNOFX 
-INfil  XsINDRX-l 
nOTn 

^10  INPIXs-INDEX 

rOMTINIIF 
I  I  -  11 ♦ ISIZF 
Ii'=I?MSlZF 
I.1=I3»ISIZE 


HTR91001 
HTR41fin^ 
HTR41  »n.i 
MTR41 bP4 

kTR4l  nn*> 
HTR4inrt6 
HTR41 on/ 
HTR4’ 

M1R41 009 
><TR4ini  II 
HTR41 «1 1 
MTR41  01  ;> 


I 


NAKSRIMODF-I )  =  «. 

J*TR4T.  1  > 

NA|fllR(MnI)E-I  )  =  n. 

VT041  ;  !  4 

noTr»  PRO 

HTRo)  lit) 

I')  II 

n  =  n>isizF 

MTR41  111  / 

17sI1 ♦ISIZE 

MTR41 nj6 

I  i=i?*isize 

HTR41  111  V 

ooTn  1.10 

HTR41  (i?0 

4  iiri 

IF  (  K1-11  ) 

4in*4ii»4in 

HTR4IU2.? 

4 1  n 

no  41?  isi.Ksize 

MTR4lO?i 

J  =  K1*  1-1 

HTR4in?4 

JlS|1*|-t 

HTR410?*> 

ill 

H( J1 1sH( J) 

:iTR410?A 

4i:i 

no  414  Jsl.NC 

MTR41 0?/ 

1=  Nr*(M0llF-1  )♦  1 

MTR410?8 

JieIK*J-1 

HTR41(i?9 

4(4 

VAI.«IF(  I  )S  H(  Jl ) 

HTR41  iiJU 

(NKIXslNnFX 

MTR4in3i 

nOTO  790 

HTR410.1? 

5*10 

M0f)FsH0nF-1 

HTR4ln34 

IF  <  MOOR  ) 

5  l0,s.f|l,»>0? 

HTR41»35 

S(l? 

IF  (  NTAPUI  ) 

»;i  <<»**04 

HTR41036 

5 II 4 

RRAII  (  NTAPIIT  ) 

(A( I ). 5^1 .M) 

HTR4 1(1.3  7 

CMI  MUIT  (A.VbCrOP 

,  H,NSIZF;»HSIZF, 

HOnR>HAXK*NAXR>H«XR  .NCd) 

MTR41  (1.39 

no  *>06  Isl.HOiJF 

MTR41  (141 

J=  (1-1  )*I‘;!7R 

1 

NTR41  (14.) 

.11  =  (1-1  )*Nr»i 

HYR410-J4 

INOEXsO 

HTR41  (lAb 

ij  '1 A 

OALI.  NORM  (H(  J),H(.‘>,NSIZF,R«FSS< 

J:.  ).  INDEX.  HAXR  .'C.i  ) 

MTR4tb46 

•>(  n 

IF  (  NTAPOT  ) 

T  i!» 

MTR41 040 

A  1  7 

MRIIF  (NIAPnr.lii  ) 

MTR41  n‘»0 

no  *>7?  isi.Monr 

MTR41  (IS  > 

.tli:  HAXSR 

HTRAIOS.^ 

.17=  NITRR(  I  )-HAXSH 

PTR4lnS3 

f(  (  J?  ) 

S(4,*JlS»Sl'. 

HTR410h4 

S  1  4 

Jl  =NI TFR(  1  > 

NTRAIOS'.^ 

.C'sll 

HTR410S6 

no  TO  *>1/ 

MTR41(IS7 

•>  1  *. 

If  (  J7-MAXnH  5 

MTR41(IS8 

•>  1  A 

MPirF  (NTAPOT. 71) 

1 

HTR410S9 

*.  1  / 

no  TO  (S)«,*>7n)  , 

NO 

HTR41  (IA(I 

•>  1  n 

WPI  IF  (NTAPOT. 1  () 

1*  VALliRd).  J1.J7,  NAKSRr  I  ),NAKUR(  1  ) 

HTR41 OA^ 

nOTO  b/? 

HTR41(iaj 

•(  'II 

.1:  7»| 

HTR41  (IAS 

MirirF  (NTAPOT. ly> 

1.  VA(OR(.l-1). 

VALIIF(J).  JT.J?.  NAK<tK(l). 

MTR4106& 

1 

NAXDRT  1  ) 

HTR410A7 

o  7 

riiNT  INOI 

MTR410A8 

.(r  Nr:  •HOOF 

WKITF  (NIAI'OI  , /»'»  ) 

roil  MPRINI  (VtriOK.NSIZfc. J.MAXR»NTAPOI) 


MTR41  liAv 
HTR41{j7U 
MTR41071 
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‘>''4  WRflF  (NTAPni.?-?)  (  (6UFSS(  I  ),  1  =  1  ,  J)  ) 

*^AH  HPRIN1  (  H.NSIZF*  J.HAXR  .NTAPOT) 

RETURN 

4"n  NtTFP(HnDE)  =  I TFRSP* I TEROR 
NAKSR(HnnE)  =  NSRAK 
NARnR(MnDB)  =  norak 
IF  (  NTAPOT  )  ^iin  • 't(jO  ,Mii 

6»f1  .1=1 
M  =4 

nOTd  (AOP, OA  ) ,  IW9 
6')?  Ilrj 

noTfi  (Sfift 
^■14  t:< 

WRITF  (NTAPOT. ?/S)  (  A  T I  TIE  ( I  ) .  f  =  J.  j  i  ) 

OflTO  ft.50 

NtTFR(MnnE)  =  !  TFHSP*  I  TERi)R 
NAKSR(MODF)  =  NSRAK 
NAKPR{M006>  =  NORAK 

(  NTAPOT  »  S.I'UOji',  Ai  1 

.1=1 

,n  =  A 

GOTO  (M  4.  M  6 , 61  >1 ) ,  IHR 
6(4  Jl5X 

GOTO  618 

6 16 

618  WRITF  (NTAPOT. P6)  ( A T ITL E ( I  ) , | s J. Jl) ,  {ATlTLE(6)) 

(ATITLECJ?)),  (ATITLF<|), 1=8,9) 

.  ‘  »  6.1ft,  6.10, 6.>X 

6/>ft  .IPs/ 

OOTO  61? 

6??  WRITE  (NTAPOT, Ifl)  MODE 
IRRsft 

J=  (NonE-j  )*n(;*maxr 

no  6?4  |SI.ISI7F 

J'  *.)♦  f 

Jp5.J+(NC-1  )»MAyR«^I 
VTCTORI J?)=u. 

6/4  VETTORf  J1  )  =  ii . 

GOTO  .1’ii 
6.1ft  IsMOHE- 

WRMF  (N|.\P0T.I6) 

GOTO  '>1111 
FNft 

<  FORTRAN  I  S  roil.  OFCK 

n  SWFFPS 

G  .S  W  F  F  P  y 

n 

r  adHMirt  tmif  mini  ahh  smhfps  ii  from  thf  matrix,  ireai  r  cipmplfxi 

IIIRIII  .  IRIIF  HIIIIAI  nillllHNS,  AS  nilMPIIIFlI.  u  =  UTIIAMIC  HA’RIX 

SFRIFS  Of  MOI)M  Ihn  HORAL  COLUMNS. FL=  COLUMN  0^  MOLNVALUFt: 
SFRIFS  or  HOlllUFf.  MORAI  ROWS  OF  U.  ‘lOtNVALUFR. 

~  MOUF  now  PF.  INO  COMPUTFD.  N  =  SIZE 

DIMFNSIONFO  number  OF  ROWS  OF  U.  US  >  H,  HTRllF 
I  IF  PROrtIFM  IS  RFAL. 

?  IF  PROHI fM  IS  COHPLFX. 

SIIRIinilTINI  SRFEPX  (HTRUF.  11.  H.  IIS.  FI.  MI.IIE,  R.  pa.  r„, 

F.P,  MSI7F.  INRlS,  IRR) 


I 


S  U  B  R  0  U  T  I  N  E 


M 

US 

MOOf- 

Ml) 

NX 


I 


niMf  NS  I  ON 


—  -  '  !»»»'•  innz 

H(1),  US(l),  u(l).  HTRUKI).  FL(1).  0(4)  iRHISd) 


MTR41  (i7  i. 
MTR41  (1/ 
PTR41  (i7(i- 
MTR4  1  078 
HTR4  I  0  ;<■* 
HTR4|08(i, 
MTR4  1  ii>;  i 
H7R41  oa  X. 
M  I  R  4  1  t- 
HTR4  1 
M  r  R  4  !  '  s  >_ 

H  7  P  -1  w  j.  x““ 

M7R41 088 
P I  R  4  1  II H  C* 
MTR4  1  |i<vi.^^ 
MTR41  lifn*" 
PIR41 
MTR41  1)9.) . 
HTR41()94'* 
PTR41  l|0‘7 
MTR4  I  1(9!:'' 
PTR4  1  119/... 
MTR'*  1  II 9  8 
;,'IR41  (19  9- 
HIR4) IftOj 
HTR4)11)1“*- 
PTR4n  (iZ 

PIR411fi,i7 
PTR41  1  044. 
M|R41  )  Ilf, 

p  IR4 )  1  n  /-p 
MIR41  1  1  I)  I 
PTR41  n  ]  ■* 
MTR41  J  1 
PTR4n  1  .1  ^ 
MTR4  1114'/' 
HTR4  1 I  1 9 
PTR4 I n  ft  ' 
PTR'il  11  7  _ 
PT»4  1  1 1 8 

HTR4n  ?ij 
PTR41 l?i 
MTR41  l?i- 
PTR41 1 ?4 


P‘TR'1  t  1?ft 
PTR4n?/ 
PTR4) )?H 
MTR411?9 
MTR4n.1U 
MTR411.71 
MTR41 
PTR41 1 
hTR41  1  .(A 
PTR41  1.18 
MTR4)  l.<6 
MTR4H  ,7/ 
HTR41  1.78 

MT04 1 ^ \0 


m 
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tNn(  X=|ARS(  INItlX) 

MTO^I  Mo 

iNftiS(Moni  )  r  mnf-x 

. 

HTR411**/ 

H=Hnnfc-i 

HTR41MB 

Ki  =H»f4n#Hr) 

HTR91 15V 

NN=NC»HD 

MtR9ll6U 

K}?  =  mS|2F-M 

HTR41161 

fr  (  M  ) 

70. ;n,oA 

HTR4116^ 

/o  OOTO  (MU,  7?)  ,Nn 

HTR4tl6o 

i">  ir  (  Msi2E-Mn  ) 

74.14(1.140 

HTR41 164 

74  = 

>*194116* 

<stf  1  *Mn«N4l 

HTR41 16b 

nn  ih  1*1, N 

>•194116/ 

K=lf-1 

HTR411 60 

1  iL-1 

HTR41169 

M ( H ) =M( L  ) 

MTR41170 

V  H  (  L  )  =  f  . 

>•1941171 

«inTft  wn 

91941172 

HP  no  n'»  I  si  ,  H 

9TR4117.) 

.i1  *M0nF4  1 

HTR41174 

yf  (Nni55(  Ji)  =  fNni<;(  I ) 

••iRniiz** 

1i*n  IF  (  INDIX  ) 

1 0?,1 «4, 1 04 

91941176 

VI  K  V  s  K  ?  ♦  1 

9194117/ 

HsH-1 

91941170 

I'M  .11  sK1  ♦K/'*NC*  < 

91941179 

J^sM  ♦HS1ZE*N(: 

• 

9lR4linu 

IF  (  J1-J?  » 

1«S,UiS,Hl7 

9T94ll«l 

111‘i  no  in#i  isji.j? 

91941102 

1116  H(I)sn. 

9T94ll«d 

1(1 7  00  10  (1  14. 1  OHI.NC 

91941184 

1i)li>  JlsKI  1 

91R4nn> 

.|?sK1  ♦K2«NC*  1 

91941186 

no  1M  l*1.K7 

9194118/ 

jt  s.n  -1 

91R41188 

J^*J'7-1 

91941189 

110  H(.ll)sH(J?> 

91941 19U 

IF  (  M  ) 

1  1  H  ,  1  )  »> ,  1  1 1 

91941191 

1 1  no  11?  i  =  i.M 

91941 192 

H<  J? )  =  n , 

91R41 194 

11?  J?sJ9^1 

91941104 

1  1 4  IF  (  M  ) 

1  in,  1  in,i2«i 

91941  MtJ 

1  />n  11  =  1 

91941 19ft 

nn  1.10  isi.H 

9194',  19/ 

1 ,71  no  .J=  1  ,  M 

91941 19K 

ji=M0nf«.i 

91941199 

IF  (  II  •  INIIIS(JI)) 

91R41 ?no 

17?  noNiiNin 

91941201 

n  =  il*i 

91R41 202 

noio  Ml 

9IR4I  ?n.j 

17.1  iNniS(Ji);'» 

91941204 

n  (  iNniS(Hon«  )>i i 

)  1i'n,M4.M4 

9TR4l?llt> 

174  !Nnis(MonM  iNniS(Mnnf)»i 

91941206 

17».  MHSI7F-1I 

9194120/ 

II  (  1  ) 

17‘>,  MO, 176 

91941208 

176  .Ms  KMHSI7I 

9TR41 209 

.\v-  K1  *MSI/F  MNO 

-1  )»M|) 

91R41  21  It 

no  MH  .1-1.1 

91941 21  1 

.11  -  .11  -1 

9194121  <; 

.JVsJ?-! 

HIR41  21.1 

H(  J'7«  1  )  =  H(.17) 

91R4121 4 

1  7«  H(J1 ♦!  )  =  H( Jl  ) 

9TR412M 

i:70  .llsKMlI 

91R41 216 
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j?-j\ ♦(Nn-»  »»Hn 

MTR-1  1  21  / 

M( Ji  )^n. 

MTP41  21  «  :• 

H( J9)=n. 

MTR41 219  ( 

1  ,in 

rnNTiNUF 

HTH41 ??0  *' 

ir  <  INDIX  ) 

1  1  1 , 1  .s  <  ,  K11 

HTR41 ??l  .. 

«  IB 

HsH^I 

HTR4  1  v?-<f  \ 

i.ii 

no  1.19  1S1,M 

MTRIV??.) 

.M*MnnE*  1 

MTR41 224 

1  >9 

tNniS(  JUsINDISC  I  ) 

HTR412?*>:* 

\ 

1  Isl 

HTR41296L 

no  1.1R  Isl.H 

HTR*!!??/ 

! 

1 

no  11.14  J=1.M 

MTR41  ??H  -- 

i 

*  1 

JlsHOOE^.I 

HTR41 22v  i 

> 

IF  (  1  1-IMI)IS(  J1  )  )  1  1.^4»1K1'?,' 

154 

HTR41  ?;i.i  *' 

1154 

rONTINIIF 

HTR41P.1)  _ 

! 

- 

PTR412V 

\ 

noTo  1.15 

MTR41  2  5)- 

t 

1 1  1% 

l 1s{NC»HSI7» -1 1»MH+MnDF4l 

4TR4 1 254 

\ 

1  .^s(MS!/E-I  I  )»NM 

HTR412.1'5-- 

im»is(ji)  =  " 

HTR4121o;, 

j 

no  i;iB  .1=1.  NO 

MTR4121/''' 

.)9sl1-(  J-l  )»Hli-1 

1*10412  W.. 

1 

IF  (  L.'-NN  > 

15R»1  J4»1.14 

PTR41 2  19 . 

1  14 

no  i.i<i  i=i,l5,nn 

r 

HTR41  240  •* 

1 

J1SL1-I. 

HTRll 24  1 

jps J1 -NN 

\ ' 

PTR4  1242'* 

1 

1 

1 

1 

1.JA 

IJS<  Jt  )sUS(  J^1 

• 

1*1941  24  i.. 

1  IsLl-Hli 

#TR41244 

1  .18 

IIS(  J?)af. 

HTR4124t>.. 

14(1 

no  149  fsl.NN 

' 

HTR4t  ?4h; 

( 

.11  SKI  ♦! 

■ 

HTR4124/** 

J 

149 

HIRUF(  1  )sH(.|1  ) 

* 

HTR41 248 

5 

IF  (  H  > 

51 .51  ,« 

MTR41249r 

J 

5 

R 

no  ?•>  Isl.M 

HIR41 

KsMOnp-I 

HTR41251 

i 

t  |sNn*MlJ*(K-i  ) 

MIR41  2S2*- 

• 

\ 

TAIL  HULI  <US(  K 

).  HTRUF.  G.  1*  NSIZP.  1 

,  Hn.  f1<  .  1  .  NC.  1  ) 

HTR41 255 

$ 

If  (  0(1  )  ) 

12,9.  1? 

8TR41  259'  ' 

0 

on  in  (ti.ii'i.NC 

H  T  R  4  1  2  5  I'l .  . 

1 

m 

IF  (  0»  ’)  ) 

12.11.1? 

81941  ?(S  1 

} 

1 1 

|RR=IRR*'^ 

«T041  /•- 

19 

CONI  INIll 

•*  1  9  4  1  2  *  .' 

OOTn  ( 1 4 . |9 ) . Nr 

81  941  2f.4  ‘ 

1  4 

If  (amsifi  (M/n.(MniiF)-i.i»)  -  fp  i  i 

‘>.15.16 

MTR4 1 262 

l»» 

Or1.,n 

MIR41  266' 

J 

OOTn  1 / 

M1H4126/.. 

Ift 

o^in  (K>-f i.(Mi)n(  ))  / 

0 

MTR4  1  26b. 

1  / 

no  1H  J=1,MSI/F 

M1R41 269  • 

1  S  1. 1  ♦  .1 

HTR41 270 

l« 

MTHOr(J)'M{|. 1-0(1 

<«HTRUE( J) 

HTR4127i'  ■ 

; 

ncTO  9'^ 

MTH-ll  2/2 

1  <i 

Ks:*«K 

H1R41 275 

l  =  :'«HOnf 

8TR4  1  2  74. 

If  (  AHS((fl  (K-1  )*F|  (J-1  )»FL(K>*f  L(j))/(FL(.|- 

1  )»*?*Fl  (J. .*9)-l.U) 

8194127  5; 

1  -FP) 

81941276" 

/n 

If  <  Aonf  (FI  (K  >*ri.(./- 

1  )-Fl  (K-1  )«FL(  J) )  /  (FL(  .1-1  )##?.f  L(  .))••?)  1 

81R41  277., 

t 

1  -IP) 

91 ,?1 ,29 

81R41  27,i* 

J. 

/I 

0(1  )  r  1  .  {( 

0  (  P  )  .  II .  ti 

001  n  9.1 

8TR41  2  79-‘ 
81R4128U 
81R412H1|' 
MTR41282!. 

^•> 

(i(  .i  )  =  o(  1  ) *n(  / }  •*:■ 

USd-*) 


0(4)s(rL(K)-rL(J))»0(l>-<ft(K-n-FL(J-1  ))»6(?) 

R(  1  )- ( (f  L(  K-' )-Fl  (  J-1  )  )*0O  )  •  (FL<K  )*Fl.(  JM»0<2  > )  /  G(.^) 

R(?)s  0(4)  /  0(.)) 

M  no  ’4  J»  =  1.  .MSIZF 
k;?sji*mh 
L*L1  Ol 
I  ?s|  ♦MO 

0(3)sHTMUF(  JM 

HTRUF(J1>  =  H(t)  ♦  G('->»*HTR0F(K2)  -G(  l  )*HTR0E(  J1 ' 
HTRUF(Ki')-  H(LV>-  0(1  )*HTRU(:(K7)-G(?)*6{.n 
y4  rnuriNUF 
CONTIHUF 

t  sO 

fAl.l.  NORM  (  HTRflF  iHTRUE.  MSIZF#  6,  1»  MO,  MC,  1) 

n  =  I 

P  r  I 
L4  r  hoof 

no  4,i  isi.msizf 

I  1  =  J1 
I  V  =  J2 
I  .^sKWl 

no  33  HMsl.MOIlF 

IF  (  I  -INnlS(MM)  )  33,39.33 

»•»  rONTINUF 

PC  3/  Js’.MSIZF 
00  MMsI.MOOE 

IF  (  .l-INniS(MM<  '  3*>,37,3‘> 

3»i  CONTINIIR 

M(Ll)  *  Md?)  -  H  (-f'US(L4) 

GOTO  (3h.3«S),MC 
34  J3?  Ll»Mn 

J4*  l2^Hn 
J'is  L4*Mn 

J6S  1.3  ♦MI)  ... 

ll(J3)sU(.M)-H(l  3)«(IS(  J**)  -  H(J^)»US(L4> 

U(L1)=  U(M)  ♦  H( J4)»US( J‘>) 

3«  nONTlNUt 

L2sL?^J 
37  I3st3^1 

RfTO  41  ..... 
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SECTION  10 


NOMENCLATURE 


NOMENCLATURE 


a  Element  of  flexibility  matrix.  In. /lb 

Generalized  amplitude  coefficient  of  rigid-body  modal  series,  in.  or  rad 
b  Reference  semichord,  ft 

Element  of  oscillatory  aerodynamic  influence  coefficient  matrix,  dimensionless 
F  Control  point  force,  lb 

g  Structural  damping  coefficient,  dimensionless 

Control  point  deflection  due  to  rigid-body  motion,  in. 

Element  in  rigid-body  modal  matrix,  in.  or  dimensionless  (see  Section  II) 
Control  point  deflection,  in. 

K  Flexibility  matrix  normalizing  constant,  dimensionless 

Reference  reduced  frequency,  dimensionless 
M  Element  of  mass  matrix,  lb. 

M  Element  of  complex  mass  matrix  (includes  aerodynamic  effects) ,  lb 

2 

m  Element  of  generalized  mass  matrix,  lb.,  in. -lb,  or  lb-in  , 

m  Element  of  sum  of  generalized  mass  and  aerodynamic  matrices,  lb,  in, -lb, 
or  Ib-in^. 

2 

Q  Element  of  generalized  aerodynamic  force  matrix,  lb,  in. -lb,  or  lb-in  . 

R  Number  of  rigid-body  modes 

s  Reference  semispan,  ft  (i.e.,  span  measured  from  root  to  tip) 

U  Element  of  dynamic  matrix,  in, 

V  Velocity,  knots 

W  Element  of  aerodynamic  weighting  matrix,  dimensionless 


189 


SYMBOLS  (continued) 


X  Eigenvalue,  ^  in* 

3 

p  Atmospheric  density,  slugs/ ft 

£  Frequency,  cps 


Matrix  Notation 
^  j Square 
I  I  Column 
f  ^^Transposed 


r ij  Unit 
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